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ABSTRACT 


A  unified  procedure  is  employed  in  the  determination  of  the  penetration  resistance  of 
continuous  strip  footings  and  foundations  buried  up  to  ten  times  their  width  in  compact  soil. 
The  quasi-static  analysis  is  based  on  the  soil  failure  theory  of  classical  soil  mechanics. 

The  theoretical  data  is  presented  by  means  of  a  two  part  additive  equation,  two 
dimensionless  factors  describing  the  effect  of  weight  and  cohesion.  These  factors  are 
functions  not  only  of  the  angle  of  internal  friction  and  the  penetration  ratio  but  also  of  the 
shaft  parameters  and  the  dimensionless  group  dependent  on  the  density,  the  cohesion  and  the 
breadth  of  the  footing  The  factors  are  calculated  for  the  same  figure  of  rupture  and  the 
practice  of  treating  a  portion  of  weight  within  the  failure  boundary  as  a  surcharge  has  been 
eliminated.  The  inherent  approximations  and  resultant  errors  in  the  computation  are  fully 
staled. 


Experimental  verification  is  given  both  of  the  penetration  resistance  and  of  the 
failure  geometry  by  model  testing  in  a  glass  sided  soil  tank.  The  penetration  resistance  data 
in  compact  soil  support  the  quantitative  theoretical  values.  The  measured  failure  patterns  are 
shown  to  be  in  agreement  with  the  predicted  shapes  only  when  the  effect  of  changing  soil 
properties  and  constraints  are  carefully  evaluated  throughout  the  duration  of  the  test. 
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OBJECT 


This  programme  of  research  is  concerned  with  the  development  of  theory  of  penetration 
failure  which  may  be  used  as  the  basis  of  a  more  complex  theory  of  si ip-sinkage. 
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1.  INTRODUCTION 


The  evaluation  of  the  stress  deformation  characteristics  of  soil  not  only  forms  the 
basis  of  foundation  engineering  but  also  represent  one  of  the  major  problems  hindering  progress 
in  soil-vehicle  relationships.  In  the  past,  this  common  interest  in  one  and  the  same  problem 
has  not  generated  viable  co-operation  between  the  two  fields.  In  fact,  rather  the  reverse  is 
true,  partly  due  to  the  divergence  in  practical  requirements. 

For  the  design  of  foundations,  complete  soil  failure  at  its  ultimate  bearing  capacity 
must  be  avoided  with  an  adequate  margin  of  safety,  and  the  total  and  relative  settlements  must 
be  less  than  the  limits  tolerated  by  the  superstructure.  The  point  of  complete  failure  is 
determined  by  the  maximum  stress  in  soil,  whereas  the  settlement  is  due  to  soil  consolidation 
under  the  weight  of  the  structure,  and  depends  on  smaller  stresses  which  can  be  described  by 
elastic  analysis.  Thus,  bearing  capacity  theories  are  based  on  mathematical  analyses  for 
ideal  soils,  and  minor  deviations  from  realistic  conditions  are  adequately  safeguarded  by  the 
adoption  of  a  factor  of  safety  of  three. 

In  mobility  investigations,  the  theories  of  penetration  resistance  must  be  more  precise 
for  application  in  the  more  complex  theories  of  motion  resistance,  slip  sinkage  and  the 
ultimate  immobilisation  of  the  vehicle.  Greater  emphasis  has  been  placed  on  empirical  rela¬ 
tions  and  similarity  techniques  to  overcome  the  problems  of  describing  accurately  the 
influence  of  stress  strain  characteristics,  compressibility  and  layered  conditions  in  the 
disturbed  top  soil.  Although  it  is  doubtful  whether  these  factors  will  ever  be  predicted 
accurately  in  a  rigorous  analysis,  any  empirical  analysis  must  be  firmly  based  on  scl1 
mechanics  theory.  The  solution  must  be  capable  of  modification  and  extension  to  include  the 
influence  of  combined  loading  and  vehicle  attitude  on  the  mode  of  failure  and  the  load  bearing 
of  the  soil. 

1 . 1  Analytical  Techniques 

The  diverse  approaches  employed  to  determine  the  stress  deformation  characteristics 
of  the  soil  may  he  divided  into  the  following  five  categories:- 

(1)  Static  analysis,  a  rigorous  mathematical  analysis  based  on  ,he  stress 
conditions  at  failure  in  a  rigid  heavy  Coulomb  material. 

(2)  Plasticity  theory ,  the  association  of  static  analysis  with  suitable  stress  strain 
relationships. 

(3)  Soil  mechanics  theory,  an  approach  based  on  1  or  2  above,  but  involving 
simplifying  assumptions  in  order  to  obtain  reasonably  economical  solutions  to  real 
problems. 

(4)  A  semi-empirical  analysis  employing  measured  soil  stress  deformation 
characteristics. 

.(5)  Model  experiments  guided  by  dimensional  analysis. 

The  main  difficulty  with  the  various  types  of  rigorous  analysis  — for  example  static  analysis, 
plasticity  theory  and  limit  analysis — is  that  they  are  not  relevant  to  real  soils,  in  order  to 
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obtain  immediate  solutions  to  practical  problems,  soil  mechanics  theory  has  incorporated  much 
of  this  more  fundamental  work  without  fully  defining  all  the  assumptions  and  approximations. 
The  result  is  that  particular  problems  have  numerous  solutions  which  have  little  point  of 
contact  with  each  other.  Both  rigorous  and  theoretical  analysis  are  combined  with  empiricism 
and  dimensional  theory  to  yield  a  wealth  of  possibilities.  The  potential  and  limitations  of 
these  evaluation  techniques  are  illustrated  in  this  review  of  previous  work. 

1.2  Static  Analysis 


Perhaps  the  most  relevant  approach  to  soil  mechanics  is  through  the  method  of  static 
analysis,  particularly  associated  with  the  work  of  Sokolovski1.  It  is  based  on  the  properties 
of  a  rigid  Coulomb  material  and  combines  the  known  relationships  between  the  stresses  at 
failure  with  the  partial  differential  equations  of  equilibrium. 

The  simplest  way  of  describing  stress  conditions  at  failure  is  by  the  Coulomb 
equation2, 


r  =  c  +  a  tan  0 


1.2.1. 


where  r  represents  the  maximum  shear  stress  along  any  plane  of  incipient  failure,  c  the 
cohesion,  a  the  normal  stress  on  the  plane  and  0  the  angle  of  internal  shearing  resistance. 
A  more  general  version  of  this  yield  condition  showing  the  relation  between  the  normal  and 
shear  stresses  on  any  two  mutually  perpendicular  planes  at  a  point  of  incipient  failure  is> 
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The  equations  of  static  equilibrium  are:- 

8  °x  +  fjxy  = 

fix  fiy 
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1.2.3. 


The  total  distribution  of  stress  throughout  a  mass  of  failing  soil  can  be  obtained  by 
integration  of  these  equations  starting  from  known  boundary  stresses.  This  can  be  done  by  a 
step  by  step  numerical  method  using  more  convenient  versions  of  equations  1.2.2.  and  1.2.3. 
Sokolovski  has  published  an  account  of  the  applications  of  this  method  to  many  rather  simple 
problems.  Prager 3  developed  an  alternative  graphical  technique  called  the  Pole  Trail  method 
which  was  used  by  Josselin  de  Jong4  to  obtain  the  bearing  capacity  of  a  surface  footing  on 
a  frictional  material  with  weight  and  surcharge. 

The  most  fundamental  proposition  that  emerges  from  this  analysis  is  that  slip  planes 
occur  in  pairs  inclined  at  an  angle  of  n/2  +  0  to  each  other  at  all  points  in  a  mass  of  soil  at 
failure.  This  can  be  seen  immediately  from  the  Mohr  diagram  in  Fig  1.2.1. 
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In  genera),  the  slip  line  field  will  consist  of  two  sets  of  curved  lines  intersecting  at 
the  correct  angle  of  -n/2  +  (?.  Where  a  discontinuity  occurs  in  the  boundary  stress  conditions 
(as  at  the  corner  of  a  footing)  two  different  slip  line  fields  meet  each  other.  In  order  to  obtain 
a  continuous  slip  line  field,  a  third  fan  shaped  zone  is  introduced  with  one  set  of  slip  lines 
emanating  from  the  point  of  discontinuity  or  singular  point  (Fig  1.2.2.) 

When  a  uniform  loading  i.s  applied  parallel  or  at  right  angles  to  the  gravitional  force, 
the  curved  slip  lines  are  stra’ghtr  t.ed  to  form  stress  fields  in  the  well  known  acti  ve  and 
passive  Rankine  stales,  respectively.  In  addition,  if  a  singular  point  occurs  between  two 
straight  slip  line  fields,  then  all  the  lines  emanating  from  thm  point  are  straight  radii,  and  the 
other  family  are  inclined  at  an  angle  0  to  the  circumferential  direction  to  satisfy  the  condition 
of  equilibrium  so  that  the  slip  lines  are  at  an  angle  of  v  /2  +  0.  The  increment  of  length  of 
the  curved  slip  line  isl¬ 


and  after  integration 


r  d  f  sec  0  -  dr  cosec  0 
r.  =  r^  e  6  tan  0 


.  1.2.4. 


where  rc  is  the  radius  for  i  -  0  (Fig  1.2.3.).  Thus  the  circular  failure  plane  is  replaced  by 
a  logarithmic  spiral  when  friction  is  present  in  the  radial  fan. 


Static  analysis  can  only  be  applied  to  problems  where  the  soil  mass  at  a  state  of 
failure  is  completely  isolated  by  boundaries  of  known  Stresses  and  known  failure  surfaces 
within  the  soil.  Examples  are  the  surface  foundation  and  retaining  wall  problems  illustrated 
in  bigs  1.2. ‘la.  and  1 lb.  Note  that,  in  these  cases,  ail  the*.. soil  within  the  boundaries 
ABCDF  and  A13CDFFGI1  is  definitely  faded. 

Unfortunately,  the  techniques  cannot  yield  a  complete  solution  to  the  deep  foundation 
problem  show n  in  big  1.2. a.  I  his  is  because  a  large  part  o!  the  zone  Mi) EG  consists  of  soil 
which  is  definitely  not  in  a  state  of  failure.  It  is  possible  to  obtain  the  zone  ABC  by  static 
analysis  and  to  see  that  there  must  be  a  fan  with  initially  straight  radii  emanating  from  B.  The 
other  limits  of  the  fan  MOD  would  normally  be  obtained  from  the  slip  line  field  in  the  region 
BDFG  by  integration  inwards  from  known  boundary  stresses  along  DC  and  CL.  This  cannot  be 
done  because  the  soil  is  not  in  a  state  of  failure  along  these  boundaries. 

A  rigorous  solution  to  this  problem  would  require  a  knowledge  of  the  stress  strain 
relation  forthe  soil  extending  over  the  whole  elastic  plastic  range.  This  would  provide  a 
separate  set  of  kinematic  constraints  from  which  extra  equations  could  be  established.  Such  a 
system  cf  static  and  kinematic  requirements  constitutes  a  theory  of  plasticity. 

1.3.  Theories  of  Plasticity 

The  combination  of  static  requirements  based  on  force  equilibrium  and  a  failure 
criterion,  and  kinematic  requirements  based  on  stress  strain  relations  is  the  basic  feature  of 
any  plastic  theory.  The  most  successful  solutions  have  been  developed  for  incompress;blt, 
frictioniess  materials  falling  according  to  the  Tresca  enterion:- 


T 


=  k 


.  1.3.1. 
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Normally  (e.g.  for  metals)  K  is  far  greater  than  the  gravitional  stresses  which  may  then  be 
neglected,  resulting  in  the  assumption  of  a  weightless  material.  The  way  in  which  stress 
strain  relation  can  be  described  in  a  plastic  material  is  highly  complex.  Let  it  suffice  to  say 
that  it  is  derived  from  the  concept  of  a  yield  function  which,  for  a  stable  material,  requires 
that  the  strain  rale  vector  is  always  normal  to  the  yield  surface  when  axes  of  plastic  strain 
are  superimposed  on  the  axes  of  stress  in  Mohr’s  circle  construction.  The  ratio  of  plastic 
extension  to  plastic  shear  rate  is  uniquely  determined  by  the  shape  of  the  yield  surface. 

This  theory  of  plasticity  has  been  used  with  limited  success  for  certain  simple 
problems  concerning  the  working  of  metals  or  of  relatively  strong  clays.  When  more  complex 
boundary  conditions  are  introduced,  the  complicated  mathematical  expressions  for  the 
curvature  of  the  slip  lines  are  not  always  capable  of  yielding  a  general  solution. 

Even  in  a  perfectly  plastic  weightless  material  with  0  =  0,  few  problems  have  been 
solved.  Prandtl5  achieved  an  exact  solution  for  the  identation  of  a  flat  surface  by  a  long  flat 
punch  (Fig  1.3.1a.).  The  regions  ABC  and  BDE  are  regions  of  constant  stress  whilst  BCD 
is  a  zone  of  radial  shear.  However,  Hill6  proposed  an  alternative  stress  distribution  for  the 
same  problem  which  gives  the  same  bearing  capacity  (Fig  1.3.1b.).  The  velocity  fields  for 
the  two  solutions  were  added  by  Shield  7 .  In  Prandtl’s  solution,  the  rupture  distance  BE  is 
twice  that  in  Mill’s  solution,  but  the  velocities  are  halved  so  that  the  volume  of  material 
raised  above  the  undisturbed  surface  after  the  punch  has  penetrated  a  small  distance  is  the 
same  in  both  cases.  It  has  been  proposed  that  the  double  wedge  is  applicable  to  a  smooth 
base  and  the  single  wedge  to  a  rough  base  when  friction  is  present.  This  distinction  is  later 
used  in  soil  mechanics. 

Considering  the  penetration  of  the  load  into  an  ideally  plastic  material,  the  surface 
bearing  capacity  is  not  the  maximum  base  penetration  resistance,  and  it  increases  with  depth 
to  some  steady  value.  At  this  ultimate  bearing  state,  Jaky8  has  postulated  the  slip  line 
field  for  the  case  of  plain  strain  (Fig  1.3.2.).  In  the  zone  of  radial  shear  BCD,  the  slip  lines 
turn  through  an  angle  of  zr radians  instead  of  7r/2  radians  in  Prandtl’s  solution.  The  two 
regions  of  constant  stress  ABC  and  BDE  remain  intact,  but  the  rupture  surface  CDE  breaks 
out  at  the  vertical  free  surface  as  opposed  to  the  horizontal  free  surface  in  the  surface  case. 

In  view  of  the  difficulty  of  obtaining  a  precise  solution,  an  approximate  method  known 
as  limit  analysis  has  been  developed  within  this  theory  of  plasticity.  The  collapse  or  limit 
design  theorems  were  first  introduced  by  Drucker  and  Prager9and  refer  to  the  constant  loads 
if  the  accompanying  changes  in  geometry  are  disregarded.  In  other  words,  the  equilibrium 
conditions  may  be  set  up  in  the  undeformed  body  and  represent  a  great  simplification  because 
the  deformation  constitutes  one  of  the  unknowns  in  the  collapse  problem. 

The  two  theorems  are':- 

(a)  Collapse  will  not  occur  if  any  state  of  stress  can  be  found  which  satisfies  the 
equations  of  equilibrium  and  stress  boundary  conditions. 

(b)  .  Collapse  must  occur  if,  for  any  compatible  flow  pattern  considered  as  plastic 
only,  the  rate  at  which  the  external  forces  do  work  on  the  body  equals  or  exceeds  the 
rate  of  internal  dissipation. 
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Based  on  the  first  theorem,  a  perfectly  plastic  material  will  adjust  itself  to  carry  the 
load  if  there  is  a  possible  way.  This  statically  admissible  state  of  stress  gives  a  lower 
bound  solution  for  the  equilibrium  stress  conditions  at  every  point  of  a  plastic  zone.  From  the 
second  theorem,  the  soil  will  fail  if  a  path  of  failure  exists.  This  kinematically  admissable 
velocity  field  gives  an  upper  bound  solution  which  is  only  valid  for  perfectly  elastic  materials. 
Thus  a  lower  bound  solution  is  maximised  whilst  an  upper  bound  solution  is  minimised.  The 
actual  solution  lies  at  the  point  where  these  two  approaches  have  the  same  value  for  the 
applied  force.  The  use  of  these  theorems  is  tedious  and  there  is  no  guarantee  that  the  upper 
and  lower  bound  solution  will  converge  accurately  on  the  same  mean  value. 

Attempts  have  been  made  to  develop  a  theory  of  plasticity  for  a  material  with  weight 
and  internal  friction,  which  would  be  appropriate  to  real  soils.  The  Tresca  yield  criterion  has 
been  replaced  by  the  Von  Mises  -  Coulomb  criterion  and  the  material  assumed  to  be  non-work 
hardening.  This  activity  is  well  represented  by  the  recent  work  of  Drucker10, 
Haythornthwaite11  and  Shield12.  This  attempt  has  failed,  however;  the  resulting  material  has 
the  property  of  enormous  dilation  with  shearing  strain  and  is,  therefore,  quite  unlike  soil13. 

One  negative  outcome  of  the  work  is  the  realization  that  the  limit  analysis  theorems  need 
extensive  modification  before  they  can  be  applied  to  frictional  materials14. 

Frequent  reference  to  the  theory  of  plasticity  is  made  in  published  work  in  soil 
mechanics.  Terzaghi’s  work15  on  bearing  capacity  is  based  on  Prandtl’s  solution  and 
Meyerhof’s16  on  Jaky's.  In  the  particular  case  of  strong  frictionless  clays,  there  is  some 
justification  for  this,  due  allowance  being  made  for  the  controversial  nature  of  the  plasticity 
solutions  themselves.  However,  the  minimisation  procedure  which  are  often  used  in  soil 
mechanics  and  justified  by  vague  references  to  limit  analysis  do  not  appear  to  be  correct. 

1.4  Shallow  Foundation  Theory 

Soil  mechanics  commenced  as  a  serious  theoretical  study  with  the  publication  of 
‘Theoretical  Soil  Mechanics’  by  Terzaghi  in  1943  1S.  The  basis  for  the  theories,  developed 
therein,  is  clearly  the  theory  of  plasticity  for  frictionless  weightless  materials.  Modifications 
to  suit  heavy  frictional  soils  have  been  attempted  by  numerous  workers,  each  applying 
different  assumptions  and  simplifications.  No  major  effort  is  made  to  explain  why  rigorous 
solutions  of  real  problems  are  not  possible,  nor  to  develop  suitable  principles  for  non-rigorous 
methods.  This  unsatisfactory  position  in  current  soil  mechanics  can  best  be  appreciated  by  a 
critical  examination  of  the  existing  theories  of  bearing  capacity. 

The  failure  geometry  assumed  by  Terzaghi  for  a  surface  foundation  is  shown  in 
Fig  1.4.1a.  From  the  conditions  of  stress  beneath  the  base  of  the  footing,  he  postulated  the 
existence  of  a  wedge  shaped  body  of  soil  ABC,  which  remains  permanently  in  a  state  of 
elastic  equilibrium  and  behaves  as  part  of  the  sinking  footing.  Furthermore,  the  shape  of  the 
wedge  has  dependent  on  the  base  roughness  such  that  the  base  angle  0  could  take  any  value 
between  ^/4  +  0/2  for  a  perfectly  smooth  base  and  0  for  a  rough  base.  The  latter  value  was 
adopted  for  two  reasons,  first  because  he  found  that  the  angle  of  base  friction  was  always  less 
than  the  angle  of  shearing  resistance,  and  secondly  because  he  assumed  that  the  footing  could 
sink  only  if  the  soil  immediately  below  the  wedge  tip  C  moved  vertically  downwards.  For  this 
condition,  the  stress  trajectories  must  emanate  from  a  vertical  tangent.  This  gives  rise  to  an 
impossible  stress  condition  at  point  C,  where  four  failure  planes  intersect  in  a  Coulomb 
material  at  angles  of  other  than  the  correct  values  of  7r/2  +  O'  and  7r/2  -  0  radians. 
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For  a  surface  foundation,  the  region  BPK  is  uninfluenced  by  the  structure,  ami  is 
dependent  only  on  the  weight  of  the  soil  below  the  horizontal  soil  surface  plus  any  uniform 
surcharge  layer,  1  Ivse  known  stress  conditions  yield  a  plane  shear  zone  in  the  passive 
llankine  .'.(ate,  '1  he  maximum  extent  of  this  region  is  dotei mined  by  the  outer  boundary  of  the 
fan  shap'd  zone  which  emanates  from  the  lip  of  the  wedge  beneath  the  fooling. 

The  calculation  of  the  resultant  bearing  pressure  was  derived  from  a  two  pa:t  calcula¬ 
tion  and  presented  in  terms  of  an  additive  equation. - 

p  cNc  +  qNq  +  MyliN-y  . 1.4.1. 

where  Nc  ,  N,.  and  Ny  are  dimensionless  factors  d'-p.-nding  only  on  0.  The  effect  of  weight 
wits  obi. lined  for  a  purely  frictional  material  and  that  of  cohesion  and  surcharge  for  a  weight¬ 
less  m  penal  w  ith,  Irietion.  1  he  former  is  reasonable  since  it  is  the  condition  foriliv  sand, 

I""  the  .second  is  completely  unrealistic  and  was  chosen  only  because  it  greatly  simplified 
llii'  compnln li on  of  the  bearing  capacity. 

Based  on  an  extension  of  l'randll’s  ana  lysis  s,  Tei/.agiii  assumed  thai  the  pole  of  the 
logarithmic  spiral  is  at  point  A  for  the  caicul.it ion  of  the  cohesion  and  surcharge  factors.  This 
is  entirely  sal  i  d  notary  since  it  gives  the.  correct  orientation  of  stress  alone  face  AC,  and 
is  tangential  to  the  spiral  CD’.  For  calculating  the  weight  factor,  T'T/.aghi  adopted  a 
semi-graphical  treatment  based  on  the  work  of  Ohdo  1 '.  The  pole  of  the  spiral  is  restrained  to 
move  anywhere  along  BD  or  that  line  produced.  The  actual  point  is  chosen  such  that  the 
lesul taut  bearing  capacity  is  a  minimum.  I  nh  s  tin  j.oie  of  the  spiral  is  ai  lie.-  (oruei  of  the 
footing,  the  stress  direction;,  are  inmn-eei  vie  re  tney  ian-rsecl  the  weJg<  boundary  ! it ;  ami 
tin1  failure  surfaces  CIl’lF  are  dissimilar. 

In  a  weightless  media,  the  fixed  .soil  body  and  the  passive  llankine  zone  are  se|xi rated 
by  a  zone  of  radial  shear,  comprising  radii  and  logarithmic  spirals  whole  pole  lies  at  the 
corner  of  the  fooling.  Ohdo  promised  that  these  constraints  were  excessive  for  a  heavy  Coulomb 
material  when  existence  of  a  doubly  curved  slip  line  field  was  known.  For  matheniatiea i 
convenience,  he  retained  the  concept  of  a  zone  of  radial  shear,  but  relaxed  lli"  requit'-mcm  that 
the  region  must  radiate  from  the  corner  of  the  fooling.  By  restraining  t tie-  pob-  to  the  boundary 
of  the  passive  llankine  zone,  only  one  stress  discontinuity  occurred.  This  discontinuity  was 
present  along  the  retaining  wall  in  his  theory  and  was  ignored  because  the  obliquity  of  the  slip 
planes  was  influenced  by  the  friction  oil  the  interface.  Thus,  both  sources  of  error  were 
collected  at  one  point  and  their  effect  minimised  by  the  adoption  of  a  minimum  value  |iostulate, 
Since  the  minimising  procedure  was  not  curried  out  for  a  series  of  failure  surfaces  which  all 
satisfy  the  conditions  of  equilibrium,  thri  boundary  c  onditions  and  the  idealised  soil  properties, 
there  seems  little  logical  justification  for  using  it.  There  is  even  less  reason  for  its  application 
in  the  case  of  bearing  capacity  where  the  stress  conditions  along  the  soil  to  soil  boundary  of 
the  fixed  soil  wedge  have  already  been  positively  assumed.  If  it  is  initially  assumed  that  the 
fan  shaped  zone  is  enclosed  by  two  straight  boundaries  (the  sides  ot  ihc  wedge  and  the  passive 
llankine  zone)  meeting  at  a  singular  point,  then  it  follows  implic  itly  that  the  intermediate  slip 
line  field  comprises  radii  and  logarithmic  .spirals  with  their  pole  ,,t  tig  singular  point  i.e,  at  the 
corner  of  the  footing. 

I  he  weak  link  in  this  reasoning  lies  in  the  initial  assumption  of  the  stress  conditions 
which  exist  beneath  the  footing.  A  study  of  the  effect  of  weight  and  surcharge  on  the  ultimate 
bearing  capacity  of  a  cohesionless  material  was  first  conducted  by  l.undgren  and  Morlensen  1 
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This  analysis  was  based  on  plasticity  theory.  The  failure  patterns  which  result  from  the  three 
combinations  of  weight  and  surcharge  are  shown  in  Fig  1.4.2.  The  upper  diagram  illustrates 
the  traditional  Prandtl  failure  pattern  for  a  weightless  soil  with  surcharge.  The  inclusion  of 
soil  weight  in  the  analysis  results  in  a  curved  wedge  boundary  and  a  doubly  curved  slip  line 
field  in  the  fan.  When  both  weight  and  surcharge  are  considered,  the  failure  surface  is  unique 
and  it  has  been  postulated  that  the  ultimate  bearing  pressure  in  a  cohesionless  soil  can  be 
expressed  more  correctly  by  the  formula:- 

P  =  M  (qNqC0  +  '47BKy0)  . 1.4.2. 

where  Nq«j  is  the  surcharge  factor  for  a  weightless  soil,  Nyo  is  the  weight  factor  for  zero 
surcharge  and  g  is  a  function  of  0  and  the  surcharge  to  weight  ratio.  For  0  -  30°  the 
correcting  factor  p  is  given  in  Fig  1.4.3.  As  its  maximum  value  is  only  1.17,  the  error  made 
when  neglecting  p  is  thus  not  excessive. 

It  is  also  interesting  to  note  that  the  value  of  u  remains  positive  for  quite  large 
surcharges.  In  Terzaghi’s  theory  where  shallow  sinkage  is  considered  in  terms  of  a  surcharge 
soil  layer,  the  maximum  value  of  the  correcting  factor  when  0  =  30°  occurs  for  a  penetration 
ratio  of  0.4  and  the  correcting  factor  is  greater  than  zero  provided  the  penetration  ratio  is  less 
than  ten.  This  lends  proof  to  Terzaghi’s  original  argument  that  by  neglecting  the  shear 
strength  of  the  overburden  and  discarding  a  correcting  factor,  the  computation  is  simplified 
and  any  error  is  located  on  the  safe  side.  Unfortunately,  Lundgren’s  value  of  Nyo  corresponds 
to  r  solution  which  is  statically  correct  but  kinematically  i.  y  s'ible  and  is  now  thought  to  be 
over  conservative. 

An  alternative  solution  was  developed  by  Gorbunov-Passadov1  9  who  based  his 
analysis  on  the  assumption  that  the  velocities  in  the  plastic  zones  have  the  same  direction 
as  the  active  slip  surfaces.  Although  this  assumption  is  open  to  criticism,  he  did  provide 
some  experimental  substantiation  from  a  photographic  analysis  of  the  particle  paths  at  failure. 
From  his  calculations,  there  results  a  two  part  wedge  (Fig  1.4.4.).  The  inner  elastic  zone 
ABC  is  almost  the  same  shape  ns  that  found  by  l.undgren.  From  the  outer  boundary  of  the, 
compacted  core  ALU),  the  slip  lines  curve  to  make  vertical  contact  with  the  boundary  of  the 
inner  elastic  core.  With  this  type  of  failure  mechanism,  the  value  of  Ny  is  much  larger  than 
those  given  by  the  classical  method  (F  ig  2.10.1.).  It  is,  however,  a  valuable  attempt  e,  in, a 
a  statically  correct  solution  based  on  the  kinematic  evidence  obtained  from  tost  data. 

The  two  preceding  theories  l8’  19  suffer  from  this  disadvantage  that  cohesion  is  not 
considered.  A  less  vigorous  but  more  general  analysis  has  been  proposed  by  1 1 u 20  in  his 
Variable  Factor-Theory’.  1  he  slip  surface  comprises  a  logarithmic  spiral  with  its  pole 
constrained  at  the  corner  of  the  footing  and  a  straight  sided  wedge  under  the  foundation  whose 
variable  base  angle  0  is  determined  by  the  condition  for  which  the  ultimate  bearing  capacity 
is  a  minimum  (big  1.4.5.).  Hence,  the  minimum  value  postulate  depending  on  the  position  of 
the  spiral  pole  is  replaced  by  one  depending  on  the  base  angle  of  the  wedge  using  the  same 
moment  equation  as  the  basis  of  the  analysis.  It  therefore  suffers  from  the  same  criticism  as 
Terzaglii’s  small  wedge,  namely,  that  the  failure  planes  at  the  wedge  tip  do  not  meet  and  the 
correct  angles  in  a^oulomb  material.  It  will  later  be  shown  that  the  minimum  value  postulate 
wholly  depends  on  the  mathematical  constraints  of  the  problem  when  dissociated  from  the 
correct  stress  conditions  and  not  on  the  mechanism  of  failure  (c.f.  Section  5). 

Howe's  solution21  combines  the  procedure  of  slices  as  described  by  Bishop  22  and 
Janbu23  with  the  stress  dilatancy  relations.  In  the  failure  geometry,  the  slip  line  is  replaced 
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by  the  a  -line  where  a  is  the  average  inclination  of  the  major  plane  at  which  the  particles  are 
interlocked  (Fig  1.4.6.).  The  angle  of  rotation  of  the  principle  stresses  when  measured  from 
the  vertical  is  denoted  by  in  the  figure  and  when  measured  from  the  horizontal  by  Cp.  This 
analysis  has  several  advantages.  The  stress  boundaries  are  nowhere  violated,  the  analysis  can 
be  simplified  by  approximate  methods  without  involving  errors  of  practical  importance,  and  the 
solutions  are  identical  with  those  based  on  the  Mohr-Coulomb  criterion  when  the  dilatancy  is 
zero. 

The  solution  seems  to  provide  a  suitable  compromise  or  correlation  between  the  two 
■conflicting  groups  adopting  the  straight  and  curved  wedge  boundaries.. 

1.5  'Deep  Foundation  Theory 

Terzaghi  applied  the  term  ‘deep  foundation’  to  one  whose  sinkage  is  greater  than  the 
width.  The  first  comprehensive  study  was  conducted  by  Meyerhof16.  He  extended  Terzaghi's 
bearing  capacity  theory  to  deal  with  deep  foundations,  and  made  a  major  contribution  in  that 
the  failure  geometry  was  not  only  dependent  on  0  but  also  on  the  sinkage.  A  failure  pattern 
was  proposed  with  a  smooth  transition  from  the  surface  case  to  that  of  a  deep  foundation 
(Fig  1.5.1.).  The  zone  of  mobilised  soil  is  again  divided  into  three  regions  -  the  wedge,  the 
plane  shear  zone  and  the  zone  of  radial  shear.  As  sinkage  increases  the  segment  if  the  spiral 
is  extended,  and  the  plane  shear  zone  is  reduced  in  size. 

Meyerhof  assumed  that  the  wedge  angle  must  always  have  a  base  angle  of  tt/4  + 

0/2  regardless  of  the  roughness  of  the  footing,  but  in  accordance  with  placticity  theory,  stated 
that  beneath  a  perfectly  smooth  base,  two  wedges  were  formed  instead  of  one  which  effectively 
halves  the  total  bearing  capacity  (c.f.  Figs  1.3.1a.  and  1.3.1b.).  In  practice,  however,  a  smooth 
base  is  never  found,  and  may  be  discounted. 

For  the  determination  of  the  failure  geometry  and  the  calculation  of  the  bearing  capacity 
factors,  Meyerhof  was  obliged  to  assume  a  certain  pressure  distribution  along  BC,  El)  and  DB 
(Fig  1.5.1a.).  Along  the  wedge  boundary  BC,  he  assumed  that  the  weight-frictional  component 
of  pressure  would  always  have  its  centre  of  pressure  two  thirds  of  the  way  from  B  to  C.  From 
the  analogy  of  hydrostatic  pressure,  this  would  be  correct  for  a  surface  foundation,  but  as  the 
sinkage  increased,  the  centre  of  pressure  would  move  towards  the  centre  of  BC. 

The  resultant  forces  acting  on  the  foundation  shaft  BG,  and  the  weight  of  the  adjacent 
soil  wedge  BEG  are  replaced  by  'equivalent  stresses’  normal  and  tangential  to  the  plane  BE, 
which  may  then  be  treated  as  an  ‘equivalent  free  surface’.  It  is  not  state.d  however  that  these 
stresses  are  uniformly  distributed  along  BE,  nor  is  it  made  clear  that  the  pressure  due  to  both 
cohesion  and  weight  are  equally  dL.  .buted  on  all  sides  of  the  triangle  BDE.  The  startling 
assumption  comes  from  using  Mohr’s  circle  to  determine  the  equilibrium  of  the  triangle  BDE 
upon  which  the  geometry  depends.  The  Mohr  circle,  however,  is  a  construction  that  is  correct 
for  the  stresses  at  a  point  where  uniform  distribution  of  stress  is  a  reasonable  assumption. 

The  bearing  pressure  was  represented  by  the  equalion;- 

P  =  cNc  +  pcNq  +  VtrBNy  . 1.5.1. 

where  p0  is  the  normal  equivalent  free  surface  stress  along  BE.  The  calculation  of  the 
cohesion  and  surcharge  factors  is  again  based  on  an  extension  of  Prandtl’s  work  in  a  weight¬ 
less  material,  and  the  pole  of  the  radial  shear  zone  lies  at  the  corner  of  the  footing.  For  the 
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calculation  of  the  weight  factor,  Meyerhof  minimised  with  even  greater  abandonment  than 
Terzaghi  and  has  attempted  to  justify  this  in  published  work24  on  surface  foundations.  Even 
further  discontinuities  are  introduced  by  the  complete  relaxation  of  the  restraints  on  the 
spiral  centre  (Fig  1.5.1.).  Unfortunately  the  unrestrained  pole  necessitates  discontinuous 
slip  planes  along  the  L  jndaries  of  the  wedge  BC,  and  the  passive  Ranking  zone,  BI).  It  is 
also  unacceptable  on  the  grounds  that  the  failure  surface  CDE  has  a  sharp  comer  in  it  at 
point  D  which  is  unlikely  where  soil  movement  is  occurring.  Unless  these  discrepancies  are 
small,  the/e  remains  little  support  for  describing  the  shape  of  the  curve  by  the  equation  of  the 
logarithmic  spiral,  apart  from  the  advantage  that  of  the  frictional  forces  pass  through  the  pole. 
Compared  with  Terzaghi’s  N'y  factor  for  the  same  wedge  angle  of  77/4  +  0/2  he  did  obtain 
values  of  up  to  10/c  lower  at  the  larger  values  of  0,  but  only  by  sacrificing  the  more  accept¬ 
able  stress  field.  With  the  inclusion  of  the  sinkage  parameter,  the  freedom  of  movement  of 
the  pole  is  much  more  considerable  and  impossible  stress  field  can  be  derived.  In  the  presence 
of  the  surcharge,  zone  BEG  is  also  neglected  so  that  the  geometry  reduces  to  that  of  a  footing 
on  a  slope.  Meyerhof  also  shows  that  there  is  little  difference  in  the  value  of  N  y  if  the  plane 
shear  zone  is  replaced  by  an  extension  of  the  radial  shear  zone,  but  rotational  equilibrium  of 
the  forces  was  not  examined  in  either  case.  It  is  hardly  surprising,  therelore,  that  the  failure 
geometry  for  the  N  y  calculation  is  much  smaller  titan  that  for  the  weightless  material;  when 
the  plane  shear  zone  is  absent,  the  former  may  Le  half  the  other. 

In  the  generalisation  procedure,  Meyerhof  confuses  the  concept  of  a  surcharge  with 
that  of  a  lorce  due  to  soil  weight  and  cohesion.  He  considers  triangle  BEG  as  a  surcharge. 
There  seems  no  justification  for  this  and  it  is  perfectly  possible  to  include  the  effect  of  this 
triangle  in  the  weight  term.  This  leaves  open  the  possibility  of  considering  any  surcharge 
applied  to  the  surface  EG. 


In  fact,  Meyerhof  did  attempt  to  express  the  resultant  bearing  capacity  by  the  eqimtion:- 


P  =  ci\  +  ‘Ay  BN 

cq  yq 


.  1.5.2. 


where  one  term  represents  the  influence  of  cohesion  and  the  other  the  influence  of  weigh;. 

The  factors  Ncq  (depending  on  Nc  and  Nq)  and  Nq  (depending  on  Nyand  \  J  are  functions 
of  0,  the  penetration  ratio  and  the  soil-structure  properties  along  the  shaft.  Equation  1.5.2. 
gives  the  base  resistance  of  the  foundation  and  any  skin  friction  along  the  shaft  must  be 
added  to  obtain  the  totai  bearing  capacity.  The  equation  may  be  used  for  only  two  special 
cases.  The  composite  factors  both  depend  on  Nq  which  in  turn  is  governed  by  the  interaction 
of  cohesion  and  weight.  For  a  general  solution,  it  is  only  possible  to  determine  Nc_  and 
Nyq  when  one  or  other  is  zero  or  may  be  neglected;  Ncq  can  Le  found  for  a  purely  cohesive 
material  when  0  -  0,  and  Nyq  for  a  purely  frictional  soil.  No  solution  is  possible  fora 
frictional  soil  with  cohesion  unless  some  arbitrary  division  of  N'q  is  assumed.  In  other  words, 
the  dependence  of  the  bearing  capacity  factors  on  the  cohesion  to  weight  ratio  is  hidden  in 
the  calculation. 


Balia25  considered  the  same  problem  of  a  deep  foundation  but  recognised  that  the 
bearing  factors  depended  on  the  cohesion  to  weight  ratio  in  addition  to  0  and  the  penetration 
ratio.  The  same  sliding  surface  was  used  to  calculate  all  the  N  factors  and  is  shown  in 
big  1.5.2.  lie  observed  that  lerzaghi’s  proposal  for  two  curved  failure  surfaces  to  be 
parallel  and  coincident  w*ith  each  other  at  point  C  is  impossible,  and  opens  them  out  to  the 
correct  angle.  Unfortunately  he  retains  Terzaghi’s  small  wedge  which  still  leave  angle  ABC 
incorrect,  and  an  intersection  of  four  failure  planes  at  point  C  which  is  not  possible  for  a  rigid 
Coulomb  Material. 
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According  to  Balia,  the  failure  surface  CDE  is  made  up  of  a  circle  CD  and  a  straight 
line  DE.  It  is  not  possible  for  a  rigid  body  bounded  by  a  circle  and  a  straight  line  to  move 
without  opening  up  large  cavities  along  the  boundary  CDE.  Neither  is  it  possible  for  the 
circular  segment  OCD  to  rotate  pushing  the  triangle  ODE  upwards  without  cavities  appearing 
along  the  line  OD. 


Another  serious  criticism  is  that  the  theory  has  been  restricted  to  a  maximum  penetra¬ 
tion  ratio  of  1.5,  without  any  possibility  of  an  extension  to  greater  depths.  This  limitation 
has  arisen  from  the  assumption  that  the  plane  boundary  DE  must  always  make  an  angle  of 
77/4  -  0/2  with  the  horizontal  free  surface.  This  is  a  valid  statement  in  itself  but  it  has  two 
repercussions:  first,  that  the  increasing  effect  with  depth  of  the  foundation  side  thrust  on  the 
base  resistance  cannot  be  evaluated  and  secondly,  that  it  allows  for  no  minor  re-orienlation 
of  the  failuie  surface.  For  example,  an  inflexion  in  the  curved  boundary,  or  a  slight  angular 
displacement  of  the  plane  shear  zone  boundary  suggested  by  Meyerhof16  would  satisfy  this 
condition  without  giving  additional  complexity  to  the  basic  theory.  On  Balia's  rigid 
interpretation  of  this  condition,  the  rupture  distance  must  increase  rapidly  with  depth  and 
ultimately  tend  to  infinity. 


The  anal  ysis  is  based  on  Kotter’s  equations  of  static  equilibrium  and  the  moment 
referring  to  the  centre  of  the  circular  arc  is  also  zero.  Even  though  a  surcharge  zone  is  not 
delineated,  it  is  noted  that  the  final  solution  is  presented  in  the  usual  form  of  a  three  part 
additive  equation,  instead  of  combining  all  the  forces  due  to  weight  in  a  single  N  factor 
(Eqn  1.4.1.). 


The  N  factors  according  to  this  theory  are  markedly  different  from  all  previous  values. 
The  Ny  factor  is  much  greater  and  therefore  the  breadth  of  the  foundation  has  more  influence 
on  the  value  of  the  bearing  capacity  (Fig  2.10.1).  The  increase  due  to  depth,  however,  is  much 
less  than  in  Meyerhof’s  theory,  but  the  latter  is  hardly  a  proper  basis  for  comparison  because 
of  the  anomalies  previously  reviewed.  Balia  also  stated  that  the  theory  is  best  applied  to 
granular  soils  or  to  those  with  little  cohesion,  so  that  N'c  is  not  of  great  significance. 

1.6  Empirical  Factors  in  the  Basic  Equation 

The  theoretical  analysis  may  be  applied  only  to  the  idealised  situation  of  a  centrally 
loaded  continuous  footing  in  an  idealised  material.  For  more  practical  problems  in  real  soils, 
many  other  effects  must  be  considered  by  empirical  chances  to  the  basic  equation. 

The  influence  of  soil  compressibility  and  the  state  of  drainage  presents  an  added 
complication  to  in-situ  testing.  In  a  compressible  material,  the  practical  application  of  any 
theory  calculated  for  soils  at  their  minimum  density  will  load  to  sources  of  error.  Tcrzaghi15 
overcame  this  problem  by  suggesting  an  empirical  reduction  in  tan  0  to  two  thirds  of  its 
maximum  value,  whereas  Meyerhof16  introduced  a  compressibility  factor  which  varied  from 
tan  0  for  shallow  footings  to  0.85  tan  0  for  deep  footings.  Similar  factors  weie  evolved  for  the 
influence  of  ground  water  conditions  but  are  not  of  relevance  to  the  present  study. 

This  practice  was  extended  to  evaluate  other  modifying  influences  such  as  shape  and 
depth  of  the  foundation  and  inclination  or  eccentricity  of  the  load.  Denoting  the  shape  factor 
s,  the  depth  factor  by  d,  and  the  inclination  factor  by  i,  Brinch  Hansen  2  6  generalised 
Terzaghi  s  formula  by  multiplying  each  term  with  a  set  of  factors:- 


cNcscdcic 


qNqSqdqjq 


JdTBNySydyiy 


.  1.6.1. 
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Using  the  known  relationship  between  Terzaghi's  Nc  and  Nq  factors,  the  number  of 
parameters  may  be  reduced  by  rearranging  equation  1.6.1.  in  the  form:- 

p  =  (c  +  q  tan  0)  Ncscdcic  +  'A  yBNySyd^i  y  .  1.6.2. 

However,  this  method  of  overcoming  the  limitations  of  a  basic  theory,  which  itself 
contains  many  simplifications,  can  yield  some  anomalous  results.  For  example,  d-y  must 
always  by  equal  to  unity  because  N-.  is  calculated  for  c  =  O  in  the  basic  equation,  and  the 
inclinations  of  the  load  would  tend  to  cause  assymmetrical  failure  instead  of  the  symmetrical 
failure  assumed  in  the  original  theory.  Thus,  both  depth  and  inclination  factors  must  bear  no 
relation  to  the  original  theory  to  be  of  practical  significance. 

On  the  other  hand,  it  is  realistic  to  adopt  an  empirical  shape  factor  to  modify  the 
bearing  capacity  of  a  continuous  strip  footing  to  that  of  a  rectangle  on  the  basis  of  similarity 
of  the  soil  movement.  The  factor  is  then  only  accounting  for  the  end  effect;  the  technique 
becomes  increasingly  suspect  however,  as  the  strip  approaches  a  square  or  circle.  Conversely, 
the  aspect  ratio  L/B  is  important  in  the  selection  of  the  length  cf  footing  to  represent  Iwo 
dimensional  failure  in  model  studies.  Skempton27  investigated  the  influence  of  shape'  for  the 
special  case  of  clay  and  proposed  t h a t : - 

Nc  (rectangle)  =  1  +  0.2  ^ 

Meyerhof16  presented  his  experimental  data  in  the  form  of  a  graph  and  suggested  that  the 
combined  shape  factor  is  a  function  of  density,  0  and  the  penetration  ratio.  Briiicii  Hansen26 
combined  these  results  with  Terzaghi’s  and  formulated  further  shape  factors  independent  of 
depth.  For  small  aspect  ratios,  cohesion  and  surcharge  become  increasingly  more  important 
than  the  weight  of  the  material 

Nc  (rectangle)  -  1  +  (0.2  +  tan6  0)'|“  Xc  (strip) 

Ny  (rectangle)  -  I  -0.5  (0.2  +  tan^'p"  Ny  (strip) 

There  is  ample  experimental  evidence  to  corroborate  these  empirical  relations  for 
shape  and  they  are  later  used  to  determine  the  aspect  ratio  of  the  footings  used  in  the  current 
tests. 

1  -7  Semi  Lfmpiricai  Analysis 

The  problems  involved  in  soil  failure  beneath  vehicles  are  similar  in  nature  to  those 
of  civil  engineering  earthworks  but  are  generally  more  complex.  All  the  early  attempts  to  solve 
them  have  been  based  on  describing  the  soil  by  measuring  its  response  to  certain  simple  load¬ 
ing  tests.  hmpirical  curve  fitting  equations  were  used  to  describe  the  experimental  results, 
and  these  used  as  a  basis  of  mathematical  analysis  of  the  more  complicated  vehicle  situations. 

The  general  pressure  sinkage  equation,  which  forms  the  basis  of  the  current  empirical 
approach  is  attributed  to  Bernstein28  and  Letoshnev29  and  is  of  the  form:- 

P  =  . 1.7.1. 


Nc  (strip) 
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where  p  is  the  pressure  beneath  the  footing.  The  values  of  the  sinkage  modulus  k  and  the 
exponent  n  were  obtained  by  plotting  experimental  pressure  sinkage  results  on  logarithmic 
axes.  The  exponent  was  found  to  be  fairly  independent  of  the  soil  properties,  but  it  was  later 
found  that  the  sinkage  modulus  was  a  function  of  the  size  and  shape  of  the  footing. 

In  an  attempt  to  improve  Bernstein's  equation  to  cover  the  complete  soil  range, 
Bekker30  proposed  that  the  relationship  may  be  expressed  by 


where  kc  and  k0  are  soil  constants  which  may  be  determined  front  ^ests  for  two  different  plate 
sizes.  Although  kcis  predominant  in  cohesive  soils  and  k0in  frictional  soils,  they  are  not 
true  soil  constants  and  were  still  found  to  vary  considerably  with  plate. size  and  shape. 


A  more  recent  change  to  the  pressure  sinkage  relationship  has  been  suggested  by 
Reece31.  Influenced  by  Meyerhof’s  bearing  capacity  theory,  he  proposed  the  equation  in  the 
form:- 


P 


(c  k'c  +  Vty B  k>0) 


. 1.7.3. 


which  is  dimensionally  more  attractive  and  introduces  the  various  parameters  in  a  more 
acceptable  way.  The  relationships  between  Meyerhof’s  bearing  capacity  factors  (equation 
1.5.2.)  and  the  constants  in  the  equation  above  are  given  by> 


Equation  1.7.4.  links  the  empirical  approach  with  theoretical  soil  mechanics  and  it 
should  be  possible  to  calculate  the  values  of  k'0  and  k'c  for  compact  soil'.  This  is,  in  fact, 
the  w-ay  in  which  this  project  originated.  Even  for  loose  soils,  an  equation  like  1.7.4.  which 
is  based  on  some  theoretical  reasoning  should  be  superior  to  others  and  a  comprehensive 
series  of  measurements  made  by  Wills32  has  endorsed  this. 


1.8  Dimensional  Analysis 


Dimensional  Analysis  has  also  been  used  to  describe  the  phenomenon  of  penetration 
failure  by  means  of  a  dimensionally  correct  equation  containing  selected  variables.  Although 
a  number  of  variables  in  the  problem  may  be  reduced  by  dimensional  reasoning,  r-<~  ’  er  the 
complete  solution  nor  the  inner  mechanism  is  revealed. 

The  most  important  variables  in  the  problem  of  penetration  failure  may  be  summarised 
by  the  following  equation 


p  =  f  (0,  y,  c,  v,  d,  z)  . 1.8.1 . 

where  v  denotes  the  velocity  of  penetration  and  d  the  characteristic  dimension  of  the  probe. 
With  the  exception  of  velocity  which  may  be  neglected33,  these  parameters  have  been  combined 
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to  form  almost  every  conceivable  dimensional  group  and  inverse  function,  some  of  which  bear 
little  relation  to  the  parallel  work  in  soil  mechanics. 


In  a  more  realistic  study  of  model  analysis  for  earthmoving  equipment  Emori  and 
Schuring3*  employed  the  following  functional  equation  to  describe  soil  deforming  processes 
in  general:- 

p/yd2  =  /((?,  c/yd,  vVgdi  .  1.8.2. 


For  the  specific  case  of  penetration  failure  in  a  composite  soil  with  both  cohesion  and 
friction,  it  was  tentatively  suggested  that  equation  1.8.2.  could  be  written  in  the  form:- 

p/yB2  -  c/7B.  fj(0)  +  /2(0)  .  L8'3- 

where  the  characteristic  dimension  B  denotes  the  breadth  of  the  fooling  and  /,  ,  and  /,  are 
functions  of  0.  As  the  comparison  of  various  footing  sizes  is  possible  only  when  there  is 
both  geometric  similarity  and  a  consistent  penetration  ratio,  it  follows  that  the  functions  /] 
and  / 2  are  also  dependent  on  the  aspect  ratio,  L/B,  and  the  penetration  ratio. 

The  final  form  of  equation  1.8.3.  is  identical  to  the  current  civil  engineering 
approach,  but  it  is  no  nearer  the  ultimate  solution;  the  main  contribution  of  the  technique  lies 
in  its  ability  to  present  the  c  omplexities  of  the  problem  in  an  organised  manner. 

1.9  Experimental  Measurements 

The  experimental  verification  of  any  failure  theory  in  real  soils  is  extremely  difficult. 
The  measurement  of  the  soil  parameters  requires  a  high  degree  of  accuracy  and  considerable 
skill.  The  soil  structure  variables,  and  the  variation  of  the  soil  parameters  with  bulk  density 
must  also  be  determined.  Osman35  summarised  the  results  obtained  by  translational, 
compressive  and  torsional  methods  showing  that  identical  values  of  cohesion  and  0  can  be 
measured  independently  by  different  methods  of  testing. 

In  an  attempt  to  verify  the  actual  .shape  of  the  failure  zone  beneath  a  foundation,  many 
workers  19-  36*39  have  employed  a  photographic  technique  through  a  glass  sided  tank.  No 
reference  in  published  literature  has  been  found  on  the  effect  of  soil  interface  fret  ion  and 
t Be  validity  of  any  photographic  analysis  is  open  to  question.  This  problem  may  be  overcome 
by  testing  the  model  in  the  centre  of  a  tank  full  of  sand  mixed  with  metallic  particles.  An 
X-ray  technique  can  then  be  adopted,  but  the  complexity  is  a  very  real  deterrent.  Alternatively 
the  sand  may  be  replaced  by  cylindrical  metal  or  wooden  rollers  of  varying  diameters37.  This, 
however,  diverges  from  the  physical  reality  of  a  natural  soil.  The  over-simplified  two 
dimensional  mechanism  of  failure  exhibited  by  the  roller  model  can  not  yield  sufficient 
corroborative  evidence  on  its  own  merits  alone  but  it  does  provide  a  simple  tool  for  demonstra 
tion  purposes  during  the  development  of  the  theory.  For  the  additional  tests,  the  most 
convenient  solution  is  to  carefully  investigate  and  minimise  the  effect  of  interface  friction 
prior  to  analysing  the  experimental  data*0. 

Both  the  force  measurements  and  the  failure  geometry  depend  on  plate  size  and  shape. 
The  influence  of  aspect  ratio  in  the  study  of  tw  o  dimensional  failure  has  already  been 
stressed.  Not  only  is  aspect  ratio  important,  but  also  the  breadth  of  the  footing.’  Dimensional 
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reasoning  con  no  longer  be  applied,  especially  in  compressible  soils,  when  the  plate  breadth 
is  small  in  relation  to  particle  size  because  lateral  compression  occurs  rather  than  general 
shear  failure.  It  has  been  suggested  that  the  minimum  breadtfy  lies  between  1  and  l!;  in  and 
will  increase  with  compressibility31.  As  this  estimate  of  the  minimum  breadth  was  obtained 
in  a  soil  which  was  below  the  maximum  bulk  density,  it  is  probably  on  the  safe  side. 
Nevertheless,  any  results  obtained  for  footings  of  less  than  1  inch  wide  must  remain  suspect. 

Tip  available  pressure  penetration  test  data  for  high  aspect  ratio  footings  are 
surprisingly  limited  compared  with  the  vast  fund  of  theoretical  work.  The  major  part  of  the 
available  experimental  work  has  been  carried  out  in  a  dry  sand  which  is  a  convenient  material 
to  use,  or  in  soils  with  little  cohesion,  lit  his  paper,  Balia25  has  summarised  the  experi¬ 
mental  data  for  the  surface  bearing  capacity  of  foundations  and  presented  a  comparison 
between  theoretical  and  test  results  in  the  form  of  a  graph.  Selig  and  McKee41  conducted 
static  loading  tests  on  small  footings  resting  on  the  soil  surface.  Meyerhof16  also  presented 
detailed  experimental  verification  oi  this  theory  for  dorp  sinkage  over  the  complete  soil  range. 
More  recently,  test  results  have  been  published  by  Mults42,  Uebeque39  and  Biarez37,  the 
latter  including  skin  friction  data.  However,  in  many  deep  loading  trials  with  strip  founda¬ 
tions,  the  effect  of  skin  friction  along  the  shaft  is  not  considered  separately.  Although  much 
more  research  on  skin  resistance  is  available  for  piies,  the  effect  of  the  small  aspect  ratio  is 
unknown.  Unless  the  exact  nature  of  all  the  test  conditions  are  fully  described,  it  is  not 
always  possible  to  apply  the  results. 

A  similar  situation  exists  in  the  field  cf  soil  vehicle  mechanics,  and  test  data  for 
footings  of  reasonable  size  and  shape  in  uniform  constant  soil  conditions  are  not  readily 
available.  Only  one  test  in  compact  sand  by  Wills 3 2  may  be  compared  with  the  present  theory 
for  incompressible  soils. 

Finally,  the  prediction  of  the  penetration  resistance  by  means  of  a  quasi-static 
analysis  depends  on  negligible  velocity  effect.  In  a  study  of  this  phenomenon,  Vesic,  Banks 
and  Woodward43  have  demonstrated  that  the  minimum  bearing  capacity  of  a  surface  footing  in 
both  dry  and  saturated  sand  occurs  at  a  penetration  velocity  of  0.25  in  min  and  is  approxi¬ 
mately  30rf  less  than  that  for  static  loading.  Thereafter,  the  bearing  capacity  increases  again 
at  a  slower  rate  with  speed.  Haley  and  llegedus33,  found  that  the  bearing  capacity  in  saturated 
sand  at  a  penetration  speed  of  1750  in  min  was  up  to  30\  higher  than  that  at  0.6  in  min.  The 
apparent  drop  in  bearing  capacity  at  low  sjieed  is  attributed  to  an  increase  in  compressibility 
when  particles  have  insufficient  time  to  reorientate  after  each  new  load  increment.  As  speed 
increases  the  loading  is  too  rapid  for  shear  failure  to  occur  along  the  paths  of  least  resistance 
and  punching  shear  replaces  general  shear  failure.  At  very  high  loading  rates  in  saturated 
sand,  apparent  cohesion  is  also  developed  by  the  negative  pore  water  pressures.  Similar 
velocity  effects  were  obtained  by  Goodman44  from  small  scale  footing  tests  on  clay.  It  may 
be  concluded,  therefore,  that  the  effect  of  velocity  is  not  insignificant.  However,  in  the 
vehicle  situation  where  high  rates  of  loading  occur,  it  is  justifiable  to  assume  that  the  penetra¬ 
tion  resistance  is  equivalent  to  the  static  bearing  capacity  by  neglecting  the  trough  in  the 
curve  of  load  versus  speed. 
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2.  THEORY 


Soil  Mechanics  cannot  be  a  rigid  discipline  and  cannot  be  treated  as  an  exercise  in 
mathematics  because  the  mathematical  properties  of  soil  are  so  complex.  An  adequate  solu¬ 
tion  to  any  problem  can  be  achieved  only  by  a  scientific  comtination  of  theory  and  experiment. 
It  has  been  proposed  that  nothing  less  than  the  following  constitutes  an  adequate  Lasis  for  a 
solution. 

Fibst,  the  properties  of  the  idealised  soil  for  the  theoretical  analysis  must  be  fully 
stated.  This  must  be  followed  by  a  careful  description  of  both  the  behaviour  of  such  soil  and 
the  simplifying  assumptions  which  are  necessary  to  overcome  mathematical  difficulties.  These 
soil  properties  and  simplifying  assumptions  must  be  combined  with  tbc  normal  conditions  of 
equilibrium  to  determine  the  shape  of  the  failure  surface  within  the  soil.  Ideally  every  part  of 
the  soil  involved  in  the  failure  should  be  in  equilibrium.  Ibis  ir.av  prove  impossible,  but  the 
departure  from  the  correct  equilibrium  conditions  must  Le  small  and  must  be  described 
quantitatively.  Furthermore-,  the  correct  stress  conditions  must  exist  at  every  point  on  the 
boundaries.  Once  again,  this  inuv  not  be  possible,  bu'  divergencies  from  the  ideal  must  be 
small,  they  must  be  shown  to  cause  only  small  changes  in  the  final  answers  and,  ideally,  they 
should  be  shown  to  represent  an  approximation  to  a  more  accurate  but  more  complex  failure 
pattern.  The  failure  pattern  so  developed  must  be  verhed  by  carefully  controlled  experiments 
on  a  limited  number  of  soils  which  cover  a  large  pan  of  the  possible  soil  range.  Once  the 
failure  pattern  has  been  established,  the  forces  involved  can  be  calculated  and  a  suitable 
computer  programme  developed.  The  results  of  these  force  calculations  must  be  experimentally 
verified  over  the  same  wide  range  of  soil  types.  Finally,  generalising  equations  must  be 
developed  to  enable  the  communication  of  the  results  oi  ihe  computer  programme  in  a  simple 
way.  The  errors  involved  in  this  generalisation  must  also  be  slated. 

The  whole  of  this  process  of  interlocking  theoretical  and  experimental  methods  is 
necessary  because  any  sir. fcle  part  is  subject  to  such  large  possibilities  of  error.  It  is  difficult 
even  to  measure  the  basic  soil  parameters  as  they  exist  in  the  problem  under  investigation, 
far  less  to  apply  them  in  a  theoretical  analysis,  and  it  is  difficult  to  repeat  experiments 
maintaining  the  same  sod  conditions.  The  process  described  above  would  require  such  a  total 
corroboration  of  theory  and  experiment  that  it  would  leave  little  room  for  doubt  that  not  only 
were  the  right  answers  achieved,  but  that  they  were  achieved  for  the  right  reasons.  ■ 

2.1.  The  Basis  of  Soil  Failure  Theory 

The  soil  is  heavy,  tfnd  it  is  assumed  to  be  homogeneous,  rigid  and  incompressible  and 
not  to  dilate  when  shearing.  It  fails  according  10  Coulomb’s  equation  from  which  it  results 
that  lines  of  incipient  failure  occur  in  pairs  making  an  angle  of  tt/2  +  0  to  one  another  in  a 
plane  strain  problem.  This  is  an  idealisation  m  that  ?al  soil  is  not  quite  incompressible  or 
rigid  and  does  dilate  slightly  as  it  shears.  It  will  be  shown  later  than  an  initially  homogeneous 
soil  can  become  completely  non-uniform  during  loading. 

\Mien  the  soil  slides  aiong  the  surface  of  the  loading  structure,  the  stresses  at  the  soil 
structure  interface  are  assumed  to  be  described  by  a  modification  of  Coulomb’s  equation.  In 
this,  the  cohesion,  c,  is  replaced  by  the  adhesion,  ca,  and  the  angle  of  internal  shearing 
resistance,  0,  by  the  angle  of  soil-structure  friction,  i.  Under  these  conditions,  the  soil- 
structure  interface  is  not  a  failure  plane,  and  if  an  interface  exists,  the  correct  orientation  of 
the  slip  planes  must  be  detennined  from  Mohr’s  circle. 
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Th1-  above  failure  criterion  together  with  known  boundary  stress  conditions  are 
sufficient  to  make  certain  rather  simple  problems  statically  determinate,  and  a  rigorous  solu¬ 
tion  may  be  obtained  by  Static  Analysis.  Frequently,  this  technique  can  be  employed  in 
special  cases  which  form  the  basis  of  a  more  general  theory,  but  when  the  stress  boundary 
conditions  cannot  be  rigorously  determined  over  the  entire  failure  surface,  the  complete  solu¬ 
tion  of  many  practical  problems  requires  the  adoption  of  further  kinematic  assumptions. 

In  most  soil  structure  problems,  the  point  of  failure  is  represented  by  the  loading  at 
which  flow  begins.  There  are  certain  characteristics  of  soil  flow  which  are  common  to  the  flow 
of  any  ^<.her  incompressible  material.  The  motion  will  arrange  itself  so  that  the  energy,  lost 
in  the  friction  generated  by  the  soil  flow,  is  minimised.  This  requires  that  the  paths  of 
individual  particles  are  as  straight  as  possible  and  only  deviate  from  straight  lines  by  making 
the  smallest  angles  of  de.flection  on  the  largest  radii.  Soil  will  not  flow  along  paths  with  sharp 
corners  except  on  a  plane  of  symmetry. 

The  kinematics  of  soil  flow  in  response  to  a  small  movement  of  the  structure  can  be 
simplified  by  defining  only  two  possible  types  of  flow.  In  the  first,  large  masses  of  soil  move 
as  rigid  bode  .  It  is  not  possible  for  soil  flow  to  result  in  the  appearance  of  large  cavities 
within  the  soil;  this  is  a  common  observation  which  depends  on  the  low  strength  to  weight 
properties  of  soil.  If  the  soil  does  not  move  as  a  rigid  body,  then  it  is  convenient  to  describe 
the  other  type  of  movement  as  shearing  flow.  Zones  of  shearing  flow  will  generally  be  bounded 
by  non-circular  curves  and  when  the  soil  exhibits  high  strength  to  weight  properties,  these 
boundaries  need  not  be  sharply  defined.  A  given  displacement  car.  be  achieved  with  a  minimum 
energy  input  by  rigid  body  motion.  This  is  therefore  the  preferred  mode  of  movement,  but  may 
not  be  possible  duo  to  the  boundary  conditions  of  the  problem. 

The  static  analysis  of  rigid  Coulomb  materials  leads  to  the  concept  of  continuous  slip 
line  fields.  In  practice,  those  are  only  continuous  in  zones  of  shearing  flow.  The  slip  line 
fields  are  non-existent  over  the  main  area  in  which  rigid  body  motion  taken  place.  The  absence 
of  com  1;o  line  fields  and  the  occurrence  of  rigid  body  motion  is  made  possible  because 

the  sc  wly  compressible.  Definite  deformations  along  the  slip  surfaces  are  necessary 

in  ordi  top  slip  lines  and  soil  failure. 

The  forces  involved  in  the  shearing  of  soil  have  two  components,  one  due  to  the 
cohesion  and  internal  friction  of  the  soil,  and  the  other  due  to  the  weight  and  the  internal 
friction.  In  practice,  these  forces  are  generated  together  simultaneously  along  the  same 
failure  plane  w  hose  path  is  dependent  on  the  angle  of  internal  shearing  resistance.  The 
relative  magnitude  01  these  two  components  can  vary,  and  the  solution  should  be  correct  when 
either  is  zero.  It  seems  logically  justifiable,  therefore,  to  divide  the  force  on  the  structure 
into  two  parts,  one  part  due  to  the  cohesive-frictional  soil  forces,  and  the  other  due  to  the 
weight-frictional  forces.  This  simplifies  the  analysis. 

The  equilibrium  of  the  soil  masses  at  failure  depends  on  the  distribution  of  the  normal 
and  shearing  stresses  on  the  straight  and  circular  failure  surfaces  of  the  rigid  body  zones. 

(This  does  not  include  the  zones  of  shearing  flow).  This  distribution  is  difficult  to  obtain  by 
rigorous  methods.  It  was  originally  proposed  to  assume  that  the  stresses  due  to  chohesive- 
fr ict; ona i  forces  were  distributed  uniformly  along  the  straight  or  circular  failure  surfaces  and 
further  to  assume  that  the  stresses  due  to  the  weight-frict  jn  1  forces  were  distributed  in  a 
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hydrostatic  manner,  namely,  that  they  were  proportional  to  the  depth,  zo,  from  the  nearest 
part  of  the  free  soil  surface  vertically  above  the  failure  plane.  These  two  assumptions  are 
described  by  the  following  relationship: 


a0c  '  70c  *  k0c 

<70  y  dZ  T0  y  £2  ^0y  Z0 


.  2.1.1. 


Along  both  the  failure  surface  boundaries  of  any  single  rigid  body  zone,  the  two  co¬ 
efficients  k0c  and  k0y  were  assumed  to  have  constant  values  which  satisfied  the  condition 
that  the  values  must  be  identical  at  the  intersection  of  the  two  boundaries.  The  determination 
of  the  two  coefficients  from  the  horizontal  ami  vertical  equilibrium  of  the  plane  shear  zone  gave 
the  correct  relationships  for  the  active  and  passive  Rankine  zones.  With  the  introduction  of 
sinkage,  an  unbalanced  couple  yielded  a  source  of  error. 


In  the  following  analysis  the  relationships  are  modified  by  considering  rotational 
equilibrium  of  the  plane  shear  zone  and  include  a  dimenionless  exponential  term  such  that: 


T0c  T0c  '*■'  0c  “H  — 


,  /z  \  n 

70y  <C  T0y  CC  K0y  <C  L \\ 


.  2.1.2. 


where  -1  <  n  <  3.  This  effectively  permits  a  variable  pressure  distribution  to  be  used  along  the 
failure  planes  to  obtain  complete  equilibrium.  This  does  not  affect  the  limiting  case  of  the 
passive  Rankine  zone  or  the  condition  that  the  values  of  the  coefficients  for  ihe  two  planes 
must  be  identical  at  the  point  of  intersection.  It  should  be  noted  that  these  relationships  are 
purely  empirical.  The  basis  for  their  selection  is  further  discussed  in  Section  2.9.  in  relation 
to  the  overall  theory  of  penetration  failure,  and  the  merits  of  the  system  may  be  assessed  from 
the  analytical  results. 

The  fi  nal  solution  to  a  soil  structure  problem  is  most  valuable  if  the  results  can  be 
extended  over  the  full  range  of  soil  and  structural  parameters  and  presented  in  a  readily  usable 
manner.  Ibis  is  particularly  true  in  the  case  uf  soil-mai  bine  mechanics  where  a  general  picture 
of  performance  over  a  wide  range  of  soil  types  is  required.  There  are  three  groups  of  variables 
involved  in  a  soil  mechan.  s  problem  — the  soil  and  soil  to  structure  properties,  the  structural 
shape  described  by  angle,  and  other  ratios  of  its  dimensions  and,  finally,  the  size  of  the 
structure  and  the  associated  soil  masses.  The  complete  solution  has  to  be  such  that  an  actual 
problem  can  be  solved  by  means  of  a  set  of  graphs  or  tables  and  a  slide  rule. 

By  considering  the  cohesive,  weight  and  surcharge  terms  separately,  the  complete 
solution  may  be  presented  in  the  form  of  in  additive  equation.  Reece47  has  proposed  that  the 
most  complete  form  of  this  equation  is: 

F  (lb)  =  d2cNc  +  d  3->N-y  +  d2qNq  +  d2c,Na  . 2.1.3. 

where  d  represents  the  characteristic  dimension  of  the  structure.  This  equation  may  be  used 
to  describe  the  force  required  to  fail  soil  by  any  kind  of  structure  whatsoever.  The  fourth  teim 
takes  into  account  the  adhesion  between  the  soil  and  the  s:ructure.  The  N  factors  are  functions 


not  only  of  0  but  also  of  the  soil  structure  friction  and  the  shape  of  the  structure.  It  should 
also  be  n'  ted  that  the  surcharge  term  must  not  be  confused  with  the  weight  of  the  scil 
involved  in  tile  soil  flow.  It  will  generally  refer  to  some  loading  on  the  soil  which  does  not 
contribute  any  shearing,  for  example  the  loading  applied  by  a  smooth  surface  foundation  of 
some  adjacent  building. 

The  additive  nature  of  this  equation  follows  from  the  addition  of  the  various 
components  of  the  soil  force,  and  the  dimensionless  parts  are  necessary  in  order  to  provide 
dimensional  homogeneity.  The  weakness  of  this  equation  is  that  the  N  factors  are  not 
independent  of  the  values  of  the  dimensionless  terms  in  the  equation.  The  use  of  the 
equation  will,  therefore,  involve  an  error,  which  can  he  computed  and  must  be  presented  in 
a  complete  analysis  of  the  problem.  Another  possibility  is  to  compute  the  N  factors  for 
different  proportions  of  the  four  sources  of  force  on  the  structure  and  it  is  shown  later  that 
this  approach  is  applicable  to  the  case  of  the  penetration  resistance  for  a  deep  footing. 

In  a  recent  analysis  of  the  retaining  wall  problem48,  which  includes  surface  bearing 
capacity,  the  following  scheme  was  adopted: 

My  calculated  with  c/yz  =  q/yz  =  ca/yz  =  0 

Nc  calculated  with  c/yz  ~  1  ,  q/yz  =  ca/yz  =0  * 

Na  calculated  with  c/yz  =  ca/yz  =  1,  q/yz  -  0 

Nq  calculated  with  c/yz  =  1  ,  c/yz  =  ca/yz  =  0 

The  My  and  Nq  factors  are,  therefore,  correct  for  a  dry  sand,  whilst  the  Nc  factor  is 
computed  for  a  soil  with  weight  but  no  adhesion  or  surcharge.  The  Na  factor  is  calculated  for 
a  soil  with  weight  and  cohesion,  which  is  reasonable  since  it  is  not  possible  to  have 
adhesion  without  cohesion.  The  actual  value  of  c/yz,  ca/yz  and  q/yz  =  1  was  chosen 

for  computational  simplicity  and  it  was  later  shown  that  the  errors  would  not  exceed  5%  over 
the  complete  range  of  likely  combinations  of  these  properties, 

Recent  work  has  cast  doubt  on  the  necessity  for  the  minimum  value  postulate  in  soil 
mechanics  although,  initially,  tremendous  effort  was  expended  in  obtaining  it  for  the  computa¬ 
tion  of  Terzaghi's  and  Meyerhof’s  bearing  capacity  factors  in  this  programme  of  research  and 
in  the  analysis  of  the  retaining  wall  problem48.  These  minimising  methods  would  seem  better 
based  upon  the  following  simple  common  sense  proposition: 

If  a  family  of  failure  surfaces  can  be  found  which  everywhere  satisfy  the  conditions 
of  equilibrium,  the  boundary  conditions  and  the  idealised  soil  properties,  then  the 
soil  will  fail  on  that  surface  which  requires  the  minimum  force. 

In  fact,  such  a  family  of  failure  surfaces  do  not  exist  and,  by  establishing  the  only  possible 
failure  surface  which  satisfies  these  conditions,  the  minimising  procedures  are  not  applicable. 

For  example,  in  the  retaining  wall  problem,  the  correct  failure  surface  may  be 
selected  by  considering  the  stress  conditions  at  the  interface.  For  a  purely  frictional  soil 
with  known  interface  friction,  the  correct  stress  conditions  for  an  element  of  soil  adjacent 
the  structure  may  be  determined  from  Mohr's  circle  (Fig  2.1.1.).  The  orientation  of  the  slip 
planes  at  the  structure,  angles  &  and  b\  may  then  be  deduced.  By  this  means,  the  correct 
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position  of  the  spiral  pole  which  gives  the  correct  stress  conditions  on  the  element  at  the 
soil-structure  interface  and  at  the  passive  Rankine  zone  may  be  located.  It  is  recognised  that 
the  location  of  the  pole  by  this  technique  does  not  eliminate  stress  discontinuities  along  the 
full  length  of  the  interface,  but  by  thi  selection  of  an  element  at  some  point  other  than  at  the 
tip  of  the  wedge,  the  discontinuities  can  be  minimised.  Alternatively,  the  correct  stresses 
along  the  full  length  of  the  interface  may  be  obtained  by  interposing  a  zone  of  plane  shear 
between  the  interface  and  the  zone  of  radial  shear  (Fig  2.1.2.)-  This  is  more  applicable  to 
the  case  of  a  wedge  shaped  foundation.  In  Mintskovskv’s  analysis49,  for  example,  this 
modification  does  not  result  in  serious  deviations  from  the  original  criteria  assumed  for  a 
sliding  surface. 

A  further  problem  arises  with  the  addition  of  cohesion  and  adhesion  because  the 
stresses  become  indeterminate  unless  a  value  for  the  normal  stress  is  assumed  (Fig  2.1.3.). 
Variation  of  the  normal  stress  in  the  range,  0,5yz  to  3yz,  makes  little  difference  to  the  angle 
between  the  soil  structure  interface  and  the  slip  plane.  For  this  reason,  an  arbitrary  value  of 
yz  was  applied  to  Mohr’s  circle  in  the  investigation.  The  positions  of  the  spiral  centre 
derived  from  the  analysis  were  compared  with  those  obtained  for  the  minimising  procedure  in 
reference  48.  The  correlation  was  very  good  indicating  that  the  minimum  value  postulate 
degenerates  into  a  very  complicated  method  of  determining  the  correct  stress  conditions. 

Adopting  this  technique  in  the  determination  of  the  failure  geometry  for  the  surface 
bearing  capacity  of  a  continuous  flat  footing,  the  correct  stress  conditions  are  satisfied  only 
by  locating  the  pole  of  the  spiral  at  the  comer  of  the  footing,  Furthermore,  it  will  be  shown 
in  a  more  detailed  investigation  of  the  unrestrained  pole  (Section  5)  that  the  nature  of  the  bear¬ 
ing  capacity  calculation  is  such  that  it  cannot  tolerate  unrestraint  and  always  tends  to  adopt 
a  pole  close  to  the  comer  of  the  footing  regardless  of  the  complexity  of  the  minimising 
procedures. 

2.2  Failure  Geometry 

The  soil  displaced  by  a  moving  structure  will  always  endeavour  to  minimise  the  energy- 
lost  in  flow.  It  is  quite  prepared  to  modify  the  shape  of  the  structure  to  make  it  more  stream¬ 
lined,  by  arranging  a  fixed  soil  body  between  the  structure  and  the  flowing  soil.  A  fixed  soil 
body  will  form,  if  the  energy  losses  are  reduced  by  a  less  angular  flow  path,  by  a  larger  radius 
of  turn  or  where  shearing  flow  can  be  replaced  by  rigid  body  motion.  Its  actual  shape  and  size 
will  be  determined  by  the  fact  that  it  is  made  of  soil  and  must  satisfy  the  requirements  of  the 
ideal  soil  propeities,  the  soil  structure  properties  and  the  conditions  of  equilibrium.  The 
existence  of  such  a  body  covering  the  structure  surface  will  greatly  simplify  the  problem  if 
for  no  other  reason  that  it  eliminates  two  of  the  variable  (i.e.  i  -  0  and  ca  =  c).  The 
determination  of  the  existence  and  shape  of  such  a  body  is  therefore  the  starting  point  for  the 
solution  of  any  problem. 

When  a  continuous  footing  is  forced  vertically  into  a  uniform  soil  under  ideal  condi¬ 
tions,  the  flow  of  the  displaced  soil  must  be  symmetrical  about  the  vertical  axis  through  the 
centre  of  the  footing  in  order  to  minimise  the  lateral  soil  disturbance.  As  two  sharp  directional 
changes  in  flow  are  imparted  to  the  soil  beneath  the  structure  when  it  is  forced  to  penetrate,  a 
body  of  soil  attached  to  the  base  of  the  footing  could  modify  the  structure  shape  and  reduce 
the  enrgy  losses  by  providing  a  less  angular  flow  path  if  the  conditions  of  equilibrium  are 
fulfilled  (Fig  2.2.1.). 
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The  shear  stress  at  any  point  on  the  base  of  a  perfectly  smooth  footing  is  zero,  and 
the  maximum  stress  must  be  normal  to  the  base  for  the  footing  to  penetrate  into  the  soil. 

Thus,  the  base  of  the  footing  is  a  major  principal  plane  with  lines  of  incipient  failure 
inclined  at  angle  of  "/ /4  +  0/2  to.  the  horizontal.  As  there  are  no  shear  stresses  on  a  vertical 
plane  through  the  axis  of  the  footing  due  to  symmetry,  the  slip  lines  ar*  also  inclined  at  angles 
of  tt/4  +  0/2  at  this  point.  Assuming  that  the  body  of  trapped  soil  is  bounded  by  straight 
failure  planes,  a  wedge  with  sides  inclined  at  angles  of  77/4  +  0/2  to  the  horizontal  is 
formed  beneath  the  footing  d  ig  2.2,1b.).  Hue  to  lateral  constraint  of  the  soil  forces,  it 
remains  permanently  in  a  state  of  elastic  equilibrium  and  is  unaltered  by  sinkage.  There  is 
no  reason  to  suppose  that  curved  flow  lines,  requiring  at  least  one  point  of  inflexion  to  join 
the  points  of  known  stress,  could  replace  the  straight  lines.  This  would  not  only  yield  more 
abrupt  changes  in  flow  but  also  impart  a  rather  peculiar  slip  line  field  to  the  adjoining  soil 
mass  (big,  2.1.1c.).  On  the  same  grounds,  two  wedges,  moving  laterally,  are  also  infcr'or 
because  the  flow  path  must  then  revert  to  one  with  sharply  angular  changes  in  direction 
(big  2.2. Id.). 

If  the  base  of  the  footing  is  rough,  shear  stresses  of  unknown  magnitude  and  distribu¬ 
tion  could  be  mobilised.  However,  if  the  wedge  is  generated  and  is  in  equilibrium  beneath 
a  perfectly  smooth  footing  without  the  development  of  any  shear  stresses,  it  is  not  conceiv¬ 
able  that  base  roughness  would  in  any  way  influence  the  shape  of  the  fixed  soil  body.  There 
is  no  need  for  shear  stresses  to  be  developed  along  the  base,  and  there  is  no  way  of  mobilis¬ 
ing  them  with  symmetrical  loading  under  ideal  conditions.  The  base  of  the  structure  is  no 
longer  important  and  is  replaced  by  the  two  soil  boundaries  of  the  wedge  both  of  which  are 
failure  planes.  The  directions  of  the  other  families  of  slip  lines  are  known  from  the  conditions 
of  stress  at  the  wedge  tip  and  must  emanate  at  angles  of  77/2  +  0  to  the  wedge  boundaries. 

As  the  geometry  is  symmetrical  about  the  central  axis  of  the  footing,  the  failure  surface  on 
only  one  side  is  subsequently  considered. 

The  particular  case  when  sinkage  in  zero  and  the  footing  smooth  is  capable  of 
rigorous  solution  by  the  method  of  Static  Analysis,  and  this  can  be  used  as  the  starting  point 
in  the  development  of  the  failure  pattern,  big  1.2.4a.  shows  Sokolovski’s1  solution  for  a 
foundation  obliged  to  fail  on  one  side  only.  The  zone  BDI£  is  always  straight  sided  (Hankine) 
as  long  as  the  surcharge  is  vertical.  The  zone  ABC  is  doubly  curved  but  as  long  as  All  is 
smooth,  angles  CAB  and  ABC  must  always  be  77/4  +  0/2.  Angle  ACB  is  77/2  -  0  under  any 
circumstances.  Although  the  case  of  symmetrical  failure  was  not  considered,  it  is  not 
difficult  to  see  that  symmetry  requires  that  C  shall  lie  below  the  mid  point  of  AB,  and  the 
angular  requirements,  mentioned  above,  force  ABC  to  become  straight  sided.  The  curved  fan 
BCD  then  becomes  a  radial  shear  zone  bounded  by  a  logarithmic  spiral  and  the  total  failure 
pattern  is  shown  in  Fig  2.2.2.  This  conclusion  is  confirmed  by  Josselin  de  Jong4.  The 
wedge  is  the  same  shape  as  had  been  deduced  already. 

For  the  more  general  case  of  shallow  and  deep  foundations,  there  must  be  a  logical 
transition  from  the  failure  geometry  of  the  limiting  surface  case.  There  are  two  possible 
modifications  which  can  be  made  to  the  failure  geometry  to  account  for  sinkage.  First  the 
shape  and  extent  of  the  radial  shear  zone  may  remain  unaltered,  and  the  boundary  DF  of  the 
plane  shear  zone  may  be  produced  to  meet  the  surface  at  the  correct  angle  of  tt/4  -  0/2 
(Fig  2.2.3.).  This  configuration  is  obviously  possible  at  very  small  sinkages  and  has  much 
to  commend  it.  For  instance,  the  correct  stresses  are  attained  along  the  free  surface  and 
along  the  smooth  foundation  shaft,  and  the  boundary  conditions  are  nowhere  violated. 
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It  does  not  apply  to  a  footing  because  there  is  no  shaft  support  to  cause  the  passive 
pressure.  It  does  not  even  apply  to  a  foundation  because  it  has  been  shown  by  experiment 
that  pressing  the  footing  down  does  not  cause  the  soil  to  exert  full  passive  pressure  against 
the  walls.  Also,  the  rupture  distance,  GE,  does  not  increase  linearly  with  depth,  and  conse¬ 
quently  this  type  of  failure  may  be  rejected.  Alternatively,  the  zone  of  shearing  flow  may  be 
extended  in  conjunction  with  a  change  in  the  shape  of  the  plane  shear  zone  (Fig  2.2.4a.).  The 
rupture  plane,  BD,  is  the  last  radial  slip  line  of  the  radial  shear  zone  BCD.  The  plane  shear 
zone  BDEG  is  delineated  on  two  sides  by  a  pair  of  failure  surfaces,  BD  snd  DE,  making  an 
angle  of  -m/2  +  0  with  each  other.  The  boundary.  DE,  no  longer  satisfies  the  stress  conditions 
at  the  surface  but  this  requirement  can  be  fulfilled  without  adding  complexity  to  the  failure 
geometry,  by  a  minor  re-orientation  of  the  slip  plane,  DE  '. 

The  failure  surface,  CDE,  in  Fig  2.2.4a.  is  bounded  by  a  logarithmic  spiral  and  a 
straight  line.  In  a  dilating  material,  it  would  probably  be  bounded  by  curves  such  as  are  shown 
by  dotted  line  CDE.  The  choice  of  the  dimensions  of  BCDEG  by  means  of  the  Coulomb 
criterion  and  static  equilibrium  ensures  that  the  surface  is  little  different  from  the  real  one; 
the  approximation  being  used  to  simplify  the  mathematics  involved. 

The  failure  geometry  can  be  extended  logically  to  deep  foundations  by  further  extension 
of  the  radial  shear  zone.  The  limiting  case  depends  on  the  nature  of  the  support  and  the  inter¬ 
face  properties  along  the  foundation  shaft  and  on  the  soil  properties.  In  very  strong  cohesive 
soils,  where  the  influence  of  soil  weight  may  be  neglected  and  the  foundation  shaft  is 
unsupported,  it  is  assumed  that  the  failure  pattern  proposed  by  Jaky8  occurs  (Fig  1.3.2.).  For 
any  other  cohesive-frictional  soil  where  weight  is  important,  this  type  of  failure  is  largely- 
discounted.  It  is  more  likely  that  the  mode  of  failure  reaches  the  limit  shown  in  F'ig  2.2.4b., 
then  changes  to  that  described  by  VcsicM  at  deep  sinkage  when  the  mass  of  soil  affected  by 
the  penetration  of  the  footing  becomes  sufficiently  great  to  ‘absorb’  the  displaced  soil  by  soil 
compaction  U  ig  2.2.5.).  Yesic  has  suggested  the  possible  upper  limit  for  general  shear 
failure  in  sand  occurs  when  the  skin  resistance  distribution  changes  from  hydrostatic  to 
uniform.  The  critical  penetration  ratio  is  dependent  on  density  and  0.  The  importance  of  this 
change  in  the  mode  of  failure  only  becomes  significant  in  the  present  study  of  incompressible 
soils  when  0  is  small.  It  is  important  to  note,  however,  that  the  effect  of  soil  compressibility- 
on  the  mode  of  failure  even  in  the  highly  frictional  materials  tested  is  to  reduce  the  penetration 
ratio  at  which  this  changeover  occurs. 

As  the  pressure  exerted  on  the  foundation  shaft  constitutes  one  of  the  unknown  para¬ 
meters,  the  problem  is  simplified  by  assuming  that  the  vertical  wall  BG  is  self-supporting  and 
represents  a  free  surface.  From  practical  observations,  this  is  quite  justifiable  over  most  of 
the  soil  range  to  depths  in  excess  of  those  investigated.  It  is  also  particularly  relevant  to 
vehicle  running  gear.  One  major  exception,  however,  it  a  dry  frictional  material.  It  is  proposed 
therefore,  to  develop  a  theory  of  penetration  failure  on  the  following  basis: 

a.  Dry  frictional  material— to  include  the  effect  of  the  shaft  support  or  of  material 

falling  into  the  shaft,  forming  a  soil  to  soil  boundary  at  the  central  axis  of  the  footing. 

b.  Cohesi ve-frictio/al  material — to  assume  the  shaft  is  unsupported,  but  to  make 

provision  for  further  extension  to  a  supported  shaft. 
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The  maximum  height  of  the  unsupported  shaft  for  a  cohesive-frictional  soil  may  be 
evaluated,  approximately,  by  considering  the  stability  of  a  vertical  bank.  It  can  be  shown  that 
the  critical  height,  iic,  is  given  by  the  equation: 

Hc  =  4(c/>)  tan  (45  +  0/2)  . 2.2.1. 

and  by  substituting  values  for  the  soil  density  and  0  of  0.06  lb  in-3  and  40°  respectively,  a 
42  in  high  bank  is  stabilised  when  the  cohesion  is  0.03  lb  in  2.  This  height  of  bank  is 
equivalent  to  a  penetration  ratio  of  3  for  a  14  in  wide  vehicle  track  arid  represents  the 
upper  limit  of  sinkage;  for  different  values  of  the  parameters  the  value  of  Hc  is  given  in 

Fig  2.2.6. 

The  zone  BDF.G  will  tend  to  behave  as  a  solid  body  (Fig  2.2.4a.).  Of  the  two  failure 
boundaries,  one  must  form  the  last  radial  slip  line  of  the  radial  shear  zone  and  therefore  must 
emanate  from  the  comer  of  the  footing;  the  other  must  make  an  angle  of  n/2  +  0  to  it  and  be 
a  tangent  to  the  curved  slip  line,  thus  forming  the  continuous  outer  boundary  of  the  failure 
surface.  These  two  plane  failure  surfaces  need  not  be  of  the  same  length.  From  the  conditions 
of  equilibrium  there  is  an  infinite  number  of  possible  configurations  for  the  zone  which  are  a 
function  of  the  sinkage  and  the  soil  parameters.  Only  one  geometrical  shape  fulfills  the 
condition  that  the  junction  of  the  two  failure  surfaces  must  lie  on  the  spiral  surface  which 
then  determines  the  failure  geometry  for  the  particular  parameters  which  have  been  selected. 

As  the  determination  of  the  failure  geometry  depends  on  the  equations  of  equilibrium,  it  is 
fully  covered  in  the  next  sub-section. 

Summarising  the  important  aspects  of  the  failure  geometry,  it  is  truly  symmetrical  about 
the  central  axis  of  the  footing  and  consists  of  only  one  figure  of  rupture  for  a  soil  with  both 
weight  and  cohesion.  It  comprises  three  different  zone  types: 

a.  The  fixed  soil  body  is  a  single  wedge  with  base  angle  of  77/4  +  0/2  regardless 
of  base  roughness. 

b.  The  zone  of  shearing  flow  is  formed  of  a  family  of  radii  and  logarithmic  spirals 
whose  pole  lies  at  the  corner  of  the  footing. 

c.  The  plane  shear  zone  is  bounded  by  one  free  surface,  the  vertical  wall  of  the 
shaft,  and  two  failure  surfaces.  These  surfaces  make  an  angle  of  rr/2  +  0  with  each 
other  and  one  determines  the  extent  of  the  zone  of  shearing  flow.  The  shape  and 
extent  of  the  plane  shear  zone  is  obtained  from  the  conditions  of  equilibrium. 

The  outer  boundary  is  a  smooth  curve  to  satisfy  the  flow  criterion.  The  size  of  the  failure 
geometry  is  determined  by  0,  the  penetration  ratio,  the  cohesion  to  weight  ratio  and  the  shaft 
parameters. 

2.3.  The  Calculation  of  Penetration  Resistance 

The  calculation  of  the  penetration  resistance  of  a  continuous  footing  may  be  considered 
in  three  parts.  First,  the  figure  of  rupture  is  obtained  by  considering  the  equilibrium  of  the 
straight  sided  zone.  The  magnitude  and  direction  of  the  soil  forces  due  to  the  weight  of  the 
soil  mass  within  the  failure  zone,  and  those  acting  on  the  boundaries  are  determined.  Finally, 
knowing  the  directions  of  the  resultant  forces,  the  penetration  resistance  is  evaluated.  The 
general  procedure,  assumptions  and  approximations  are  discussed  in  the  following  sections  and 
the  complete  theoretical  analysis  is  given  in  the  appendix. 
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It  is  assumed  that,  for  static  equilibrium  in  a  cohesionless  soil,  the  weight,  \\1 ,  of  the 
plane  shear  zone  BDEG  is  balanced  by  two  weight-frictional  boundary  forces  Fyj  and  Fy2 
on  any  pair  of  failure  planes,  ED  and  DB  respectively  and  by  a  shaft  force  Fys,  on  BG 
(Fig  2.3.1.).  The  weight  acts  vertically  downwards  the  boundary  forces  act  at  0  to  the  normal 
of  the  plane  and  the  shaft  force  acts  at  an  angle  8  to  the  normal.  The  two  boundary  forces  are 
governed  by  the  same  boundary  coefficient  of  earth  pressure,  k^y  modified  in  accordance  with 
equation  2.1.1.  and  the  normal  stresses  are  distributed  in  a  way  which  diverges  slightly  .Tom 
the  hydrdstatic  manner.  The  shaft  force  is  governed  by  the  coefficient  of  earth  pressure  kys 
which  is  assumed  to  have  a  known  value  and  by  a  hydrostatic  stress  distribution.  Considering 
any  line,  BD,  sloping  at  an  angle,  £  ,  to  the  horizontal,  the  three  equations  for  horizontal, 
vertical  and  rotational  equilibrium  are  respectively: 


F yj 

cos  £ 

-  Fy  2 

sin  (0  -  0  -  Fy  s  cos  ? 

.  2.3.1. 

Fyj 

cos  (0 

-  £)  - 

Fi  sin  £  =  Wj  +  FyS  sin  8 

.  2.3.2. 

Widj 

-  Fy, 

d3  - 

l‘y2d4  +  Fy  sdj  -  0 

.  2.3.3. 

In  simple  terms,  there  are  three  unknowns,  the  boundary  coefficient,  the  boundary- 
exponent  and  the  slope  angle.  These  are  functions  of  the  angles  of  internal  and  interface 
friction,  the  shaft  coefficient  of  earth  pressure  and  the  penetration  ratio.  However,  the 
relationships  are  extremely  complex  because  of  the  exponential  terms,  and  the  equations  can 
be  solved  only  be  substitution.  The  equations  are  written  in  general  terms  so  that  the  slope 
angle  can  take  any  value  in  the  range  -  (45  -  0/2)  <  £  90  when  the  sinkage  is  greater  th 

than  zero.  In  other  words,  the  extent  of  the  logarithmic  spiral  zone  may  vary  from  the  surface 
condition  up  to  the  full  extent  where  it  breaks  into  the  foundation  shaft.  Normally,  however, 
the  calculations  arc  terminated  when  the  value  of  the  slope  angle  reaches  0,  so  that  the 
failure  geometry  for  the  limiting  case  is  that  shown  in  l-  ig  2.2.4b. 

When  cohesion  is  added,  six  further  forces  must  be  considered  in  the  equations  of 
equilibrium.  First,  there  are  the  cohesive  forces,  Ci  and  C2,  acting  along  the  boundaries 
DE  and  BD,  respectively,  and  the  adhesive  force  along  the  shaft,  Cs-  The  magnitude  of 
these  forces  depend  on  the  length  of  the  relative  surfaces  and  the  soil  or  interface. properties. 
The  cohesive-frictional  boundary  forces,  Fcl  and  Fc2 ,  act  at  0  to  the  normal  of  the  planes 
DE  and  BD,  and  the  cohesive-frictional  shaft  fcrce,  Fcs,  acts  at  angle  8  to  the  normal. 

The  two  cohesive-frictional  boundary  forces  are  governed  by  the  same  boundary  co¬ 
efficient  of  earth  pressure  k/,  (modified  in  accordance  with  equation  2.1.2.)  and  the  forces  are 
distributed  in  a  way  which  diverges  slightly  from  a  uniform  manner.  The  cohesive-frictional 
shaft  force  is  governed  by  the  shaft  coefficient  of  earth  pressure,  kcsof  known  value  and  the 
force  is  uniformly  distributed.  Including  these  forces  in  the  equations  of  horizontal,  vertical 
and  rotational  equilibrium,  the  complete  relationships  are: 

Fy1  cos£  -  Fyj  sin  (0  -  ^)  r  Fcl  sin  (0  -  £>  = 

Fys  cos  8  +  pcs  cos  ?  +  C2  cos  £  -  Ci  sin  (0  -  £)  . 2.3.4. 

Fy2  cos  (0  -  £>  -  Fyx  sin  £  +  Fc2  cos  (0  -  ^)  -  Tcl  sin  f  = 

Wj  +  Fys  sin?  +  Fcs  sin?  +  Cs  +  Cj  cos  (0  -  £)  -  C2  sin  £ 

Ml  •  Pyid,  -  Fy2d4  +  F^s  ♦  Cj d7  -  Fcldg  -  Ec2  d9  +  Fcsd10  =  0 
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.  2.3.5. 

.  2.3.6. 


For  this  general  case,  there  are  four  unknowns — the  weight-frictional  and  cohesive- 
frictional  boundary  coefficients  of  earth  pressure,  the  boundary  exponent  and  the  slope  angle. 
These  depend  on  three  further  input  parameters,  the  cohesive-frictional  shaft  coefficient  of 
earth  pressure,  the  cohesion  to  weight  ratio  and  the  adhesion  to  weight  ratio,  in  addition  to 
those  already  listed  for  a  cohesionless  material.  The  parametric  complexity  may  be  greatly 
simplified  by  considering  an  unsupported  foundation.  It  avoids  the  problem  of  assessing 
independent  coefficients  of  earth  pressure  for  the  weight-frictional  and  cohesive-frictional 
shaft  forces  and  it  allows  the  interaction  between  the  remaining  variable  to  be  presented  more 
simply  in  graphical  fotm. 

v 

One  difficulty  remains.  In  order  to  determine  the  four  unknowns,  a  fourth  equation  must 
be  developed  to  relate  the  two  boundary  coefficients.  Alternatively,  three  of  the  unknowns  may 
be  obtained  in  terms  of  the  slope  angle  and  the  final  solution  minimised  with  respect  to  £• 

The  latter  course  was  rejected  because  of  the  excessive  amount  of  computation  which  would 
be  required.  The  following  relationship  is  adopted  to  relate  the  two  boundary  coefficients: 

k0c  =  (k0>)  M  .....2.3.7, 

This  relationship  is  correct  for  the  surface  case  and  has,  therefore,  some  theoretical  basis. 
When  the  penetration  ratio  is  greater  than  zero,  the  relationship  becomes  more  of  an  empirical 
weighting  factor  because  the  two  coefficients  are  dependent  on  the  cohesion  to  weight  ratio. 

It  is  shown  in  the  appendix  that  k^c  =  1  for  a  cohesive  soil  without  weight.  Hence,  the 
greatest  source  of  error  in  the  value  of  weight-frictional  boundary  coefficient  occurs  in  strong 
cohesive  soils  where  its  influence  is  minimal.  Similarly,  the  greatest  source  of  error  in  the 
value  of  the  cohesive-frictional  boundary  coefficient  occurs  when  the  cohesion  is  small,  so 
that  this  basis  of  the  analysis  does  not  represent  an  excessive  approximation. 

The  value  of  the  slope  angle,  determined  from  the  solution  of  equation  2.3.4.  to 
2.3.7.  inclusive,  fixes  the  shape  of  the  failure  geometry.  The  soil  weight  and  the  soil 
boundary  forces  may  be  determined  as  shown  in  the  appendix.  The  weight-frictional, 
cohesive-frictional  and  cohesive  forces  on  the  plane  1313  are  ignored  in  the  calculation  of 
the  resultant  forces  because  they  are  internal  soil  forces.  Thus  the  resultant  w  eight- 
frictional  force  acting  at  0  to  the  normal  of  the  wedge  boundary,  BC,  is  dependent  on  the 
weight  of  the  soil  within  the  failure  zone,  BCDEG,  and  on  the  weight-frictional  forces  act¬ 
ing  on  the  boundary,  CDE.  Similarly  the  resultant  cohesive-frictional  force  acting  at  0  to  the 
normal  cf  the  wedge  boundary,  BC,  is  dependent  on  the  cohesive-frictional  forces  and  the 
cohesive  forces  acting  along  the  side  of  the  wedge,  BC.  The  distribution  of  the  resultant 
forces  is  in  accordance  with  equation  2.1.1. 

The  weight-frictional  and  cohesive-frictional  forces  along  the  spiral  boundary  CD,  are 
difficult  to  determine.  As  these  forces  pass  through  the  pole  of  the  spiral  they  may  be 
eliminated  from  the  calculation  by  taking  moments  about  that  point.  The  resultant  forces  are 
given  by  the  expressions: 

Fy  =  [Wj.dj  +  W.d2  -  F7l.d3)/d6  . 2.3.8. 


Fc  =  (q.d7  -  Fcl.dg  +  Mc)/dn  . 2.3.9. 
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Where  the  foundation  shaft  is  supported  additional  terms  must  be  added  to  the 
numerators  of  equations  2.3.8.  and  2.3.9.  (c.f.  Appendix  A). 

Resolving,  vertically,  forces  Fy ,  Fc  and  the  'adhesion'  force,  C3  (all  of  which  act 
on  the  two  boundaries  of  the  wedge,  An  and  BC)  and  considering  the  weight  of  the  wedge, 
ABC,  the  resultant  penetration  resistance  of  the  foundation  is  obtained: 

F  lb  in*1  =  2Fy  cos  (</■■  -  0)  -  /j^B*  tan  41  +  2FC  cos  (\p  -  0)  +  2C3  sin  \p 

.....  2.3.10. 

This  expression  can  be  written  in  the  form  of  a  two  part  additive  equation: 

p  lb  in  -  'Ay  BNy  +  cNc  . 2.3.11. 

where: 

Ny  =  2Fy  cos  (y  -  0)/Cl'H32)  -  tan  .  2.3.12. 

Nc  -  2[Fc  cos  (v>  -  0)  +  C3  sin  y  /(cB)  . 2.3.13. 

The  dimensionless  factors,  Ny  and  Nc,  are  dependent  on  0,  the  penetration  and 
cohesion  to  weight  ratios.  The  implication  of  the  two  rat.os  is  derived  from  the  equilibrium 
of  the  plane  shear  zone,  and  for  a  buried  footing  there  is  no  surcharge  term.  Any  real 
surcharge  on  the  soil  surface  can  be  considered  separately. 

When  the  foundation  shaft  is  supported,  equation  2.3.11.  represents  the  base 
penetration  resistance  and  the  dimensionless  factors  are  further  dependent  on  the  shaft 
parameters.  The  total  penetration  resistance  is  the  summation  of  the  base  penetration 
resistance  and  the  surface  traction  on  the  foundation  shaft. 

2.4.  Computer  Programme 

The  computer  programme  is  divided  into  three  sections,  namely  for  a  cohesionless 
soil,  a  cohesive  soil  without  weight  and  a  cohesive  soil  with  weight.  Additional  procedures 
are  required  for  certain  individual  cases  such  as  when  either  0,  or  the  sinkage  or  the  slope 
angle  is  zero.  The  complete  block  diagram  is  shown  in  Fig  2.4.1. 

The  soil  parameters  (0,T,  c),  the  soil  structure  parameters  (S,  ca,  kys,  kcs)  and 
the  structure  parameters  (B,  z)  are  stated  on  the  input  data  and  are  thereafter  available 
throughout  the  programme. 

The  failure  geometry  is  determined  first.  The  three  transcendental  equations  of 
equilibrium,  which  are  functions  of  the  boundary  coefficient,  the  boundary  exponent  and  the 
slope  angle  ate  solved  by  substitution.  The  range  of  the  slope  angle  -(45-0/2)  <  £  <  0  is 
divided  into  an  arbitrary  number  of  20  steps.  For  each  value  of  the  slope  angle  the  value  of 
the  boundary  coefficient,  the  boundary  exponent  and  the  difference  or  error  between  the 
L.ll.S.  and  R.ll.S.  of  the  equation  of  rotational  equilibrium  are  determined.  When  two  values 
of  the  slope  angle  have  beer,  selected  such  that  the  error  is  zero  within  the  interval,  the 
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correct  solution  for  the  slope  angle  is  obtained  by  an  iterative  process.  The  range  is 
bisected  and,  after  each  cycle,  the  section  not  containing  the  solution  is  discarded  until  the 
final  answer  is  achieved  with  a  prescribed  accuracy  of  0.5  x  10  ‘6  after  about  twenty  iterations 

In  the  case  of  a  cohesive  soil  with  weight  the  same  procedure  is  used  to  determine  the 
value  of  the  boundary  exponent  by  substitution  for  each  value  of  the  slope  angle.  This 
additional  complexity  is  not  required  for  the  limiting  values  of  the  cohesion  to  weight  ratio. 

Once  the  failure  geometry  has  been  determined,  the  final  calculation  of  the  soil  forces 
and  the  penetration  resistance  is  straightforward.  No  minimising  techniques  are  included  in 
the  final  form  of  the  programme.  However,  provision  was  made  in  earlier  work  for  investigating 
the  effect  of  unrestrained  spiral  pole  and  of  a  variable  wedge  angle  on  the  penetration 
resistance  using  simplex52  (c.f.  Section  5.). 

2.5.  Computed  Results  for  a  Dry  Friction  Material  with  Weight 

The  dimensionless  Ny  factors  for  three  values  of  the  shaft  coefficient  of  earth 
pressure  are  shown  in  Figs  2.5.1.  -  2.5.3.  The  Ny  factor  for  each  value  of  the  shaft  co¬ 
efficient  is  completely  represented  by  two  families  of  curves,  one  for  the  perfectly  smooth 
shaft  interface  ( 8  -  0)  and  the  other  for  a  perfectly  rough  shaft  interface  (8  =  0).  The  data  in 
Fig  2.5.1.  is  included  because  the  value  of  unity  is  commonly  selected  for  the  shaft  co¬ 
efficient  to  provide  a  rapid  solution  to  practical  problems.  The  additional  data  in  Figs  2.5.2. 
and  2.5.3.  provide  the  necessary  interpolation  material  for  intermediate  values  of  the-  shaft 
coefficient. 

The  presence  of  interface  friction  is  shown  to  increase  the  base  penetration  resistance 
of  the  foundation  in  proportion  to  the  downward  thrust  acting  on  the  radial  boundary  of  the  zone 
of  shearing  flow.  1  lie  magnitude  of  this  effect  is  also  dependent  on  the  shaft  coefficient  which 
may  take  any  value  between  the  coefficient  of  active  earth  pressure,  ka,  and  the  coefficient  of 
passive  earth  pressure  kp.  In  practice,  however,  it  is  statically  impossible  to  obtain  a  solu¬ 
tion  for  a  perfectly  rough  interface  using  the  maximum  value  of  the  shaft  coefficient  because 
soil  failure  would  occur  by  lateral  compression  instead  of  lv  vertical  compression.  Hence,  the 
shaft  coefficient  was  'cstricted  in  value  to  two-thirds  of  that  for  the  coefficient  of  passive 
earth  pressure.  Increasing  the  value  of  the  shaft  coefficient  up  to  this  limit  increases  the  base 
peneiration  resistance.  For  the  smooth  shaft  interface,  however,  increasing  the  value  of  the 
shaft  coefficient  results  in  a  decrease  in  the  base  penetration  resistance.  This  occurs  because 
the  hor  zontal  shaft  force  generates  a  less  than  hydrostatic  stress  distribution  on  the 
boundaries  of  the  plane  shear  zone  and  effectively  increases  the  anti-clockwise  moment  of  the 
boundary  force,  Fyi  .  This  condition  would  not  normally  be  achieved  because  it  is  unlikely 
that  the  interface  friction  would  be  zero  except  when  the  shaft  coefficient  was  at  its  minimum 
value. 


The  rupture  distance  for  the  three  values  of  the  shaft  coefficient  are  shown  un 
Figs  2.5.4.  to  2.5.6.  As  the  rupture  distance  is  proportional  to  the  breadth  of  the  footing, 
the  dimensionless  ration  is  f/D  where  f  is  the  horizontal  distance  between  the  edge  of  tiie 
foundation  shaft  and  the  point  at  which  the  failure  plane  breaks  out  at  the  soil  surface.  1  he 
effect  of  increasing  either  the  angle  of  interface  friction  or  the  shaft  coefficient  is  to  increase 
the  rapture  distance  because  it  depresses  the  slope  angle  of  the  radial  boundary.  U hen  the 
shaft  is  smooth,  the  failure  geometry  tends  towards  the  statically  correct  stress  field  for 
the  li  miting  condition  in  which  full  passive  earth  pressure  is  developed  on  the  shaft. 
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2.6. 


Computed  Results  for  a  Frictional-Cohesive  Material  without  Weight 


The  values  of  Nc  for  a  weightless  soil  are  shown  in  Fig  2.5.7.  In  the  absence  of 
weight,  the  shear  bulb  can  be  formed  around  the  base  of  the  foundation.  Thus,  the  bearing 
capacity  reaches  a  constant  value  which  is  independent  of  any  further  sinkage.  This  explains 
why  the  curves  for  Nc  have  a  rather  unusual  form. 

The  rupture  distances  fall  to  zero  wl^p  tiie  shear  bulb  is  fullv  formed  (Fig  2.5.8.). 
Thereafter,  the  rupture  distance  f  is  taken  as  'he  vertical  distance  between  the  base  of  the 
footing  and  the  point  at  which  the  boundary  of  the  shear  bulb  breaks  into  the  vertical  shaft  of 
the  foundation. 

2.7.  Computed  Results  for  a  Frictional-Cohesive  Material  with  Weight 

The  Ny  factors,  Nc  factors  and  rupture  distances  for  the  intermediate  c/7 B  ratios 
are  presented  for  three  different  shaft  conditions.  The  relevant  data  for  the  unsupported  shaft 
are  shown  in  Figs  2.7.1.  to  2.7.3.;  those  fora  perfectly  rough  shaft  (i  =  ca  =  0)  are 
presented  in  Figs  2.7.4.  to  2.7.6.  and  those  for  a  perfectly  rougii  shaft  (S  -  0,  c  a  -  c)  are 
given  in  Figs  2.7.7.  to  2.7.9.  The  calculation  of  the  dimensionless  groups  for  the  foundation 
with  smooth  and  rougii  shafts  were  based  on  an  arbitrarily  chosen  value  of  the  shaft  coefficients 
kys  and  k^.^  of  unitv. 

The  Ny  factors  follow  the  expected  pattern.  For  the  smooth  shaft,  the  curves  for 
c/TB  =  0.1  are  very  simiiar  to  those  for  cAlJ  =  0,  with  the  curves  for  the  unsupported 
shaft  and  the  rough  shaft  lower  and  higher,  respectively.  The  genera!  effect  of  increasing  the 
c /y B  ratio  is  to  increase  the  value  of  the  Nv  factor  until  suc  h  time  as  equilibrium  can  no 
longer  be  achieved,  and  the  curves  tend  very  rapidly  to  zero. 

The  Nc  factors  show  comparable  trends.  For  the  unsupported  shaft  the  curves  for 
c/7B  =  100  are  very  similar  to  those  for  c,  /B  =  «.  As  the  c/>B  ratio  decreases  for  all 
shaft  conditions,  the  curves  for  the  Nc  factors  exhibit  either  an  initial  reduction  in  gradient 
or  a  negative  gradient  with  increasing  penetration  ratio.  This  becomes  more  pronounced  at 
the  higher  values  of  0.  Although  this  effect  may  be  explained  in  tire  case  of  the  unsupported 
foundation  shaft  as  the  negative  cohesive  force  required  to  replace  the  lateral  support  of  the 
shaft,  the  continued  presence  of  the  phenomenon  when  the  shaft  is  supported  would  appear  to 
be  due  to  the  inflexible  coupling  between  the  weight-frictional  and  cohesive-frictional 
boundary  coefficients. 

The  variation  of  the  rupture  distances  for  the  intermediate  c/ '13  ratios  is  a  reflection 
of  the  changes  in  N  factors  with  a  gradual  reduction  between  the  limits  of  those  for  a  cohesion¬ 
less  soil  and  those  for  a  weightless  soil  when  the  shaft  is  unsupported.  More  significant 
changes  occur  for  the  different  shaft  conditions.  The  influence  of  shaft  support,  which  is 
augmented  by  the  presence  of  shaft  traction  is  to  increase  tiie  rupture  distance  by  depressing 
the  slope  angle  of  the  radial  boundary  of  the  zone  of  shearing  flow. 
f 

2.8.  Interpolation 

Since  0  is  represented  by  a  family  ot  curves  increasing  in  steps  of  5°,  the  main 
interpolation  is  required  to  cover  the  N  factors  at  intermediate  values  of  0.  Interpolation  is 
also  required  for  the  influence  of  the  shaft  parameters  each  of  which  is  accounted  for  at  only 
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a  few  specified  values.  Interpolation  lias  not  been  attempted  for  the  data  presented  in 
section  2.7.  because  of  the  larger  number  of  variables  and  because  of  the  complexity  of  the 
interaction  between  the  various  parameters. 

For  any  given  penetration  ratio,  the  variation  of  N(Ny  or  Nc)  with  0  can  be 
expressed  with  considerable  accuracy  by: 

N  =  u  e  v  tan  0  .  2.8.1. 

where  u  and  v  are  constants  dependent  on  the  penetration  ratio,  f  or  small  values  of  0  radians, 
tan  0-0  (the  error  of  such  an  approximation  being  only  0.2%  for  0  =  +  5°).  lienee,  in  the 
5°  intervals  it  follows  that: 

N  =  u  e  v  ^  ....  2.0.2. 

Since  the  ordinates  in  the  N  factor  graphs  are  logarithmic  scales,  a  linear  interpola¬ 
tions  between  the  5°  curves  yields  intermediate  values  of  N  to  a  high  degree  of  accuracy.  The 
intermediate  values  of  the  rupture  distance  may  be  obtained  in  a  similar  manner. 

For  the  calculation  of  the  Ny  factors  at  intermediate  values  of  5,  the  angle  of  inter¬ 
face  friction  is  considered  as  a  fraction  of  the  angle  of  internal  friction  and  this  ratio  varies 
from  y.em  for  a  smooTi  shaft  to  unity  for  a  perfectly  rough  one.  The  variation  of  the  Ny 
factor  ith  the  S/0  ratio  is  substantially  linear  (Fig  2.8.1.).  The  following  relationships 
can  tL..  be  established: 
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wlu.ro  Nyj,  is  the  requited  value  of  Ny  and  Ny0  and  Ny^  are  the  corresponding  values  of 
Ny  at  S  0,  respectively.  A  similar  form  of  equation  2.9.3.  may  be  derived  for  the  inter¬ 
mediate  values  of  the  f/B  ratio  which  exhibit  the  same  trend  with  respect  to  the  S/0  ratio 

(Fig  2.0.2. ). 


The  Ny  factor  is  also  shown  to  exhibit  a  linear  relationship  with  variation  of  the 
shaft  coefficient  except  at  the  lower  values  of  0  which  are  not  di recti y  applicable  to  a 
cohes ionlcss  soil  (Fig.  2.8.3.).  The  form  of  the  interpolation  equation  is  identical  to  that  for 
the  5>/0  ratio. 
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wheic  Ny*  is  required  value  of  Ny  uml  Nya  ami  Nyp  are  th  .icspondng  values  of  N  y 
for  the  iiiiiiiiiri.uin  and  maximum  values  of  the  shaft  coefficient  ol  earth  pressure,  'i  hi;,  is  also 
applicable  to  the  f/B  ratio  (Fig  2.8.4.). 
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The  application  of  the  interpolation  procedures  is  presented  in  some  detail  for  the 
following  soil  and  structural  properties. 

Soil  Properties:  0  -  40.5  ,  7  -  0,06  1b  in'3 

Soil-structure  interface  variables:  8  =  32°,  kys  -  1.23 

Structural  parameters:  13  -  2  in,  z  =  6  in 

The  relevant  values  of  \y  factor  may  be  obtained  by  three  successive  interpolation 
for  kys,  8  and  lastly  for  0.  Interpolation  for  the  shau  coefficient  must  be  considered  first 
because  the  limiting  values  of  the  shaft  coeflicient  are  not  the  same  for  the  rough  and  smooth 
shaft  conditions.  Improved  accuracy  can  also  be  achieved  by  interpolating  for  a  shaft  co¬ 
efficient  of  1.23  between  the  value  of  unity  and  the  upper  limiting  value  rather  than  between 
the  value  of  the  coefficient  of  active  earth  pressure  and  the  upper  limiting  value.  The  eight 
initial  values  of  Ny  are  obtained  from  higs  2.5.1.  and  2.5.3.  and  listed  in  Table  1,  together 
with  the  interpolated  Ny  values  at  each  successive  step. 


TABLE  1 


The  di  Iterance  between  the  interpolated  N-.  factor  of  745  and  the  computed  Ny  factor 
of  736  represents  uri  error  of  approximately  lVc.  7 


One  of  the  problems  encountered  in  predicting  the  penetration  resistance  for  a  footing 
is  that  the  shaft  parameters  are  indeterminate.  In  the  absence  of  any  vertical  restraint,  the 
sand  falls  into  the  shaft  and  forms  a  vertical  sand  to  sand  barrier  along  the  axis  of  the  footing. 
As  the  sand  in  this  region  is  completely  disturbed,  it  must  be  at  the  minimum  density,  it  may 
be  assumed,  therefore,  that  the  material  is  close  to  the  active  state.  The  value  of  the  angle 
of  interface  friction  probably  lies  between  zero  and  the  value  of  0  at  minimum  density. 
Similarly,  the  value  of  the  shaft  coefficient  of  earth  pressure  is  probably  between  the  active 
value  and  unity. 

It  is  proposed,  therefore,  that  adequate  accuracy  may  be  achieved  from  the  Ny  values 
fora  shaft  coefficient  of  unity  interpolated  for  the  value  of  the  angle  of  interface  friction  equal 
to  the  value  of  0  at  minimum  density. 

2.9.  1.  herent  Errors 

Some  penalty  must  be  paid  in  fulfilling  the  requirements  of  rotational  equilibrium.  This 
jienalty,  resulting  in  some  inherent  errors  in  the  calculation,  may  be  assessed  from  the  varia¬ 
tion  of  the  values  of  the  coefficients  of  earth  pressure  and  the  stress  distribution  along  the 
sides  of  the  wedge  and  along  the  plane  shear  zone  boundaries.  For  simplicity,  the  errors  for  a 
cohcsionless  soil  and  for  a  cohesive  soil  are  analysed  sequentially. 

From  the  initial  assumption  that  the  wedge  boundary  is  straight,  it  is  logical  to  adopt 
the  hydrostatic  distribution  for  the  normal  stresses  dependent  on  the  weight  of  the  displaced 
material.  As  sinkage  increases,  the  centre  of  pressure  for  the  weight-frictional  component 
rapidly  moves  towards  the  mid-point  of  the  wedge  boundary  (Fig  2.9.1.).  The  effective  move¬ 
ment  of  the  centre  of  pressure  is  exhibited  in  the  Ny  curves  whose  gradient  increases  with 
penetration  ratio.  The  assumed  stress  distribution  on  the  wedge  boundary  suffers  from  two 
known  sources  of  error,  particularly  at  the  deeper  sinkages.  First,  the  material  has  inherent 
frictional  properties,  whereas  hydrostatics  is  applicable  to  frictionless  media.  Secondly,  the 
stress  distribution  assumed  for  the  wedge  boundary  differs  from  that  for  the  pair  of  failure 
planes.  In  both  cases,  however,  the  error  is  minimised  because  the  total  length  of  the  wedge 
boundary  is  small  in  relation  to  the  total  depth  at  which  the  discrepancies  betweer,  the  stress 
distributions  become  significant. 

The  same  hydrostatic  distribution  is  used  for  the  weight-frictional  stresses  on  tin- 
foundation  shaft.  This  seems  adequately  justified  bv  the  experimental  data  37,50  on  the  varia¬ 
tion  of  skin  resistance  with  foundation  depth  up  to  a  penetration  ratio  of  10  in  compact  soil. 

Using  a  hydrostatic  distribution  on  the  shaft  makes  it  impossible  to  use  the  same 
distribution  on  the  two  failure  boundaries  of  the  plane  shear  zone  unless  0-0,  The  introduc¬ 
tion  of  an  exponent  to  vary  the  pressure  distribution  along  these  two  boundaries  satisfies  the 
equations  of  equilibrium.  However,  the  curved  distribution  could  not  occur  along  a  straight 
failure  boundary  and  it  is  reasonable  to  assume  that  this  failure  geometry  is  an  approximation 
to  the  curved  boundaries  .hown  in  Fig  2.2.4a. 

The  rather  nebulous  errors  involved  in  the  determination  of  the  boundary  conditions  for 
the  plane  shear  zone  are  best  assessed  from  a  study  of  the  calculated  values  of  the  boundary 
coefficient,  the  boundary  exponent  and  the  slope  angle,  and  an  understanding  of  the  complex 
mechanism  for  achieving  equilibrium.  The  respective  values  of  these  parameters  over  a  range 
of  0  are  plot  ted  for  an  aruitrarv  value  of  the  shaft  coefficient  of  unity  when  the  shaft  is 
peifectly  smooth  and  perfectly  rough  (Figs  2.9.2.  to  2.9.5.).  The  equilibrium  mechanism  may 
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be  considered  in  two  parts:  first,  for  a  constant  slope  angle  and  varying  pressure  distribution, 
and  secondly  for  the  opposite  case.  Taking  a  specific  value  of  the  penetration  ratio,  the 
shaft  force  is  constant  for  any  given  set  of  shaft  parameters.  In  the  first  case,  the  magnitude  and 
and  direction  of  the  force  due  to  the  internal  weight  of  the  soil  mass  remain  constant  and  the 
directions  of  the  boundary  forces  are  defined  by  the  slope  angle.  Under  these  restraints,  the 
magnitudes  and  lines  of  action  of  the  boundary  forces  are  also  defined.  They  may  be  achieved 
by  either  increasing  the  value  of  the  boundary  coefficient  and  decreasing  the  value  of  the 
exponent  or  vice  versa.  The  stress  distribution  for  the  former  case  is  more  unitorm  and  that 
for  the  latter  increases  rapidly  with  depth.  In  general  terms,  the  uniform  distribution  applies  to 
a  smooth  shaft  and  other  to  a  rough  shaft.  In  the  second  case,  the  slope  angle  not  only 
controls  the  direction  of  the  boundary  forces  but  also  the  internal  weight  of  the  plane  shear 
zone  and  its  centre  of  gravity.  Depressing  the  slope  angle  causes  the  boundary  forces  to  move 
towards  the  intersection  of  the  planes  and  alters  the  relative  lengths  of  the  two  boundaries. 

This,  in  turn,  increases  the  internal  weight  of  the  zone  and  changes  the  position  of  the  centre 
of  gravity  such  that  the  stresses  on  the  outer  plane  boundary  are  increased.  Basically, 
depression  of  the  slope  angle  provides  a  larger  horizontal  component  of  the  weight-frictional 
forces  acting  on  the  outer  boundary  of  the  plane  shear  zone  and  explains  why  this  process  is 
involved  when  full  passive  pressure  is  assumed  to  act  on  a  perfectly  smooth  shaft. 

In  soils  with  cohesion,  the  boundary  coefficient  and  exponent  additionally  control  (in 
modified  form)  the  divergence  of  the  cohesive-frictional  stress  distribution  on  the  plane  shear 
zone  boundaries  from  that  described  as  uniform  along  the  shaft  and  wedge  boundaries  so  that 
equilibrium  conditions  may  be  satisfied.  The  variation  of  both  the  boundaiv  coefficient  and 
exponent  with  the  penetration  ratio  is  shown  in  Figs  2.9.6.  and  2.9.7.  for  various  shaft 
conditions  and  c/VR  ratios.  The  effect  of  increasing  the  c/yB  ratio  is  to  reduce  the  value  of 
the  boundary  coefficient  and  to  inciease  the  value  of  the  exponent,  l  or  both  the  unsupported 
shaft  and  the  smooth  shaft,  the  value  of  the  coefficient  tends  towards  unity  and  that  of  the 
exponent  towards  zero.  In  the  case  of  the  rough  shaft,  the  equilibrium  condition  is  more 
difficult  to  achieve  because  of  the  shaft  traction  which  is  generated  by  the  adhesive  forces  in 
addition  to  the  'vertical'  components  of  the  frictional  forces  on  the  shaft.  The  large  shift  in 
stress  distribution  on  the  boundaries  of  t he  plane  shear  zone  may  be  achieved  only  by  high 
values  of  the  exponent  which  in  turn  governs  low  values  of  the  coefficient.  This  explains  why 
equilibrium  was  not  achieved  when  the  c/'Al  ratio  exceeds  unity. 

There  is  no  rigorous  method  of  evaluating  the  errors  for  the  approximate  boundary 
conditions.  In  a  negative  sense,  the  variation  in  the  boundary  parameters  from  those  accepted 
for  zero  sinkage,  indicate  the  undisclosed  errors  of  earlier  theories  in  which  the  combination 
of  hydrostatic  and  uniform  stress  distributions  is  adopted.  Moreover,  the  vaiiation  indicates 
inconsistencies  in  the  present  theory  because  the  equilibrium  of  the  radial  shear  zone  is  still 
not  considered.  A  suitable  extension  of  the  existing  stress  distribution  on  the  plane  boundaries 
to  the  spiral  boundary  would  enable  the  only  remaining  unknown  forces  in  the  system  to  be 
determined.  1'he  equation  of  equilibrium  for  the  radial  shear  zone  could  then  i.o  calculated  and 
used  to  assess  the  existing  errors.  These  errors  couid  be  eliminated  by  changing  the  spiral 
boundary  to  some  form  of  ellipse  whose  shape  could  be  determined  by  the  angle  of  dilation. 

The  adoption  of  a  variable  angle  of  dilation  would  permit  variation  in  density  to  be  reflected 
in  the  shape  of  the  failure  boundary.  ’Ihe.se  improvements,  however,  could  lie  achieved  only 
at  the  expense  of  additional  theoretical  complexity. 


2.10.  Comparison  with  Other  Theories 

The  theory  of  penetration  resistance  is  compared  with  existing  theories  of  hearing 
capacity  for  both  a  foundation  resting  on  the  soil  surface  and  for  a  buried  foundation  in 
cohesionlcss  and  tn  cohesive  soils. 

The  values  of  the  surface  hearing  factor,  N'v,  are  shown  in  Fig  2.10.1.  1  erznghi's 
and  Meyerhof's  N-,  factors  wore  computed  by  the  author;  the  former  are  exactly  the  same  as 
those  published  by  Lotkin53  but  the  latter  differ  slightly  from  those  published  in  reference  2d, 
particularly  for  0  ~  5°,  10°,  15°,  and  50°.  Balias'  values25  are  derived  with  the  aid  of  his 
graphs  for  the  ratio  r,  and  Mizuno's  values54  are  obtained  from  his  published  graphs. 
Published  values  are  used  to  construct  the  curves  of  Caquot  and  Kcrisel’s55  and  l.undgren 
and  Mortrnsen's 1  8  Ny  factors.  As  the  current  theory  is  based  on  the  maximum  size  of  the 
wedge  with  no  minimising  procedures,  it  is  not  surprising  that  the  Ny  values  tend  to  be 
higher  than  most  other  theories. 

The  comparison  at  shallow  depths  of  penetration  is  complicated  by  the  fact  that  none 
of  the  bearing  capacity  theories  consider  shaft  support  in  a  logical  way.  Terzaghi's  theory  of 
surface  hearing  capacity  may  not  he  extrapolated  beyond  a  penetration  ratio  of  1  because  the 
soil  above  the  base  of  the  foundation  is  considered  as  a  pure  surcharge  with  no  boundary 
stresses.  Balias’  theory  is  restricted  to  a  penetration  ratio  of  1.5  because  further  extension 
would  lead  to  an  impossible  figure  of  rupture  when  shaft  support  is  ignored.  Meyerhof 
considers  shaft  support  as  a  surcharge  effect  for  a  figure  of  rupture  dissimilar  to  that  for  the 
determination  of  the  bearing  capacity  due  to  the  internal  weight  of  the  mass  of  soil  at  failure. 

The  dimensionless  graphs  p/liTB  and  p/c  are  used  in  Figs  2.10.2.  and  2.10.3  to 
account  for  the  combined  effects  of  Ny,  \q  and  Ny,  Nq,  Nc,  respectively.  The  excessively 
small  scale  of  Meyerhof's  published  graphs16  of  the  combined  factors  in  the  relevant  range 
precludes  a  direct  comparison  with  his  work.  For  expediency,  Wills  results32,  for  a 
restrained  version  of  Meyerhof's  theory,  are  presented  for  in  =  1.  From  a  few  test  computa¬ 
tions  of  Meyerhof's  theory,  it  would  appear  that  the  effect  of  restraint  on  the  pole  position  is 
to  increase  the  Nyq  factor  from  a  value  less  than  the  current  theory  to  a  value  greater  than 
the  current  theory  when  the  shaft  coefficient  of  earth  pressure  is  unity  and  the  interface 
friction  is  zero.  (The  Nyq  factor  is  not  affected  hy  restraint). 

Although  direct  equivalence  of  shaft  conditions  is  not  possible  in  the  comparison,  the 
present  theoretical  curve  exhibits  a  marked  similarity  to  the  existing  bearing  capacities 
theories  adopting  some  form  of  spiral  failure  boundary.  The  divergence  of  Balias’  results 
from  the  general  trend  is  largely  due  to  the  use  of  a  failure  boundary  which  does  not  appear 
to  have  experimental  validity. 
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3.  APPARATUS,  SOIL  AND  EXPERIMENTAL  PROCEDURE 


The  experimental  verification  of  a  pressure  penetration  theory  requires  apparatus  for 
the  measurement  of  the  forces  on  a  rigid  model  foundation  as  it  is  pushed  vertically  into  the 
soil.  The  evaluation  of  two  dimensional  soil  movement  necessitates  the  use  of  a  glass-sided 
tank,  and  in  order  to  cover  the  complete  soil  range,  tests  must  be  conducted  in  at  least  three 
soils,  a  cohesionless  soil,  one  with  cohesion  and  friction  and  a  frictionless  material.  Each 
soil  must  be  processed  to  a  stale  which  can  bo  readily  reproduced  for  a  number  of  experiments, 
and  the  soil  parameters  must  be  obtained  for  these  precise  conditions.  Finally,  the  interface 
effects  due  to  model  testing  in  a  confined  space,  the  influence  of  the  glass  interfaces  and 
sources  of  error  incurred  by  the  photographic  analysis  of  the  failure  geometry  must  be 
investigated. 

3.1.  The  Pressure  Penetration  Apparatus 

Pressure  penetration  experiments  were  conducted  on  a  large  penetrometer  developed  by 
U.S.A.T.A  C.  (Fig  3.1.1.).  This  was  used  in  order  to  obtain  the  structural  rigidity  necessary 
for  the  development  of  symmetrical  soil  failure.  The  penetration  device  is  mounted  on  the  base 
of  a  vertical  shaft  running  through  linear  bearings.  This  prevents  any  lateral  movement  of  the 
shaft.  Although  the  shaft  is  not  completely  restrained  in  rotational  movement,  it  did  not  appear 
to  occur  during  the  tests.  1  he  vertical  shaft  movement  of  approximately  42  in  is  controlled 
by  two  double  acting  hydraulic  cylinders.  The  penetration  speed  may  be  varied  from  12-60  in  min'1 
by  means  of  a  variable  flow  hydiaulic  control  unit. 

A  glass  sided  soil  tank  was  mounted  cm  a  stand  to  simplify  the  photographic  procedures. 
The  maximum  size  of  the  tank  was  limited  by  the  weight  which  could  be  conveniently  man¬ 
handled,  and  by  the  volume  of  soil  which  would  be  readily  processed  between  experiments.  The 
internal  dimensions  of  the  tank  were  45  in  long,  35  in  high  and  6  in  wide  and  it  contained 
approximately  500  lb  of  soil  to  a  total  depth  of  30  in.  As  the  extent  of  failure  is  a  function  of  0, 
provision  was  made  for  varying  the  quantity  of  soil  by  using  two  vertical  compression  plates 
which  could  be  moved  inwards  to  give  any  desired  tank  volume. 

The  base  and  ends  of  the  tank  were  of  6  in  x  2  in  channel.  A  !;  in  thick  ‘Plexiglass’ 
sheet  was  used  to  form  the  front  of  the  tank  and  a  h  in  thick  sheet  of  plywood  formed  the  back. 
The  ends  were  braced  along  the  length  of  the  tank  and  the  sides  across  the  width  of  the  tank 
by  a  series  of  removable  cross  bars  (big  3.1.2.).  (It  was  subsequently  found  that  the  2  in  x  i:  in 
crossbars  had  to  be  strengthened  by  a  factor  of  3  to  minimise  elastic  deformation  of  the  cross¬ 
bars  during  tests  in  very  compact  materials).  On  the  inside  faces  of  the  ‘Plexiglass’  front  and 
the  wooden  back  were  placed  two  sheets  of  V8  in  thick  window  glass  to  form  the  soil  glass 
interfaces.  The  inner  sheets  were  used  to  reduce  replacement  costs;  these  were  expendable 
when  scoring  occurred.  For  a  free  flowing  material  such  as  dry  sand,  the  replacement  of  the 
glass  was  only  possible  when  the  box  was  empty.  For  self  supporting  cohesive  materials,  the 
glass  could  be  replaced  without  emptying  the  soil  by  gently  moving  the  cross  bars  after 
processing. 

The  penetration  device  was  designed  to  measure  both  the  base  resistance  and  the  shaft 
forces.  The  outer  shell  of  the  penetration  device  was  a  channel  section  measuring  30  in  Jong, 

6  in  deep  and  1*  in  wide  (Fig  3.1.3.).  This  housed  the  loading  shaft  which  had  a  60001 L 
integral  ring  dynamometer  at  the  lower  end  to  measure  the  base  resistance.  The  horizontal  and 
vertical  side  plate  loads  were  transmitted  to  five  cantilevers  on  the  loading  shaft  through 
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spherical  rod  end  bearings  and  small  chain  links.  Any  lateral  movement  of  the  side  plate  was 
constrained  by  ball  bearings  recessed  in  grooves  in  the  strengthening  bars  on  the  side  plate 
and  impinging  on  the  channel  flanges.  By  measuring  the  horizontal  loads  with  two  pairs  of 
cantilevers,  each  with  a  design  load  of  100  lb,  the  total  thrust  and  its  point  of  application 
could  be  determined.  The  vertical  thrust  was  measured  by  a  single  cantilever  with  a  design 
load  of  300  lb. 

t 

The  base  resistance  dynamometer  was  calibrated  against  a  standard  load  cell 
(Fig  3.1 .4.5.  A  linear  relationship  between  strain  rate  and  load  was  obtained  with  no 
hysteresis  effects. 

The  side  plate  cantilevers  were  calibrated  by  dead  weights.  As  the  individual  calibra¬ 
tions  of  the  four  cantilevers  measuring  the  horizontal  thrust  were  very  uniform,  the  two  signals 
from  the  two  lower  cantilevers,  and  those  from  the  two  upper  cantilevers  were  added 
electrically  to  give  only  two  outputs  instead  of  four.  On  assembling  the  complete  side  plate, 
however,  problems  were  experienced  with  the  calibration  of  the  cantilevers.  It  was  anticipated 
that  some  interaction  would  occur  between  the  horizontal  and  vertical  force  measuring  units, 
but  unforunately  this  was  combined  with  excessive  hysteresis.  This  was  attributed  to  some 
major  misalignments  in  the  assembly  of  the  cantilever  linkage.  Although  some  modification 
did  reduce  the  friction,  there  was  insufficient  time  to  redesign  and  manufacture  a  linkage  to 
eliminate  the  problem  and  the  existing  rig  had  to  be  used  for  the  current  series  of  experiments. 
Owing  to  these  factors,  a  comprehensive  series  of  calibrations  were  conducted  on  the  side 
plate.  The  penetration  device  was  mounted  in  a  Hartford  chuck  and  measurements  were  made 
for  a  series  of  weights  placed  at  different  points  along  the  length  of  the  side  [date  and  three 
different  angles  of  the  side  plate—  horizontal,  45°  and  vertical— were  also  used  to  assess  the 
interaction  of  combined  loading  (Fig  3,1.5.). 

The  cantilever  calibration  curves  are  shown  in  Figs  3.1.6.  to  3.1.0.  The  results  for 
horizontal  loading  do  not  exhibit  too  much  hysteresis,  and  by  combining  the  upper  and  lower 
calibration  curves,  the  resultant  calibration  curve  is  identical  and  independent  of  the  position 
of  the  load  (Fig  3.1.7.),  the  interaction  obtained  in  horizontal  loading  was  found  to  be 
relatively  constant  at  9T  of  the  horizontal  load  (Fig  3.1.8.).  This  high  figure  is  the  result  of 
minimising  hysteresis  by  increasing  the  linkage  clearances  which  then  introduced  considerable 
slack.  An  angular  movement  of  only  t  4°  from  the  vertical  could  account  for  this  interact  inn. 
This  slack  also  caused  ‘weight  transfer'  to  occur  between  the  upper  and  lower  cantilevers  when 
when  the  side  plate  was  tested  at  an  angle  to  the  horizontal.  It  was  therefore  decided  that 
the  total  force  calibrations  could  be  used  but  that  the  location  of  the  point  of  application  of 
the  total  force  was  subject  to  an  excessive  degree  of  error  The  horizontal  component  of  the 
force  applied  at  45°  was  identical  to  that  found  for  horizontal  loading;  the  vertical  component 
was  some  8fc  higher  which  agrees  well  with  the  interaction  figure  obtained  earlier  (Fig  3.1.9.). 

The  base  resistance  was  recorded  against  sinkage  on  an  XY  plotter  and  the  three 
force  measurements  on  the  side  plate  were  also  recorded  against  sinkage  on  a  multi-pen 
recorder.  The  two  sinkage  measurements  were  each  obtained  from  hall  bridge  incorporating  a 
rotary  potentiometer  which  was  forced  to  rotate  when  penetration  occurred. 

Two  further  refinements  were  made  to  the  basic  rig.  Rough  or  smooth  side  plates  and 
ba  sc  plates  could  be  fitted  to  investigate  the  effect  of  interface  friction.  The  aspect  ratio  of 
the  penetration  device  may  be  increased  from  3:1  to  12:1  by  attaching  *  wooden  block  on  each 
side  of  the  central  measuring  unit  (Fig  3.1.10).  This  was  used  to  determine  the  effect  of 
the  glass  interfaces. 
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Some  of  the  experimental  results  used  in  this  report  were  obtained  from  an  earlier 
series  of  tests  which  were  conducted  on  a  different  apparatus  using  footings  without  side- 
walls.  This  apparatus  is  fully  described  in  reference  40  together  with  the  photographic 
tecliniques  and  the  methods  employed  to  reduce  the  effect  of  the  glass  interfaces. 

3.2  Soil  Preparation  and  Properties 

Tests  were  conducted  in  a  roller  model  soil,  two  dry  sands  (Oh io  and  Leighton  Buzzard), 
damp  sand  and  clay.  The  rollers  were  processed  manually  by  rearranging  their  positions  until 
any  visible  failure  lines  from  the  previous  tests  had  b  en  removed  and  some  degree  of  compac* 
tion  had  been  achieved  throughout  the  mass.  Care  was  also  taken  to  ensure  that  completely 
uniform  packing  did  not  occur. 

The  dry  sand  was  first  completely  cultivated  by  emptying  and  refilling  the  tank,  then 
compacted  by  an  electro-magnetic  vibrator  until  maximum  density  was  achieved  and  the 
surface  levelled.  In  the  larger  tank,  a  vibrating  probe  was  also  used  to  ensure  even  compac¬ 
tion. 


The  loose  damp  sand  was  spread  evenly  in  2  in  layers  and  compacted  by  means  of  a 
falling  weight  compactor  to  maximum  density.  The  surface  of  each  successive  layer  was 
disturbed  before  the  addition  of  fresh  material  and  final  layer  levelled  to  a  set  depth  by  a 
scraper  running  along  the  sides  of  the  tank. 

Whereas  the  tank  was  completely  emptied  of  sand  between  each  test,  the  clay  was  not 
disturbed.  It  was  processed  by  filling  in  the  cavity  left  by  the  fooling  and  kneading  the  mass 
to  obliterate  any  rupture  lines.  The  surface  was  levelled  to  a  set  depth  using  the  same  volume 
of  clay  as  a  check  for  the  presence  of  any  air  pockets. 

The  diameters  of  the  wooden  rollers  were  5/8  in  and  %  in.  At  a  mean  density  of 
0,02  lb/in'3  the  angle  of  internal  friction  was  28°. 

The  Ohio  sand  contained  rounded  particles  with  a  good  grain  size  distribution 
(Fig  3.2.1.).  The  angle  of  internal  friction  was  40.5°  in  a  density  of  0.065  lb  in'3.  The  angle 
of  repose  was  32°  at  a  minimum  density  of  0.59  lb  in'3  . 

The  Leighton  Buzzard  sand  was  a  sieved  washed  sand  again  with  rounded  particles 
but  with  a  very  uniform  large  grain  size  laying  between  12  mesh  (0.0553  in)  and  24  mesh 
(0.0236).  At  a  maximum  density  cf  0.0605  lb  in'3  the  angle  of  internal  friction  was  45.5°.  At 
the  minimum  density  of  0.054  lb  in'3,  0  was  32°.  Using  these  two  values  of  the  density,  the 
maximum  dilation  is  10%. 

The  damp  sand  was  a  fine,  washed  ‘plastering’  sand  with  sharp  grains  and  had  a 
moisture  content  of  4.53%.  At  a  maximum  density  of  0,0672  11)  in'3  the  soil  properties  were 
(/  =  43.4°  and  c  ::  0.3  lb  in*3.  At  minimum  density  of  0.047  lb  in'3,  0  “  33°  and 
c  -  0.2  lb  in"2. 

The  clay  had  a  cohesion  of  1.8  lb  in*^  and  a  negligible  friction  angle  of  between  2.5° 
and  4°.  The  moisture  content  was  32.8%  (D.B.)  and  the  density  0.0625  lb  in'3. 
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4.  ANALYSIS  OF  RESULTS 


The  measured  pressures  as  a  function  of  the  penetration  ratio  gave  good  agreement 
with  the  theory  in  ail  materials.  The  photographs  showing  the  failure  boundaries  in  the  roller 
model  were  almost  identical  with  the  theory  and  those  in  dr/ and  damp  sand  were  encouraging 
in  that  they  were  of  the  same  fomr  as  the  predicted  lioundaries.  Unfortunately,  the  base  wedge 
was  much  flattened  and  the  overall  width  of  failure  in  sand  was  only  about  hall  as  big  as  the 
theory  predicts.  In  the  clay,  the  failure  was  markedly  different  to  that  predicted. 

4.1  I Vnctra t ion  F ailuro  in  the  Holler  Model 

The  theoretical  determination  of  the  base  resistance  and  the  failure  geometry  depends 
on  the  coefficient  of  earth  pressure  and  the  angle  of  interface  friction  on  the  shaft  which  wen' 
found  by  experiment  to  be  0.45  and  12°  ,  respectively.  The  experimental  jx'itetration 
resistance  together  with  the  predicted  curve  for  kys  =  0.45  and  S  r  12°  is  shown  in 
big  4.1.1. 

The  failure  geometry  is  shown  in  a  series  of  photographs  of  a  2  in  wide  sand  coated 
foundation  driven  progressively  into  the  roller  model  from  the  surface  to  a  penetration  ratio 
of  10  (1  igs  4.1.2a.  to  4.1.2].).  Mobilisation  of  the  material  did  not  always  occur  on  both 
sides  of  the  foundation  shaft  at  the  same  instant  even  though  great  care  was  taken  to  ensure 
that  the  penetration  device  was  vertical.  A  particularly  good  example  of  symmetrical  failure 
is  shown  in  1  ig  4.1 .2j.  After  a  number  of  tests,  the  boundary  of  the  failure  /.one  was  clearly 
visible  by  still  photography  (Fig  3.1.1.).  This  was  made  possible  by  the  mixing  action  of 
successive  penetrations  within  the  mobilised  material.  Rollers  were  carried  down  in  the-  \w  dec 
attached  to  the  base  of  the  footing.  (Figs  4.1.3.).  The  base  wedge  angle  varied  between  57" 
and  60"  beneath  a  rough  footing.  However,  it  is  not  possible  to  lay  too  much  stress  on  tlmse 
figures  because  there  was  a  tendency  for  the  rollers  to  be  arranged  in  the  ideal  packing  forma¬ 
tion  by  displacing  the  larger  rollers  (l  ie  4.1.4.).  This  results  in  a  base  wedge  ancle  of  (>(i" 
which  does  not  necessarily  reflect  on  the  actual  value  of  0  or  the  correct  stress  field  when 
both  sizes  of  rollers  are  present. 

A  comparison  between  the  experimental  and  theoretical  failure  geometry  is. shown  in 
Fig  4.1.6.  for  three  penetration  ratios.  Considering  the  experimental  variation  between  tests, 
the  relationships  are  very  satisfactory.  At  the  largest  penetration  ratio,  the  photograph  of  lire 
roller  movement  (big  4.1.2i.)  clearly  shows  two  roller  streams,  one  moving  in  a  direction 
similar  to  that  predicted,  and  the  other  moving  outwards  and  upwards  following  the  line  of 
least  resistance  for  uniformly  packed  rollers.  The  movement  of  both  streams  is  largely 
dissipated  before  they  reach  the  surface. 

4.2  Failure  under  Hough  Strip  Foundations  and  Footings  in  l)rv  Sand 

Preparatory  tests  with  the  penetration  device  were  conducted  in  a  large  mass  of  dry 
sand  to  determine  the  effect  of  density  on  the  base  resistance  and  on  the  appearance  of 
rupture  planes  at  the  sand  surface.  The  influence  is  demonstrated  by  the  two  pressure  penetra¬ 
tion  records  for  Ohio  Sand  (Fig  4.2.1.).  In  soil  at  the  minimum  density  the  trace  is  a  smooth 
curve  which.is  the  result  of  localised  compaction  failure  within  the  soil  mass.  Almost  no 
surface  upheaval  occurred  at  the  larger  penetration  ratios.  When  the  compressibility  effect  is 
eliminated,  the  experimental  penetration  resistance  is  much  less  uniform  with  each  successive 
peak  relating  to  the  development  of  a  clearly  defined  failure  boundary.  The  pressure  ttouglis 
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represent  the  continued  failure  along  an  existing  failure  plane.  Thus  any  predicted  curve 
should  pass  through  the  pressure  peaks.  The  theoretical  penetration  resistance  in  Fig  4.2.1. 
is  based  on  a  measured  coefficient  of  earth  pressure  on  the  shaft  of  1.23  and  S  was  assumed 
to  be  equal  to  0  at  the  minimum  density  of  the  sand. 

The  base  penetration  resistance  for  a  long  strip  footing  (aspect  ratio  of  6:1)  in  compact 
Leighton  Buzzard  sand  was  used  as  a  control  test  to  determine  the  effect  of  the  glass  in tc r- 
face.  As  the  shaft  parameters  for  a  footing  are  indeterminate,  the  penetration  resistance  is 
predicted  in  accordance  with  the  assumptions  given  in  Section  2.8.  The  experimental  and 
theoretical  curves  shown  in  Fig  4.2.2.  exhibit  a  close  similarity.  The  comparison  between 
these  results  and  those  shown  in  Fig  4.3.2.  indicates  that  two  dimensional  failure  may  be 
simulated  in  a  glass  sided  soil  tank. 

4.3  Failure  under  Rough  Footings  in  Dry  Sand 

A  series  of  photographs  record  the  soil  movement  which  was  caused  by  a  2  in  wide 
sand  coated  footing  driven  progressively  ;nto  dry  sand  from  the  surface  to  a  penetration 
ratio  of  3  (Figs  4.3.1a.  to  4. 3. If.).  The  location  of  each  of  the  succession  of  photographs, 
together  with  the  theoretical  base  resistance  is  shown  on  the  XV  plotter  record  for  this 
particular  experiment  (Fig  4.3.2.).  Photographs  (a)  and  (b)  show  the  initial  development  of  a 
long  sharp  wedge  and  sand  particle  movement  extending  very  widely  on  both  sides  of  the 
footing.  Photographs  (c)  and  (d)  are  taken  just  after  the  first  soil  failure  and  shew  a  much 
blunter  wedge  and  smaller  zone  of  moving  sand.  Photograph  (c-)  again  shows  a  sharper  wedge 
and  a  wider  zone.  Notice  that  this  coincides  with  the  second  pressure  peak  on  the  XV  plotter 
record.  1  his  cycle  is  repeated  throughout  the  test  and  is  accompanied  by  a  general  widening 
of  the  failure  zone  with  increasing  depth  as  shown  in  photograph  (f).  The  existence  and  the 
variation  in  the  size  of  t he  wedge  is  shown  more  clearly  in  photographs  taken  with  the  camera 
attached  to  the  footing  (Figs  4.3.3a.  and  4.3.3b.).  The  maximum  pressure  tends  to  occur 
when  the  wedge  is  at  its  sharpest  and  the  peak  values  correspond  well  with  the  theory. 

The  vaiiatic-n  in  wedge  shape  and  the  smaller  failure  boundaries  are  attributed  to  tin- 
variation  in  0  with  changes  in  density.  The  dilation  of  the  sand  due  to  the  first  footing  failure 
was  measured  with  a  planimeter  from  photographs  4.3.1b.  amd  4. 3. Id.  In  both  cases  the 
volume  of  sand  lias  increased  and  corresponds  to  an  V/.i  expansion  of  the  sand  lying  above 
the  soil  failure  Loundary.  Photograph..  4.3. id.  also  shows  two  edges  to  the  surface  heave. 

The  inner  edge  corresponds  with  the  intersection  between  the  main  failure  boundary  and  the 
surface  of  the  sand.  The  outer  edge  is  much  further  out  and  represents  the  much  large  volume 
of  sand  mobilised  in.  photographs  4.3.1a.  and  4.3.1b. 

As  the  footing  begins  to  sink,  the  sand  deforms.  Since  it  is  originally  at  its  maximum 
density,  it  dilates  and  the  value  of  0  diminishes.  This  process  tends  to  spread  from  the 
edges  of  the  footing,  and  before  the  complete  failure  zone  develops,  the  footing  is  operating 
in  a  sand  with  0  =  32°  instead  of  40°.  For  such  a  material,  the  failure  surface  is  only  half 
as  large.  A  secondary  failure  therefore  occurs  entirely  within  the  soil  loosened  bv  the  initial 
penetration  of  the  footing.  It  is  shown  in  Fig  4.3.4.  that  the  measured  failure  surfaces 
correspond  very  closely  to  the  theoretical  one  for  0  =  32°.  This  comparison  is  based  not 

on  the  same  penetration  ratio  for  theory  and  experiment,  but  on  the  same  position  of  the  wedge 
tip.  As  the  original  wxdge  tip  constitutes  a  pattern  of  weakness,  it  tends  to  remain  in  its 
initial  position.  The  wedge  angle  is  reduced  within  the  loosened  material  by  the  sinkage  of 
the  footing  which  is  required  to  produce  the  first  failure. 
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With  further  penetration,  the  footing  nears  the  original  wedge  lip  and  forms  a  now  sharp 
wedge.  Front  this  wedge,  the  soil  is  mobilised  outwards  and  upwards  to  start  a  new  cycle  of 
soil  loosening,  weakening  and  failing. 

The  magnitude  of  the  recorded  penetration  resistance  confirms  this  explanation.  It  is 
only  thi'  peak  values  of  the  pressure  that  equal  the.  theory.  In  between,  the  pressure  falls 
towards  the  value  appropriate  to  the  lower  value  of  0.  There  is  a  tendency  for  only  the  first 
peak  to  reach  tin*  predicted  vaine  and  all  the  subsequent  peaks  to  be  lower.  This  is  because 
the  subsequent  failures  have  previously  loosened  sand  txjbove  the  newly  generated  boundary. 


4.4  Failure  under  Smooth  Footings  in  Dry  Sand 


Photographs  4.4.1a.  to  4. 4. Id.  are  of  a  lubricated  glass  fooling  taken  with  a  camera 
rigidly  mounted  to  the  footing  and  show  the  fixed  wedge  and,  less  clearly,  the  failure  zone 
boundary.  The  experimental  penetration  resistance  is  shown  in  Fig  4.4.2.  The  photographs 
(and  others  which  were  taken)  show  that  the  failure  zones  are  of  the  same  magnitude  as  those 
foi  the  rough  footing  and  have  a  wedge  behaving  in  the  same  cyclic  manner.  The  only 
difference  is  that  the  wedge  can  disappear  altogether  and,  on  occasion,  a  double  wo  dee  can 
formed.  This  latter  phenomenon!  together  with  the  maximum  size  of  the  wedge  at  the  surface 
arc  illustrated  in  b  ig  4.4.3. 


The  cycle  variation  of  the  wedge  beneath  a  rough  and  smooth  footing  is  shown 
diagrammatical  ly  in  bigs  4.4.4a.  and  4.4.4b.  Note  that  the  tip  of  the  original  wedge,  resting 
on  the  compact  soil,  remains  constant  througiiout  one  cycle.  The  minimum  size  of  the  wedge 
depends  on  the  mobilisation  of  shear  stresses  on  die  base  of  the  footing.  Beneath  a  rough 
fooling,  the  wedge  never  completely  disappears  (2nd  stage  of  Fig  4.4.4a.).  A  larger  wrtlire 
then  forms  to  start  a  new  cycle  but  it  never  reaches  the  maximum  size  of  the  wedge  present 
during  the  initial  process  because  the  upper  half  of  the  new  wedge  is  developing  in  previously 
loosened  material. 


Beneath,  a  smooth  footing,  no  shear  stresses  can  be  developed  and  the  complete  cych- 
includes  two  further  stages.  In  some  cases,  a  mere  vestige  of  the  wedge  is  retained  t3rd  stage 
of  Fig  4.4.4b.).  Finally,  the  footing  may  penetrate  until  it  reaches  the  original  curved 
boundaries  emanating  from  th"  tip  to  form  two  planes  of  weakness  at  an  angle  of  r.; 2  +  Jl/2 

to  each  other  joining  at  the  c  ntre  of  the  base.  Thus,  they  form  the  inner  boundaries  of  two 
wedges  such,  that  the  soil  flows  horizontally  outwards  from  the  centre  of  the  footing  (4th  stage 
of  F ig  4.4.4b.).  Thereafter,  the  cycle  is  repeated.  The  new  wedge  is  almost  as  large  as  the 
original  because  it  is  being  developed  in  almost  totally  undisturbed  material.  This  explains 
why  all  the  peaks  in  the  pressure  penetration  record  are  of  the  same  amplitude  as  the  first 
peak  and  why  they  occur  less  frequently  than  in  the  record  for  the  rough  footing. 


There  is  no  theoretical  explanation  for  the  cyclic  variation  of  wedge  chape  based  on 
the  assumption  of  a  rigid  footing  penetrating  a  homogeneous  material.  The  observed  tendency 
for  assymmetrical  failure  after  the  initial  mobilisation  of  the  soil  imposes  considerable  lateral 
thrust  on  the  wedge  and  can  cause  its  partial  destruction  (Fig  4.4.5.).  The  mode  and  extent 
of  the  destruction  is  then  determined  by  the  interface  properties  along  the  base  of  the  footing. 
In  the  presence  of  interface  friction,  the  lateral  thrust  is  imparted  to  the  footing,  causing  it  to 
‘weave’  during  penetration.  Slices  of  the  wedge  are  removed  by  a  shearing  action  until  the 
ba  se  angle  readies  the  value  of  0.  At  this  point,  the  lateral  thrust  is  eliminated  because  the 
resultant  forces,  acting  at  an  angle  of  0  to  the  normal  on  the  sides  of  the  wedge,  are  vertical 
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and  maximum  shear  stresses  can  be  developed  along  the  base  to  prevent  a  further  reduction 
in  the  wedge  size.  The  footing  must  now  displace  the  dense  sand  beneath  the  blunt  wedge. 
Applying  the  principle  of  conservation  of  energy,  the  sharp  wedge  is  formed  when  the  penetra¬ 
tion  pressure  for  the  development  cf  a  new  failure  boundary  through  the  undisturbed  soil  is 
less  than  that  for  the  blunt  wedge.  The  process  is  similar  for  a  smooth  footing  except  that 
the.  absence  of  shear  stresses  permits  the  complete  disappearance  of  the  wedge  without 
imposing  any  lateral  thrust  on  the  footing  itself. 

4.5  Penetration  Failure  in  Damp  Sand 


Preliminary  experiments  were  carried  out  to  assess  the  difference  between  a  rough 
and  smooth  footings.  No  difference  in  either  the  failure  patiems  or  the  pressures  could  be 
detected.  This  may  well  be  attributed  to  the  difficulty  of  eliminating  soil  to  structure 
adhesion.  The  wedge  beneath  a  lubricated  glass  footing  (  big  4.5.1.)  may  be  compared  with 
the  wedges  for  a  sand  coated  base  (f  ig  4.5.7a.). 

The  main  experiments  were  with  buried  footings  at  initial  penetration  ratios  of  0,  1, 

2  and  3.  The  fooiir.g  was  driven  until  complete  failure  was  observed  in  each  case.  The 
failure  patterns  arc  shown  in  rig  4.5.2.  to  4.5.5. 

Failure  always  occurred  on  one  side  of  the  footing  first,  with  a  relatively  sharp  wedge 
(Figs  4.5.2.,  4.5.4a.  and  4.5.5a.).  Further  sinkage  was  necessary  to  produce  failure  on  the 
other  side,  and  by  this  time  the  wedge  was  smaller  (Figs  4.5.3.,  4.5.4b.  and  4.5.5b.).  The 
photographs  show  the  same  general  patterns  as  in  dry  sand.  The  failure  surfaces  are  of  the 
correct  form  but  are  much  smaller  than  the  theory  predicts.  The  predicted  and  experimental 
failure  patterns  are  compared  for  penetration  ratios  of  0  and  2  in  Fig  4.5.6. 

The  development  of  the  wedge  beneath  a  rough  footing  is  illustrated  in  Fig  4.5.7. 
Initially,  the  theoretical  sharply  pointed  wedge  develops  but  later  there  is  a  tendency  to  form 
a  double  wedge  on  which  appear  the  hall  sized  ilili  type  failure  surfaces. 

The  soil  failure  seemed  to  be  different  to  that  occurring  in  dry  sand  where  there  was  a 
general  disturbance  throughout  the  soil  progressing  from  the  corner  of  the  footing.  In  damp 
sand,  it  appeared  as  if  large  masses  moved  bodily.  It  is  clear  that  large  cracks  could  open  no 
to  separate  the  various  parts  of  the  mobilised  sand  at  laiiun  .  Furthermore,  the  first  signs  of 
failure  occurred  in  the  deepest  layer,  but  before  they  were  fully  developed,  much  smaller 
failure  surfaces  were  generated  at  the  mid  point  of  the  original  wedge  and  Lroke  to  the  surface 
first. 


The  successive  development  of  the  failure  surfaces  for  a  driven  footing  are  shown  in 
Fig.  4.5.8.  In  this  case,  there  is  an  obvious  tendency  for  the  successive  failures  to  rise 
sharply  upwards  and  coalesce  with  previous  failures,  thus  producing  a  very  smell  rupture 
distance. 


The  force  measurements  for  both  buried  and  driven  footings  are  shown  in  Fig  4.5.9. 
together  with  the  predicted  values.  Once  again,  there  is  good  agreement  considering  that  the 
interface  effects  were  higher  in  this  material.  (The  adhesion  and  interface  friction  were  not 
reduced  by  the  lubricant  to  the  same  degree  as  in  dry  sand.) 

Once  again,  the  measured  value  of  0  for  the  sand  in  its  dense  state  gives  the  correct 
theoietical  value  for  the  base  resistance  but  the  faiha*  surface  i:  r-jch  smaller  .red 
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corresponds  to  the  lower  value  of  0.  The  explanation  used  for  the  dry  sand  does  not  seem 
applicable  here.  The  photographs  do  not  exhibit  much  evidence  of  general  dilation  throughout 
the  soil  mass.  It  seems  more  likely  that  the  frictional  resistance  was  reduced  and  the  cohesion 
almost  eliminated  only  along  the  actual  lines  of  rupture.  The  problem  of  soil  compressibility 
was  only  overcome  at  the  expense  of  over-consolidation  and  the  resultant  material  exhibited 
some  recovery  after  deformation.  This,  combined  with  the  evidence  that  parts  of  the  soil  were 
in  tension,  could  yield  the  observed  failure  geometry  and  the  rapid  reduction  in  the  penetration 
resistance  after  the  ctitical  point  of  failure. 

A  further  noteworthy  observation  was  the  marked  tendency  to  unilateral  soil  failure. 
Perhaps  this  situation  should  be  investigated  theoretically.  It  seems  possible  that  an 
alternating  succession  of  one  sided  failures  is  kinematically  feasible  and  might  lead  to  a 
reduced  size  of  failure  pattern. 

4.6  Penetration  Failure  in  Clay 

Entirely  different  trends  were  noted  in  clay.  The  pressure  penetration  traces  for  buried 
footings  driven  from  penetration  ratios  of  0,  1,  2  and  3  were  all  smooth  curves  (Fig  4.6.1.). 

In  each  case,  the  gradient  of  the  curves  rapidly  reduced  to  a  small  constant  value  after  the 
footing  displacement  exceeded  a  distance  equal  to  the  width  of  the  footing.  An  example  of  the 
total  deformation  of  the  soil  during  this  footing  displacement  is  illustrated  in  the  superposed 
photograph  (Fig  4.6.2.).  The  magnitude  and  direction  of  the  actual  soil  deformations  were 
obtained  from  the  two  images  of  each  grid  intersection  and  presented  in  a  displacement 
diagram  (Fig  4.6.3.).  Note  the  vertical  displacement  of  the  soil  together  with  the  complete 
absence  of  any  bulging  of  the  unsupported  walls  of  the  excavation.  It  was  also  found  that  the 
extent  of  failure  increased  in  proportion  to  the  depth  of  penetration. 

This  type  of  failure  was  rather  unexpected.  It  was  anticipated  that  the  soil,  which  was 
displaced  by  the  passage  of  the  footing,  would  flow  into  the  excavation  and  rapidly  yield  some 
constant  maximum  value  for  the  base  resistance  independent  of  depth  in  accordance  with 
Jaky’s  theory8.  Although  Jaky’s  theory  is  for  c/7B  =  00 ,  it  does  compare  well  with  the 
experimental  base  resistance  for  c/7B  =  15,  but  not  with  the  failure  geometry.  The 
observed  failure  patterns  are  attributed  to  two  factors,  the  effect  of  adhesion  and  the  effect  of 
weight.  First,  the  failure  boundaries  shown  in  Fig  4.6.3.  can  be  reproduced  theoretically  in 
the  absence  of  weight  only  when  full  passive  earth  pressure  is  acting  along  the  walls  of  the 
excavation.  Although  no  side  thrust  was  used  in  the  experiments,  a  small  calculation  suffices 
to  show  that  the  adhesion  of  the  glass  interfaces  is  more  than  adequate  to  create  the  necessary 
side  thrust  for  this  type  of  failure.  The  introduction  of  gravitational  stresses  also  tends  to 
increase  the  size  of  the  zone  of  mobilised  soil.  Although  some  tentative  conclusions  have  be 
been  drawn,  confirmation  of  failure  mechanism  in  clay  depends  on  the  complete  elimination  of 
the  adhesion  by  a  suitable  interface  film  which  controls  the  adhesion  without  also  affecting 
the  cohesion. 

4.7  Comparison  with  Other  Experimental  Research 

The  current  theory  tends  to  underestimate  the  surface  bearing  capacity  in  sand  but  to 
a  lesser  extent  compared  with  other  theories  with  the  exception  of  Balia’s  theory.  For  example, 
Selig  and  McKee41  measured  an  average  base  resistance  of  22.2  lb  in’2  using  3  in  wide  rect¬ 
angular  footings  with  an  aspect  ratio  varying  from  5:1  to  7:1.  The  density  of  sand  was 
0.065  lb  in’3  and  0  was  ‘'•round'  41°  -  to  quote  the  author  -  and  yield  a  theoretical  base 
resistance  of  19.5  lb  in'2. 
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Comparisons  between  drivrn  or  buried  footing  tests  results  and  thf  tluorrti.  ,>l  \..!ues 
is  difficult  because  the  sand  fulling  into  t tie  cavity  left  by  the  footing  present  indeterminate 
shaft  conditions,  An  indication  of  comparative  trends  may  be  obtained  by  assuming  the  valuer, 
of  kyS  and  ?  to  be  1  and  0  res j>ect ively  in  accordance  with  the  reasoning  stated  in  section 
2.8. 


Leboque39  conducted  driven  footing  tests  in  a  poorly  graded,  round  gunned  sand  at 
three  different  densities  ranging  from  1.125  gut  cm"'1  to  1.57  gm  cm'3  (O.Ooio  lb  in*'  and 
0.0566  ii,  hf3  )  foi  i.lilcli  C  »ai  red  fiou.  22°  to  38''.  As  tin-  properties  of  the  sand  show  a 
marked  similarity  to  the  Leighton  Buzzard  sand  used  in  the  current  series  of  tests,  it  is 
suspected  that  maximum  density  was  not  achieved.  Thus,  compressibility  and  dilation  effects 
may  be  present,  and  only  the  results  for  the  dense  sand  are  considered  (Fig.  4.7.1.).  The 
discrepancy  between  theory  and  practice  is  greater  near  the  surface.  This  may  well  be  due  to 
compressibility,  which  would  yield  a  higher  experimental  bearing  capacity  because  compac¬ 
tion  under  the  footing  increases  the  density  of  the  sand  prior  to  failure.  The  improved 
correlation  with  dep.h  may  be  attributed  to  the  fact  that  the  footing  penetration  required  to 
cause  failure  becomes  less  significant  in  terms  of  the  total  penetration  depth  for  the  larger 
penetration  ratios. 

The  correlation  between  the  theoretical  penetration  resistance  and  the  experimental 
data  obtained  by  Wilis  in  very  compact  Ohio  sand  show  an  entirely  opposite  effect 
(big  4.7.2.).  There  is  close  agreement  up  to  a  penetration  ratio  of  2  and  thereafter  the 
experimental  values  increase  more  rapidly  than'the  predict?!  values.  Inaccuracies  in  the 
measurement  of  0  may  well  explain  this  divergence.  Although  '.lie  penetration  test  was 
conducted  in  sand  at  a  den-ilv  cf  0.067  lb  in"',  the  initial  density  of  the  triaxial  shear  lest 
samples  was  0.059  lb  in'3.  Using  confining  pressures  of  up  to  150  psi,  Wills  obtain'*.-;  0  r  40 3 
whereas  the  author  obtained  0  -  40.5  in  the  same  sand  at  a  density  of  0.065  lb  in'3  (r.f. 
section  3.2.).  The  values  of  M-,  are  extremely  sensitive  to  changes  in  0  when  0  is  large. 

For  example,  a  1°  change  in  0  from  40°  to  41n  increased  lire  value  of  \y  by  approximately 
2o?£.  1  his  increase  would  be  sufficient  to  yield  a  predicted  curve  almost  identical  with  the 
mean  experimental  curve  throughout  its  length. 

With  this  general  tendency  towards  theoretical  underestimation, _  it  is  rather  suprisinp 
to  find  that  the  theory  overestimates  the  test  data  obtained  by  Biarez  3'  in  the  roller  mode!. 

At  a  penetration  ratio  of  2  the  theoretical  value  :  s  29.6  kg  cm’1  as  opposed  ic  an  average 
experimental  value  of  22.5  kg  cm",  l'his  tends  to  indicate  that  the  model  is  suspect. 

Biarez  found  that  the  value  of  the  coefficient  of  earth  pressure  k-/S,  on  the  rough  sices 
of  the  foundation  was  2.6  for  0  -  26°  .  This  is  also  surprisingly  high.  From  an  analysis  of 
skin  friction  measurements  on  piles,  Meyerhof  found  that  kvs  varied  from  0.5  for  loose  sand 
to  1  for  dense  sand.  Similar  results  to  those  of  Meyerhof  are  quoted  by  Vesic  for  long  rect¬ 
angular  foundations. 
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s.  i  I'Kiiiru  Tiiri.iur. i icai.  invfstigations 


Putinr  tin-  iIi-m  of  tho  tliooty  of  penetration  resistance,  ancillary  invest!  pit  ions 

yicl, It'll  i . u -  iiit i' n- si i it|?  results. 

5.1  Th rjl ini ads in ’  |-,nll.icv 


Tin-  bearing  capacity  of  n  foundation  is  frequently  derived  hy  taking  moments  of  all 
llie  soil  lot re«  about  sspsfi..:  c-;.v>  jiieiit  pci;;;.  V.licn  the  curved  boundary  o f  the  /.one  of  shearing 
flow  is  di'sciibrd  by  the  equal  ion  to  the  logarillimic  spiral,  the  calculation  is  simplified  by 
taking  moments  about  the  pole  of  the  spiral  because  the  frictional  forces  acting  on  the.  spiral 
boundary  pass  through  the  pole  of  the  spiral  and  may  be  neglected  in  the  moment  equation.  In 
order  to  overcome  the  problem  I  hat  the  stress  field  composed  of  radii  and  logarithmic  spirals  is 
true  only  in  the  case  of  a  weightless  soil,  a  smaller  failure  zone  with  a  modified  shape  is 
achieved  for  a  soil  with  weight  by  relaxing  the  constraints  on  the  position  of  the  pole.  The 
pole  does  not  necessarily  lie  at  the  corner  of  the  footing,  hut  at  some  point  close  to  it  such 
that  the  resultant  bearing  capacity  factor,  Ny,  is  a  minimum.  As  has  been  described  previously, 
this  view  is  rejected  because  the  resulting  inconsistencies  in  the  calculation  appear  to  be 
greater  than  the  original  problem.  Nevertheless  a  considerable  amount  of  work  was  expended  in 
investigating  the  effect  of  the  wandering  pole. 

If  the  pole  is  restrained  to  move  along  the  boundary  between  the  plane  and  radial  shear 
zones  or  that  line  extended  (line  OBD,  Fig  1.4.1a.),  the  variation  in  the  value  of  the  resultant 
frictional  force  on  the  contact  face,  DC,  is  dependent  on  the  distance,  OB,  defined  as  A,  winch 
defines  die  distance  of  the  pole  from  the  corner  of  the  footing  and  is  taken  as  positive  when  the 
pole  lies  outside  the  failure  zone.  Lotkin53  published  the  general  shape  of  the  curve  satisfy¬ 
ing  ll;e  equation:- 

Ny  =  f  (A  )  . 5.1.1. 

The  particular  curve  for  0  -  30°  is  shown  in  Fig  5.1.1. 

It  is  observed  that  the  resultant  factor  can  become  infinite  when  the  pole  is  located  at 
quite  small  distances  from  the  position  of  Am[n  at  which  Ny  is  a  minimum.  This  is  simply 
explained  by  the  fact  that  when  the  pole  is  located  on  the  line  of  action  of  the  resultant  force, 
the  moment  arm  becomes  zero.  Hence,  the  resultant  force  and  N-y  tend  to  infinity.  Similarly, 
if  the  pole  is  located  immediately  beyond  this  limiting  value  of  A,  the  value  of  Ny  tends  to 
a  negative  value  of  infinity,  even  though  the  failure  geometry  is  entirely  reasonable  and  the 
change  in  shape  is  almost  imperceptible. 

Although  the  increasing  values  of  N-y  for  negative  values  of  A  are  related  to  increasing 
size  of  the  failure  zone,  the  minimum  value  of  Ny  is  usually  located  at  a  positive  value  of  A. 
Hence,  the  minimum  value  of  Ny  is  governed,  to  a  large  extent,  by> 

(a)  The  point  of  application  of  tne  resultant  force  for  soil  weight  on  the  contact 

face  or  wedge  boundary  and  its  direction. 

(b)  The  variation  in  the  extent  of  failure. 

(c)  The  angle  of  inclination  of  the  contact  face  with  respect  to  the  horizontal 

(d)  The  degree  of  constraint  placed  on  the  location  of  the  pole. 
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As  the  point  of  application  of  the  force  is  normally  governed  by  hydrostatics  and  the  force  is 
assumed  to  act  at  an  angle  0  to  the  normal,  the  effect  of  these  parameters  can  be  considered 
constant.  The  location  of  the  pole  does  affect  the  extent  of  failure,  but  the  variation  is  not 
systematic  -  the  changes  in  extent  of  failure  and  in  values  of  Ny  are  considered  in  more  detail 
for  a  fixed  pole  (cf  section  5.2.).  The  most  significant  effect  on  the  bearing  capacity  is 
achieved  by  increasing  the  backward  rake  of  the  contact  face  (wedge  boundary)  or  by  decreas¬ 
ing  the  angle  of  inclination  of  the  line  along  which  the  pole  is  restrained.  Separately  or 
jointly,  these’ changes  serve  to  create  greater  constraint  on  the  possible  positions  of  the  pole 
and,  hence,  increase  the  minimum  value  of  Ny. 

In  perspective,  then,  Glide  originally  devised  the  minimising  procedure  for  retaining 
walls  for  which  the  calculation  was  an  acceptable  way  of  optimising  for  the  interface  proper¬ 
ties  along  the  contact  face.  Terzaghi  increased  the  constraints  on  the  pole  by  considering  a 
backward  raked  contact  face.  Meyerhof  relaxed  the  constraints  again  by  adopting  the 
unrestrained  pole  no  longer  confined  to  a  line,  but  free  to  take  any  position  in  two  dimensional 
space.  Fortuitously,  however,  this  freedom  of  movement  was  largely  inhibited  and  governed  by 
the  same  limits  determined  for  the  restrained  pole  i.e.  the  line  of  action  of  the  resultant  force 
effectively  halving  two  dimensional  space  and  an  excessively  large  failure  geometry.  The 
location  of  the  unrestrained  pole  is  thus  never  far  removed  from  the  location  of  the  restrained 
pole.  Wills,  modifying  Meyerhof’s  theory,  again  restrained  the  pole  to  the  boundary  BD  or  to 
that  line  extended  (Fig  1.5.1.);  here  the  angle  of  inclination  of  the  boundary  to  the  horizontal 
decreases  with  sinkage.  Hence  the  constraint  is  more  severe  and  die  final  answer  is  nearer  to 
that  achieved  in  the  absence  of  any  minimising  procedure. 

The  effect  of  the  minimising  techniques  on  the  value  of  the  surface  bearing  capacity 
factor  N-y  may  be  summarised,  for  0  =  4CP  and  </»  =  45  +  0/2,  as  follows; - 

No  minimising  (corner  pole)  163 

Line  restraint  129 

Unrestrained  119 

At  a  penetration  ratio  of  2.5  the  approximate  values  of  the  Ny  factor  for  the  pole  restrained  to 
a  line  (Wills)  and  the  unrestrained  pole  (Meyerhof)  are  640  and  500  respectively.  The  trend 
indicated  by  these  figures  is  similar  to  that  for  the  surface  case. 

It  has  been  demonstrated  that  different  pole  locations  can  provide  a  large  number  of 
solutions  for  the  resultant  forces  generated  on  the  contact  face  by  the  soil  mass  during  failure. 
This  is  plainly  impossible.  This  absurdity  results  f:om  using  a  method  which  does  not  relate 
the  frictional  forces  on  the  spiral  boundary  to  both  the  soil  properties  and  the  dimensions  of 
the  failure  surface. 

The  apparently  cunning  method,  whereby  knowledge  of  the  frictional  force  on  the  spiral 
boundary  is  not  required,  is  fallacious.  The  resultant  force  or.  the  wedge  boundary  can  be 
calculated  only  if  all  the  forces  are  known,  are  shown  to  be  in  equilibrium  and  are  compatible 
with  the  soil  properties.  Such  knowledge  can  only  be  gained  by  a  process  of  rigorous  analysis 
but  owing  to  the  complexity  of  the  problem  has  not  been  fully  achieved.  A  reasonable 
alternative  for  immediate  applications  is  to  accept  a  plausible  failure  surface  with  some 
rigorous  foundation  and  eliminate  the  dangers  inherent  in  the  minimising  procedures. 
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5.2  The  Effect  of  the  Wedge  Angle 


The  comparative  effect  of  the  wedge  angle,  yp ,  on  the  surface  bearing  capacity  is  more 
apparent  when  minimising  techniques  are  eliminated.  For  the  pole  of  the  spipil  located  at  the 
corner  of  the  footing,  the  factors  Nc,  Nq  and  N-y  are  given  by  the  following  equations:- 

Nc  =  tan  yp  +  cos  {yp  -  0) 
sin  0  cos  yp 


a 


a2*U 


+  sin  0) 


-Cl 


5.2.1. 
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where  a  Q  -  e  ^  ^  ^  tan  ^  and  the  range  of  values  for  the  wedge  angle  is 

0  yp  4  7r/4  +  0/2.  In  order  to  complete  the  comparison,  the  rupture  distance,  f,  measured 
horizontally  from  the  corner  of  the  footing  to  the  point  at  which  the  failure  boundary  joins  the 
soil  surface,  can  be  used  in  dimensionless  form,  f/B,  to  give  the  extent  of  failure:- 


f/B 


a^cos  e/cos  yp 


>.2.4. 


The  ruptuie  distance  and  the  bearing  capacity  factors  are  calculated  for  the  two  limiting 
values  of  the  wedge  angle  over  the  complete  range  of  0  (Figs  5.2.1a.  and  5.2.1b.).  The 
variation  in  the  same  functions  with  respect  ot  yp  is  also  given  for  a  constant  value  of  0 
(Figs  5.2. Id.  and  5.2.1e.).  Although  the  rupture  distances  increase  exponentially  with  respect 
to  yp,  the  curves  for  the  Nc  and  factors  always  exhibit  a  negative  gradient.  The  curve  for 
the  N-y  factor  resembles  a  parabola  reaching  some  minimum  value  with  respect  to  yp  between 
almost  identical  values  of  the  Ny  factor  at  yp  =  0  and  yp  -  v/i  +  0/2.  Thus,  the  value  of 
the  bearing  capacity  factor  is  not  related  to  the  extent  of  failure. 

This  is  completely  unrealistic  and  is  determined  by  the  mathematical  constraints 
imposed  on  the  free  body  of  mobilised  soil  rather  than  to  the  mechanism  of  failure.  Any  system 
of  forces  can  be  resolved  into  the  resultant  horizontal  and  vertical  forces  acting  at  any  point 
in  the  plane  plus  a  couple.  However,  the  force  required  to  balance  the  couple  depends  on  the 
length  of  the  moment  arm  which  in  turn  depends  on  the  line  of  action  of  the  resultant  force 
relative  to  the  corner  of  the  footing.  The  resultant  frictional  forces  act  vertically  when  the 
wedge  angle  assumes  the  minimum-value  of  0,  whereas  they  require  resolution  to  obtain  the 
vertical  component  of  the  frictional  force  for  any  value  of  the  wedge  angle  which  is  greater 
than  the  minimurp.  In  addition,  the  moment  arm  for  the  small  wedge  is  less  than  that  for  the 
large  wedge.  These  combined  characteristics  more  than  counteract  for  the  greater  extent  of  the 
failure  zone  in  the  determination  for  the  larger  values  of  Nc  and  Nq.  The  parabolic  form  of  the 
equation  for  is  probably  due  to  the  cubic  exponential  function  of  the  term  a#,  resulting 
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from  the  monte  rl  product  of  ;nu  and  Ji.slaiuc.  It  it,  puieiy  uinai  s'  ii  lli.it  this  1  ■  ’ nil  oflsets 
the  influence  of  the  wedge  angle  on  the  length  of  the  moment  anil  and  on  the  direction  of  the 
resultant  forces. 

In  view  of  tin  -;c  facts,  the  adoption  of  tin-  correct  values  for  tin:  base  angle  of  the 
wedge  forms  an  essential  ['art  of  the  healing  capacity  calculation  and  must  be  determined  from 
the  boundary  stress  conditions.  C.nnscqucnt  iv  ,  the  determination  of  the  bearing  capacity  Ly 
minimising’ with  respect  to  the  wedge  angle"'  is  of  doubtful  validity.  It  is  also  interesting  to 
note  that  Terzaghi  adopts  a  wedge  angle  which  contravenes  the  stress  boundary  conditions 
to  yield  a  high  value  for  the  Nc  and  \\,  factors  and  yet  went  to  tremendous  lengths  of  minimis¬ 
ing  the  My  factor. 
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6.  CONCLUSIONS 


6.1  A  general  theory  of  penetration  failure  has  been  developed  for  strip  footings  or  founda¬ 
tions  both  lying  on  the  surface  and  buried  up  to  depths  of  not  greater  than  ten  times  their 
width.  This  is  an  improvement  on  previous  theories  in  the  following  ways':- 

(a)  The  stress  field  has  fewer  inconsistencies.  1 

* 

(b)  Only  ore  set  of  simple  assumptions  are  consistently  employed  to  describe  the 
distribution  of  stresses  along  a  failure  boundary. 

(c)  New,  and  quite  different  conclusions  are  made  on  the  effect  of  roughness  of  • 

the  base.  It  is  shown  that  rough  and  smooth  foundations  must  have  exactly  the  same 
penetration  resistance  under  ideal  conditions.  : 

i 

(d)  The  theory  treats  all  the  weight  of  the  soil  within  the  figure  of  rupture  as  a  : 

complete  entity  and  not  partly  as  a  surcharge.  i 

j 

t 

(e)  The  theory  goes  some  way  towards  a  more  realistic  combination  of  cohesion  and  ,  J 

weight,  utilisijg  one  figure  of  rupture.  j 

(f)  The  dimensionless  factors  are  functions  not  only  of  0  and  the  c/TB  and  j 

penetration  ratios,  but  also  of  the  shaft  parameters  i 

t 

(g)  The  computed  results  may  be  applied  more  easily  to  practical  problems 

6.2  The  theory  still  Ins  certain  differences':- 

i 

(a)  Consideration  of  the  effects  of  dilation  had  to  be  sacrificed  for  the  sake  of  f 

simplicity  and  generality  over  a  wide  range  of  parameters. 

(b)  Computations  of  the  dimensionless  factors  proved  impossible  at  the  higher 
penetration  ratios  with  low  values  of  0  because  of  the  excessive  tendency  for  the 
failure  boundary  to  close  into  the  shaft  without  achieving  equilibrium. 

6.3  The  experimental  penetration  resistance  data  support  the  quantitative  theoretical 
values. 

6.4  The  measured  failure  patterns  are  shown  to  be  in  agreement  with  the  predicted 
shapes  only  when  the  effect  changing  soil  properties  and  constraints  are  carefully  evaluated 
throughout  the  duration  of  the  test. 
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7.  IttXOMMKMUTlONS 


The  new  theory  of  penetration  failure  is  given  in  a  form  which  is  amenable  to  further 
extension  in  the  development  of  a  theory  of  slip-.sinkage. 
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APPENDIX  A  -  THEORETICAL  ANALYSIS 

The  complete  theoretical  analysis  for  the  determination  of  bearing  capacity  is  described. 
The  failure  geometry  is  derived  by  considering  the  equilibrium  of  the  plane  shear  zone  and  the 
boundary  conditions.  The  magnitude  and  directions  of  all  the  forces  and  the  moment  arms  are  de¬ 
termined.  By  taking  moments  about  the  pole  of  the  spiral,  the  bearing  capacity  is  evaluated. 

A.  1 .  THE  WEIGHT  AND  WEIGHT- FRICTIONAL  FORCES 

IN  THE  PLANE  SHEAR  ZONE 

The  plane  shear  zone,  BDEG,  is  bounded  by  the  foundation  shaft,  BG,  the  horizontal 
soil  surface,  GE,  and  the  two  failure  planes,  BD  and  DE,  along  which  the  shear  stresses  are 
fully  mobilised  (Fig.  A.  1.  1.  ).  Using  co-ordinate  geometry,  the  weight  of  the  zone,  BDEG,  is 
given  by  the  equation  : 


Wj  =51  L  Pi  (z2  -  Z:  )  +  Pa  (7.3  -  7.1  )  ] 


When  the  failure  boundary  breaks  into  the  horizontal  free  surface,  the  following  relationships 
apply  : 


z,  =  0 
z2  *  z  -  f!  sin 


ri  cos  h  +  z  sin  (c  -  C ) 
cos  (v  -  ') 


A.  1 . 2a. 
A.  1.2b. 
A.  1.2c. 


A.  1.  2d. 


p«  =  r,  cos  C 


A.  l.2e. 


Substituting  for  the  co-ordinates  in  equation  A.  1. 1.  ; 


Wl= 


A.  1.  3. 


where 


Aj  =  z~  sin  (T  -  C)  +  2zr!  cos  «p  -  r{  sin  cos  ^ 


The  weight-frictional  forces  on  the  two  failure  planes  are  a  function  of  the  depth  from  the 
nearest  horizontal  free  surface: 


fm  I  vk 


Integrating  : 


C§0  cTsff 


v  kcv  V 


A.  1.  5. 


2(iwi)B  cos  (•;  -  ;) 


•32 


where : 

Similarly : 


and  after  Integration: 


where : 


a(”1) 


Za 


sin  Z 


dz„ 


n  k 


<PY 


V6 


2  (n  +  1 )  BSt  sin  £ 


z2 


The  weight-frictional  side  thrust  is  : 

Fy.  = 

and  after  integration : 

Fy. 


>kY. 

-  j  z0  d  z0 

COS  i  z, 


v  k  z? 

_ T* _ 

2  cos  6 


A. 1.6. 


A.  1. 7. 


A.1.8. 


A.  1.9. 


A.  2.  THE  MOMENT  ARMS  FOR  THE  WEIGHT  AND  WEIGHT- FRICTIONAL 

FORCES  IN  THE  PLANE  SHEAR  ZONE 


The  moment  arm  for  the  force  due  to  the  weight  of  the  soil  in  the  plane  shear  zone  may 
be  obtained  by  co-ordinate  geometry  : 


d,  =  COS(CP-C)  [  pf  (Zi-Zx)  +  Pi  pa  (Zs  -  Z,)  +  pi  (Zj-  Zi)  ] 


and  substituting  for  the  co-ordinates  from  equation  A.  1.2.  : 


3  Ai  cos(o-  £) 


A.  2. 1 . 


where : 


A*  =  z3  sin2  C)  +  3z'  rj  sin  (--  :)  +  3  zr\  cos'  * 


-  rf  cos  +  sin  C  [cos  '  cos  (:?-£)  +  cos  vj 

.  J 


A. 2. 2. 
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The  moment  arms,  ds  and  d* ,  for  the  weight-frictional  forces,  Fv,  and  F.,  ,  may  be 


found  by  integration; 

j 

r  r 

Vk+v 

2:. 

zo 

z:  d  z;  > 

j  ZE  ’  z. 

cos2  (■--:)  1 

C  B  ) 

y 

cos(^-r) 

1 

5,. 

r!n 

z_  d  zc 

<  ! 

cos  (  +  -■' ) 

✓ 

j  2(n+l )  V->  )  cos  ¥ 

~  - -  ) 

|  (2  n+3)  V,  )  cos (a-') 

where  Vs  is  given  by  equation  A.  1. 6.  and  : 


V,  =  z, (»■*«) 

Similarly : 

,,  .  I  „  2  (n+1 )  Ve  |  cos  •+ 

d4  -  Za  -  -  - 

y  (2  n+3)  VB  j  sin  ' 
where  Vc.  is  given  by  equation  A.  1.8.  and: 

Ve  =  Z,  (»*♦•)  -  Z2^L'3) 

Finally,  the  moment  arm  for  the  side  thrust  may  be  found  by  hydrostatics 


A.  2.  3. 


A. 2. 4. 


A.  2.  5. 


A. 2. 6. 


d.-  -  Z  COS  i 

3 


A.  2.  7. 


A.  3.  THE  COHESIVE,  COHESrVE- FRICTIONAL  AND  ADHESIVE  FORCES 

IN  THE  PLANE  SHEAR  ZONE 

The  cohesive  forces  are  proportional  to  the  relative  lengths  of  the  failure 

boundaries  : 


Cj  =  c  DE  -  c(z  -  r;  sin  '  ) 
cos  (¥-  ') 

C3  =  c  BD  =  cri 

The  cohesive-frictional  force,  Fcl,  is  given  by  the  equation  : 


A.  3. 1. 

A. 3. 2. 
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r  vc  dzo 

"  Zj  cos(v-  C) 


‘VC  V; 


(n+1)  BT  cos  ( v  -  :) 


A. 3. 3. 


where : 


V  \  =  z£(r''’)-Zl(='0 


A.  3.4. 


Similarly : 


(n+1)  B11  sin  C 


A. 3. 5. 


where : 


(»*>). za(»  +  0 


Va  =  Z.'"1  >-Za 


A.3.6. 


The  cohesive-frictional  side  force  depends  on  the  cohesive  component  of  the  coefficient 
of  earth  pressure : 


c  kc,  z 


A.  3.  7. 


The  adhesive  force  on  the  shaft  is  : 


C,  =  c.  z 


A.  3.  8. 


THE  MOMENT  ARMS  FOR  THE  COHESIVE  AND  COHESIVE- 
FRICTIONAL  FORCES  IN  THE  PLANE  SHEAR  ZONE 

The  moment  arm,  d7  ,  of  the  cohesive  force,  Ci  ,  is: 


d7  -  r!  cos  v 


A.  4. 1. 


The  moment  arms  of  the  cohesive  force,  Ca,  and  of  the  adhesive  force,  C, ,  are  zero  as 
the  lines  of  action  pass  through  the  corner  of  the  footing. 

The  moment  arms  for  the  cohesive-frictional  forces  are  found  in  a  similar  manner  co 
their  weight-frictional  counterparts : 


^  I  z>  -  (n+1)  V3  \  cos  <p 
(  (n+2)  Vi  j  cos(r+-0 


A.  4.2. 


where  V,  is  given  by  equation  A.  3.  4,  and  : 


V.i  =  -  z,(: 

A.  4.  3 

Similarly,  the  moment  arm  for  Fc;  is  : 

1  Za  -  (n+1)  Vi  | 

d  9  \ 

COS  V 

A.  4.  4 

\  (n+2)  V2  ! 

sin  C 

where  Vn  is  given  by  equation  A.  3.6.  and: 


V*  =  -  z A. 4.  5. 


The  side  thrust  is  assumed  to  be  uniformly  distributed  and  thus  : 


d:;  =  Z  COS  i 

5 


A. 4. 6. 


A.  5.  SPECIAL  RELATIONSHIP  FOR  '  =  0 

WhenC=0,  the  failure  plane,  BD,  is  horizontal.  Thus,  z2  =  Zj  and  the  forces 
FYa  and  F;S  are  indeterminate.  In  this  case,  the  forces  are  given  by  the  following  equations  : 


F 


>2 


Y  k  r,  z 


A. 5.  1. 


F=2  -  c  kcpc  riC"By  A.  5.2. 

Both  forces  are  uniformly  distributed.  Hence  : 

d4  =  d9  =  1  r,  cos  y  A.  5.  3. 


A.  6.  THE  EQUATIONS  OF  EQUILIBRIUM  FOR  THE 

PLANE  SHEAR  ZONE 

The  general  equations  of  equilibrium  are  presented  in  section  2.3.  By  expanding  and 
rearranging  the  terms,  equations  3.  4.  and  2.  3.  5.  may  be  derived  in  a  form  which  is  more 
amenable  for  solution.  The  equation  for  horizontal  equilibrium  is : 


F  cos  ■'  -  F,  sin(vp-C)  +  Fei  cos  '  -  F02  sin  (V.-  ) 

’i  l  Vi 

=  F  cos  +  Fc)  cos  6  +  Ci  cos  :  -  Ci  sin  (  r  -  C) 

y  a 
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Multiply  each  term  by  B/vB.  Thus: 


where : 


Equation  A.  6. 1 


where : 


When 
required : 


but  when  Q  -  0 : 


Thus  : 


k..^  [  Vs  sin  '  cos  '  -  V6  sin  (  +  -  ?)  cos  (^  -  C)  j 

+  2k..  S,  1 )  [V,  sin  '  cos  C  -  VB  sin(x-')cos (+-C)  j 
V  c 


=  [  Ao  +  Si  A*  I  2  (n+1)  B‘"  sin  C  cos  (c #-  ■[)  A.  6. 1. 

Aa  =  |  k  z2  A. 6. 2. 

1  f 

A 4  =  B  [ri  cos  '  +  k  ,  z  -  (z-n  sin  7)  tan  (cp  -  C)  ]  A.  6.  3. 

S,  =  c /  y  B  A.  6.  4. 

.  may  be  written  in  the  form  : 

kcpY  Rl  +  kec  Ra  *  P<>  =  0  A-6-5- 

R!  =  Vs  sin  '  cos  '  -  Vesin  (v-C)cos  (cc  -  t)  A. 6. 6. 

Ri  =  2  S,  B^*1)  [VjSin  ■:  cos  '•-V2stn(v-')cos('- A. 6. 7. 
Ra  =  [Aa  +Si  A,]2(n^l)  Ba”  sin  C  cos  (tf  -  :  )  A.6.8. 


=  0,  special  forms  of  equations  A.  6. 1.  ,  A. 6. 6.,  A. 6. 7.  and  A.6.8.  are 
k  [  V=  -  2  (n+1)  r,z^3“1'1)  sin  +>cos  +■] 

v?  Y 

+  2  k„  S,  B^'"1)  [  Vi  -  (n+1)  rj  z"  sin  cc  cos  ¥  j 

v  C 

=  [Aa  +  Si  A,]  2  (n+1)  B.Sc  cos  ce  A. 6.  la. 

V,  =  z^*1) 

VE  =  z<a"  +  3) 

Ri=  [z-  2  (n+1)  r!  sin  c?  cos  co  ]  A.  6.6a. 

Ri  =  2  Si  b(“+1)  z  ‘  [z  -  (n+1)  r,  sln+-  cos  •v  ]  A.  6.  7a. 


R3  =  (As  +  S,  A*]  2  (n+1)  B3"  cos  +> 


A.  6.  8a. 


The  equation  for  the  vertical  equilibrium  of  the  plane  shear  zone  may  be  similarly  re¬ 
arranged  : 

F  c.os(e-')-F  sin  '  +  F-aCOS  (s- :)  -  Fci  sin  C 
-  u 

=  W,  +  F  sin  -  +  FCJ  sin  i  *  C,  cos  (V  -  -  C2sin  '  +  C, 

'  •  • 

Multiply  each  term  by  B/\  B  : 

k  [  V6  cos2  (y-  ')  -  V-  sin'  •;  I 
▼  : 

+  2  k  S,  B^*1)  [  Vs  cos' ( v  -  ')-  V>  sin2  C 

-  ..  As  +  As  +  S2  A?  -  2  (n+1 )  B' r‘  sin  ■!*  cos  (cp  -  ~  )  A.  6.  9. 

where : 

.  k  za  tan  & 

As  =  - -  +  -ii -  A.  6. 10. 

2  cos ( v -  :)  2 

A£  =  B  [z  (1  +  k*.  tan  5)  -  2  r,  sin  -  ]  A.  6.11. 

A->  =  Bz  A.  6. 12. 

S2  =  c./'.B  A. 6. 13. 

Expressing  equation  A.  6.  9.  in  the  form  : 

k  R4  +  k .  R5  -  R6  =  0  A.  6. 14. 

c-  >  v  c 

the  equations  for  the  constants  are : 

R*  =  CV6  cos2  (c  -  :)  -  Vs  sin2  C  [  A.  6. 15. 

Rs  =  2  Si  B  1  ^  [Va  cos2  («?  -  C)  -  Vj  sin2  "  j  A.  6. 16. 

Re  =  -As  +  Si  As  +  S3  A  7  j  2  (n+1 )  B2”  sin  [cos(c-C)  A. 6. 17. 

When  •:  =  0  ,  a  special  form  of  equation  A,  6.  9.  is  again  required  : 

k^  ^  ri  cos  i  +  k^^  Si  r,  z'cos-v- 

=  [A=  +  S]  Ae  +  Sa  A->  ]  B2"  A.  6.  Sa. 

For  this  special  case,  the  equations  for  the  constants  are : 
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H*  =  r j  1  ^  cos  i 

A.  6. 15a. 

Rs  =  Sj  B  1 )  r>  z“  cos  cr 

A.  6.16a. 

R*  =  (A,  +  SL  A*  +  S2  A i  j  Bac 

A.  6. 17a. 

A.  7. 

THE  SOLUTION  FOR  kr-  ,  n  AND  C  IN  A  FRICTIONAL 
MATERIAL  WITH  WEIGHT 

In  a  non-coheslve  soil,  the  values  of  ks%  and  n  may  be  obtained  in  terms  of  C  from  the 
equations  of  horizontal  and  vertical  equilibrium.  Selecting  the  relevant  terms  due  to  weight  in 
equations  A. 6.1.  and  A.  6.9.  respectively: 


and : 


2  Aa  (n+1)  B2r'  sin  r  (y  -  f) _ 

Vs  sin  '  cos  '  -  V«sin(<?- C)cos(<t  *  O 

2  As  (n+1)  Ba‘  sin  C  cosjy-  :) 

V6  cos2  (i'- ' )  -  VF  sin“C 


A.  7.1. 


A.  7.2. 


Combining  equations  A.  7.1.  and  A.  7. 2.,  and  expanding  the  terms  Vs  and  V6  according  to 
equations  A.  1.6.  and  A.  1.8.  : 


where : 


Thus : 


A,  =  '  (A,  +  A5 ) 

A<5  =  cos  (c  -  ;)  lAa  cos(v  -  ')  +  AEsin(«  -  C) ' 
As  =  Sin  f  [Aa  ‘  C  +  Ae  cos  c  ‘ 

n  =  log  (Ag  +  A9)  -  log  A„  _  j 
2 (log  Za  -  log  Zz) 


A.  7.3. 

A.  7.4. 
A.  7.  5. 

A. 7. 6. 


Similarly  the  expressions  for  kSv  and  n  may  be  determined  from  equations  A.  6.  la.  and 
A.  6.  9a.  when  Z  -  0: 


2  Aa  (n+1)  B2t  cos  ¥ 


and : 


Thus  : 


Vr  -  2  (n+1)  Tj  )  sin  v  cos  cp 

2  Ar  (rf+1)  B2t  cos 
<CPY  2  (n+1)  r,  z<s~l)  cos2  cp 

z  A« 


2  rt  A. 


-  1 


A.  7.  la. 


A.  7. 2a. 


A.  7.6a. 


The  values  of  n  and  k-  ■,  may  be  determined  for  any  value  of  and  the  input  parameters. 
The  value  of  '  satisfying  all  three  equations  may  be  found  by  substitution  in  the  equation  of  rota¬ 
tional  equilibrium  (equation  2.3.5.  ). 


A.  8.  THE  SOLUTION  FOR  k.  c,  n  AND  '  IN  A  FRICTIONAL-COHESIVE 

MATERIAL  WITHOUT  WEIGHT 

In  a  frictional-cohesive  material  without  weight,  the  foundation  shaft  is  assumed  to  be 
unsupported.  Using  the  same  analytical  method  as  in  section  A.  7.,  equations  A.  6.1.  and  A.  6.  9., 
respectively,  may  be  presented  in  the  form  : 


k 


'•’C 


[r,  cos  '  -  (z  -  ra  sin  ')  tan  (c  -  ')j(n+l)  B  r’  sin  '  cos  (■-  -  '  ) 
V..  sin  '  cos  C  -  V=  sin  (  v  -  ■:)  cos  («  -  C) 


and 


s’C  = 


:  z  -  2  r,  sin  -r  J  (n+1)  Br  sin  cos  (v  -  '  ) 
Vj,  cosr  (v  -  C)  -  V)  sin"  ' 


A.  8.  1. 


A.  8.  2. 


Thus : 


n  =  log  [zg  sin  '  +  rj  cos"  (~  -  C)]-  log  [  cos'  -  ')  ] 

- 1 


log  Za  -  log  z. 


A. 8. 3. 


It  may  be  shown  that  rotational  equilibrium  is  achieved  only  when  the  two  failure  boundaries, 
BD  and  DE,  are  of  equal  length  which  gives  the  relationship: 

z  =  r,  sin  '  +r,  cos  (■¥  -  ■') 

Using  this  equilibrium  condition  in  equation  A.  8.  1.  and  A.  8.  3.  ,  k  =  1  and  n  -  0  for  any 
value  of  '  . 


A.  9.  THE  SOLUTION  FOR  k®  c ,  n  AND  C  IN  A  FRICTIONAL-COHESIVE 

MATERIAL  WITH  WEIGHT 

When  both  weight  and  cohesion  are  present,  the  coefficient  kc,'C  may  be  determined  from 
the  quadratic  equation  A.  6.  9.  in  terms  of  n  and  "  : 

i 

k  =  -  R -  *  (Rj  +  4R,  Re)a  A.  9.1. 

H  C  _ _ 

2  R< 

One  root,  derived  by  subtracting  the  discriminator,  is  valid  when  7  <  0  and  the  other  root  is 
valid  when  r  -  0  . 

The  values  of  n  and  "  may  be  found  by  substitution  in  the  equations  for  the  horizontal 
and  rotational  equilibrium  of  the  plane  shear  zone. 

6  0 


A.  10. 


THE  FORCES  IN  THE  RADIAL  SHEAR  ZONE 


The  area  of  the  spiral  sector,  BCD,  is  determined  first : 

=  f '  1  r?  d  ', 


1 


4  tan 


(«V~  -  rcs  ] 


A.  10. 


Thus,  the  weight  of  the  zone  is  : 


W 


A.  10.: 


In  order  to  calculate  the  moment  of  cohesion  about  the  pole  of  the  spiral,  consider  a 
small  element  of  length  5  1  (Fig.  A.  10. 1. ).  Taking  moments  about  B  : 


But : 


d  Me  =  c  dl  r  cos  V 


dl  =  r  d* 


cos  o 


Integrating : 


dM5  =  c  r  di 


M.  =  •  cr!  d) 


(»Y  -  r/ ) 


For  the  special  case  when 


2  tan  c? 

0,  the  spiral  reduces  to  a  circle: 

^  0  r  0  :j 


Mc  =  c  r‘ 


A.  10.  3. 


A.  10.  la. 
A.  10.3a. 


A.  11.  THE  MOMENT  ARM  FOR  THE  WEIGHT  FORCE  EN  THE  RADIAL 

SHEAR  ZONE 

The  determination  of  the  moment  arm,  d*  ,  for  the  weight  force,  W,  of  the  spiral  sector, 
BCD,  is  not  possible  directly.  The  moment  of  area  for  the  spiral  sector  with  respect  to  the  pole 
of  the  spiral  ts  determined  instead  (Fig.  A.  11. 1.  ).  In  the  limit,  as  a  -  0,  the  area  of  the  spiral 
sector  approximates  to  that  of  a  triangle  and  the  centroid  is  then  approximately  that  of  a  triangle. 
Taking  one  axis,  r. ,  along  rQ  and  the  other  axis,  C  ,  normal  to  It  and  passing  through  the  pole  of 
the  spiral,  the  cartesian  co-ordinates  of  the  centroid  are  (n0 ,  C0)  where: 


n  =  r  2/3 


r  c  os  m-  da 


and : 


C0  =  ^/3  r  cos  *  da 


The  moment  of  area,  ±c  C0,  is  given  by. the  equation: 


M.  =  !.  '  r‘  d-i  .  ^'3  :  r  sin  r  da 


V3  r-2  e  ^  ,  tan  '* 


sin  ■)  da 


Integrating  by  parts,  the  final  expression  is  in  the  form: 

M.  =  v  [r,3  (u  sin  j  -  cos  j  )  +  rQ3  3 


where  : 


u  =  o  tan 


3  (liu2) 

The  moment  of  area,  '.c  r,0 .  may  be  derived  in  a  similar  way  : 

M_  =  V3  rc5  'r  e  3  ,tan  v  cos  i  da 


After  integration : 


=  v  [r.3  (sin  y  +  u  cos  j)-urc3  ’ 


A.  11.  1. 

A.  11.2. 
A.  11.3. 


A. 11. 4. 


Transforming  the  ”  ,  axes  at  the  pole  of  the  spiral  lo  x,  y  axes,  the  following  rela¬ 
tionships  are  derived  (Fig.  A.  11.2.  ): 


xc  cos  .  -  -  (rD+  yc  sin  v  ) 
xc  sin  v  =  ~  y,.  cos  1 

The  moments  of  area,  z,0  x0 ,  about  the  pole  of  the  spiral  can  be  written  ir.  the  form  : 

±c  xc  =  ic.  x_(cosp  ,  sin"  i  ) 

=  0C  cos  v  (x0  cos  ■ )  +  .‘.0sin  v  (xcsin  , ) 

Substituting  for  x0  cos  i  and  x0  sin  v  from  equations  A.  11.5.  and  A.  13.6.  ,  the  following 
expression  for  the  moment  of  area  is  obtained  : 

/  x.  -  M,  sin  1  - 

C  t- 

t.2 


A. 11. 5. 
A. 11. 6. 


M  cos  . 


A.  1 1 . 7. 


A. 11.8 


cb  =  x0  =  M,-  sill  >  -  Mr-  cos  ♦ 


The  moment  arm,  d£ ,  is  derived  in  this  general  form  to  be  applicable  to  bearing  capacity 
analyses  in  which  the  pole  of  the  spiral  is  not  necessarily  located  at  the  corner  of  the  footing. 


A.  12.  THE  MOMENT  ARMS  FOR  THE  RESULTANT  FORCES 

It  is  assumed  that  the  resultant  weight-frictional  component  of  pressure  is  distributed 
hydrostatically.  Considering  an  inclined  plane  surface  of  unit  width  immersed  in  a  homogeneous 
material,  the  centre  oi  pressure  may  be  determined  by  hydrostatics  (Fig.  A.  12.  1.). 

L  =  k’d 

where  K  is  the  radius  of  gyration  and  is  given  by  the  equation : 

k?  -  !2  +  d? 


But,  tt  can  be  Ghown  that : 


d  -  z,  +  Zj 


1  •-  z*  -  z, 


Combining  these  relationships,  the  distance,  L,  of  the  centre  of  pressure  from  the  surface  along 
an  inclined  plane  is  : 


L  =  2  (Zj2  +  zaa  +  z,  z«) 
3  (z,  +  z£)  sin  C. 


A.  12.  1. 


The  moment  arm  for  the  resultant  weight-frictional  force  on  the  wedge  boundary  may  be 
found  by  substituting  the  values  : 


z£  =  z  +  V  B  tan  ; 


in  equation  A.  12.  1.  This  yields  the  equation: 


dr  =  B(3z  cos  v  +  B  sin  v)  cos  v 
3(4z  cos  ,  +Bsin  .  )  cos  i 


A.  12.2. 


As  the  resultant  cohesive-frictional  force  is  assumed  to  act  at  the  midpoint  of  the  wedge 


boundary,  the  moment  is: 


dn  =  B  cos  'V 


A.  12.  3. 


4  cos  (45  +  v/2) 


A.  13.  THE  DETERMINATION  OF  THE  RESULTANT  BEARING 

CAPACITY 

The  complete  system  of  forces  and  moment  arms  is  shown  in  Fig.  A.  13.  1.  By  taking 
moments  about  the  pole  of  the  spiral,  the  resultant  forces,  Fy  and  Fc ,  are  determined  for  the 
general  case  when  a  side  thrust  is  applied  to  the  wall  of  the  excavation: 


F  =  W,  d,  +  Wd2  -  Fy,  cb  +  Fy,  cb  A.  13.1. 

do 


Ft  -  C,  d?  -  F„i  do  +  Ft,  die  +  Mc  A.  13. 2. 


d„ 


The  total  force  on  the  base  of  the  footing  due  to  the  weight  of  the  soil  mass  at  failure  is 
obtained  from  the  vertical  component  of  the  force,  Fy  ,  on  both  boundaries  of  the  wedge,  less  the 
weight  of  the  soil  wedge;  the  total  force  due  to  cohesion  is  obtained  from  the  vertical  components 
of  the  forces  Fc  and  C3  acting  on  the  wedge  boundaries  (Fig.  A.  13.2.  ): 


F  lb/ln  =  2F^  cos  (;  -  i)  -  rBs  tan  *  +  2FC  cos  (»  -  «p)  -t  2  C3  sin  , 

4 


This  expression  for  the  base  resistance  may  be  written  in  the  form  of  a  two  part  additive  equat.on: 


F  lb/in  *  |  >BS  Nv  +  cBNc 

where : 

=  4  F-,  cos  (v  -  y)  __ 

vBr 

Nc  =  2  [F.  cos  (v  -  ®)  +  C3  sin  l  * 


A.  13.  3. 

A.  13.4. 

A.  13.5. 


c  B 

The  dimensionless  bearing  capacity  factors,  Ny  and  N„ ,  are  a  function  of  c,  fc,  c/vB,  z/B, 
kyf  and  kc, . 
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FOR  THE  CURVED  SUP  LINES  IN  A  RADIAL  FAN. 


(a)  THE  ASYMMETRICAL  FAILURE  FOR  A  SURFACE  FOUNDATION. 

(b)  FAILURE  FOR  A  BROKEN  BACK  RETAINING  WALL. 


FiG.  1.2.5.  THE  INCOMPLETE  SLIP  LINE  FIELD  FOR  SYMMETRICAL 
SOIL  FAILURE  BENEATH  A  BURIED  FOUNDATION. 


(a)  PRANDTL 


(b)  HILL 


FIG.  1.3.1  SOLUTIONS  FOR  THE  STRESS  FIELDS  IN  AN 


IDEALLY  PLASTIC  METAL  DEFORMED  BY  A  FLAT  PUNCH. 


FOR  THE  MAXIMUM  BASE  RESISTANCE  OF  A 


DEEP  PILE  IN  AN  IDEALLY  PLASTIC  CLAY. 


E 


( c )  y  B  *  O 

FIG.  1-4-2.  FAILURE  PATTERNS  UNDER  A  SHALLOW  FOUNDATION 
FOR  c  =Q  AND  0  =  30°  (AFTER  LUNDGREN  AND  MORTENSEN) 


7B 

FIG  1.4.3.  CORRECTION  FACTOR  FOR  THE  BEARING  CAPACITY 
OF  A  SHALLOW  FOUNDATION 


8  max 


1 


s 

I 

r 


FIG.  1.4.4.  THE  FAILURE  PATTERN  BY  THE  THEORIES  OF 
PLASTICITY  AND  ELASTICITY  (AFTER  GORBUNOV - 
POSSADOV) 


BC-  BOUNDARY  OF  THE  ELASTIC  PART  OF  THE 
COMPACTED  CORE. 

BC  -  BOUNDARY  OF  COMPLETE  CORE. 


FIG.  1.4.5.  THE  FAILURE  PATTERN  BY  THE  VARIABLE 
FACTORS  THEORY  (AFTER  HU) 


FIG.  1.4,6.  THE  FAILURE  PATTERN  BY  THE  STRESS  - 
DILATANCY  THEORY  (AFTER  ROWE) 
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FIG.I.5.1.  FOUNDATION  FAILURE  AFTER  MEYERHOF.  FOR 

COHESION.  SURCHARGE  AND  SKIN  FRICTION  (LHS\ 
AND  WEIGHT  (RHSl).~ 


*1 


FIG. 2.1.1.  THE  ORIENTATION  OF  SLIP  PLANES  AT  THE 


STRUCTURE  FOR  KNOWN  INTERFACE  FRICTION. 


FIG-2  12.  THE  STRESS  FIELD  FOR  A  WEDGE-SHAPED 
FOUNDATION  AFTER  M1NTSKOVSKY. 


FIG. 2.1,3.  THE  ORIENTATION  OF  THE  SLIP  PLANES  AT  THE 
STRUCTURE  FOR  KNOWN  INTERFACE  FRICTION  AND 
ADHESION. 


45*0/2^ 


J  !  I 


FIG.  2.2.1.  THE  SOIL  FLOW  PATH  BENEATH  A 
PERFECTLY  SMOOTH  FOOTING  MODIFIED  BY 
POSSIBLE  CONFIGURATION  OF  THE  FIXED  SOIL  BODY. 


FOR  A  DEEP  FOUNDATION  WHEN  PASSIVE 


PRESSURE  ACTS  ALONG  WALL  BG. 
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13.  ABSTRACT 

unified  procedure  is  employed  in  the  determination  of  the  penetration  resistance 
f  continuous  strip  footings  and  foundations  buried  up  to  ten  times  their  width  in 
ompact  soil.  The  quasi-static  analysis  is  based  on  the  soil  failure  theory  of 
lassical  soil  mechanics. 

he  theoretical  data  is  presented  by  means  of  a  two  part  additive  equation,  two 
imengionless  factors  describing  the  effect  of  weight  and  cohesion.  These  factors 
re  functions  not  only  of  the  angle  of  internal  friction  and  the  penetration  ratio  but 
Iso  of  the  shaft  parameters  and  the  dimensionless  group  dependent  on  the  density, 
the  cohesion  and  the  breadth  of  the  footing.  The  factors  are  calculated  for  the  sami 
igure  of  rupture  and  the  practice  of  treating  a  portion  of  weight  .within  the  failure 
oundary  as  a  surcharge  has  been  eliminated.  The  inherent  approximations  and 
resultant  errors  in  the  computation  are  fully  stated. 

xperimental  verification  is  given  both  of  the  penetration  resistance  and  of  the 
ailure  geometry  by  model  testing  in  a  glass  sided  soil  tank.  The  penetration 
esistance  data  in  compact  soil  support  the  quantitative  theoretical  values.  The 
measured  failure  patterns  are  shown  to  be  in  agreement  with  the  predicted  shapes 
■only  when  the  effect  of  changing  soil  properties  and  constraints  are  carefully 
'valuated  throughout  the  duration  of  the  test.  (  ; 
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